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Voorwoord 


Programmeringstechnieken en in het bijzonder de z.g. lineaire pro- 
grammering, die de eenvoudigste vertegenwoordigster van dit ge- 
zelschap vormt, hebben zich in enkele jaren ontwikkeld van theore- 
tisch interessante methoden tot practische werktuigen, die door 
een toenemend aantal bedrijven worden toegepast voor de oplossing 
van een groot aantal zeer uiteenlopende vraagstukken op het terrein 
van productie, voorraadvorming, afzet e.d. De verdere ontwikke- 
ling van deze toepassingen in ons land wordt echter geremd, omdat 
men niet gemakkelijk de kennis kan vergaren die nodig is om met 
de bedoelde technieken te kunnen werken. Eenvoudige uiteerizet- 
tingen over hun aard en werkwijze zijn wel. voorhanden, zodat de 
belangstellende administrateur zich een beeld kan vormen van de 
draagwijdte der vraagstukken die met de programmeringstechnie- 
ken kunnen. worden aangepakt, maar een handleiding die hen in 
staat stelt zich de technieken zó eigen te maken dat ze ze zelfstandig 
kunnen toepassen, ontbreekt. 

Het boekje dat de heer van der Zijpp over deze materie heeft samen- 
gesteld, beoogt deze lacune op te vullen. Zich baserend op een 
grondige kennis van de theorie van en zijn practische ervaringen 
mét de lineaire programmering heeft hij een leerboek geschreven 
voor de toepassing van deze methode op bedrijfseconomische vraag- 
stukken. 

Hij heeft daarbij niet getracht de moeilijkheden uit de weg te gaan. 
Wil men toch de techniek zelfstandig toepassen, dan is het niet 
voldoende enkele rekenschema’s van buiten te leren met het daar- 
aan verbonden risico van onjuiste toepassingen, maar dan is het 
nodig ook de achtergronden te doorzien. 

Indachtig hieraan heeft de schrijver zich dan ook zeer zorgvuldig 
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afgevraagd welke wiskundige hulpmiddelen voor zijn doel onont- 
beerlijk zijn. Deze heeft hij op duidelijke en overzichtelijke wijze be- 
handeld, met weglating van al het overige. 

Zijn boek is daarom geen litteratuur voor diegenen, die met de wis- 
kunde op zeer gespannen voet staan. Het is echter toegankelijk 
voor iedere bedrijfsman, die enige belangstelling voor dit vak bezit 
en die nog enig begrip heeft van de schoolalgebra. Kennis van de 
hogere wiskunde is niet vereist. Enkele onderwerpen uit de hogere 
algebra, die onmisbaar zijn, worden stap voor stap en op bevat- 
telijke wijze ingevoerd. Voorts is waar mogelijk aangegeven welke 
delen bij eerste lezing overgeslagen kunnen worden. 

Naast de theoretische uiteenzettingen is ook een ruime plaats 
toegekend aan de practische toepassingen. Reeds in de inleiding 
worden enkele bedrijfsproblemen ten tonele gevoerd, die de lezer 
een beeld willen geven van de gevarieerdheid der toepassingsmoge- 
lijkheden. In het laatste deel worden verschillende practische ge- 
vallen uitvoerig besproken. In het bijzonder wordt veel aandacht 
besteed aan een concrete uitwerking van een gecombineerde pro- 
duktie- en afzetplanning, waarbij een aantal belangrijke aspecten 
van het bedrijfsgebeuren ter sprake komen, die de lezer ongetwijfeld 
zullen stimuleren het bestudeerde toe te passen op vele verwante 
problemen. 

Ik spreek daarom gaarne de hoop uit dat het boekje de belangstel- 
ling zal verwerven waarop het mi. ten volle recht heeft. 


Heemstede, Pror. P. DE WOLFF 





Inhoud 


II. 


II. 


IV. 


Inleiding 
$ 1 Enkele voorbeelden 
$2 Probleemstelling 


Vectoren en lineaire vergelijkingen 

$ 1 Vectoren en matrices 

$ 2 Determinanten en inverse matrices 
$3 Verzamelingen van vectoren 

$ 4 Lineaire vergelijkingen 


Lineaire programmering als wiskundig probleem 
$ 1 De oorspronkelijke en de herleide vorm 

$ 2 Bruikbare en optimale oplossingen 

$3 Substitutie van vectoren in een basis 

$ 4 De simplexmethode 

$5 Degeneratie 


De toepassing in het bedrijf 

$ 1 Het produktiemodel 

$ 2 Het verkoopmodel 

$3 Verkoop- en produktieplanning 
$4 Transportplanning 

$5 Dualiteit 


Literatuur 





11 
130 
147 
165 


175 
































|. Inleiding 


$ |. Enkele voorbeelden 


Ll. Efficiënt werken is doorgaans alleen mogelijk, wanneer men 
volgens een zorgvuldig opgesteld plan handelt. Dit geldt voor een 
bedrijf als geheel; het geldt ook voor-elke afdeling afzonderlijk. 
In het algemeen zal men er daarom naar streven zoveel mogelijk 
vooraf vast te stellen, welke prestaties moeten worden geleverd, 
welke werkzaamheden hiervoor nodig zijn en wanneer en door wie 
deze werkzaamheden moeten worden uitgevoerd. 

In vele gevallen is de planning betrekkelijk eenvoudig, omdat de te 
leveren prestaties een gegeven vormen en er in de loop van de tijd 
een vaste procedure voor de uitvoering der werkzaamheden is 
ontstaan. Het komt echter ook voor, dat de planning zeer ingewik- 
keld is doordat vele alternatieve doeleinden en werkmethoden denk- 
baar zijn. Als de keuze van de te volgen politiek dan tevens een 
belangrijke invloed op het ondernemingsresultaat uitoefent, kun- 
nen voor de planning een veelheid van analyses en berekeningen 
nodig zijn. 


2. Bij een omvangrijke planning zal men in de eerste plaats de in de 
toekomst als gegeven te beschouwen omstandigheden zo goed moge- 
lijk trachten te ramen. Deze omstandigheden kunnen bijv. zijn de 
vraag naar een bepaald artikel, de concurrentie positie, het aanbod 
van grondstoffen, de arbeidsverhoudingen en de stand van de 
techniek. 

In de tweede plaats zal men nagaan, welke maatregelen binnen het 
raam der gegeven factoren zouden kunnen worden genomen, en 
welke consequenties elk der alternatieve mogelijkheden zou hebben. 
Men kan bijv. nagaan welk assortiment op een bepaalde markt 
kan worden aangeboden, welke prijzen kunnen worden gevraagd, 
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in welke machines men zou kunnen investeren enz. De consequenties 
van elke maatregel zullen in het algemeen talrijk zijn. Ten dele 
zullen deze gevolgen te kwantificeren zijn als kosten en opbrengsten. 
Vrijwel elke maatregel heeft echter ook een aantal niet-kwantifi- 
ceerbare gevolgen welke voor het bedrijf van groot belang kunnen 
zijn. 

De laatste fase in de mien is steeds de vergelijke van de 
alternatieven en het doen vän een keuze uit de verschillende moge- 
lijkheden. Hierbij zal men primair nagaan, welke gedragslijn de 
grootste winst oplevert. Op de beslissing kunnen evenwel ook ver- 
ML schillende niet-meetbare overwegingen van invloed zijn. 


3. Vanzelfsprekend zal een uitgebreide planning vaak verschillende 
p 5 statistische en boekhoudkundige problemen met zich brengen. Zon- 
der de toepassingsmogelijkheden van de wiskunde in het zaken- 
leven te overschatten kan men zeggen, dat het bovendien soms 
u voordeel oplevert de planningvraagstukken als wiskundige pro- 
blemen te formuleren. 
Het laatste kan in de eerste plaats voortvloeien uit het feit, dat het 
aantal te beschouwen alternatieven zeer groot is. Hiervan is bijv. 
sprake als- de te overwegen maatregel bestaat uit het vaststellen 
van een-te verkopen of te produceren hoeveelheid, een verkoop- 
prijs, een bestelhoeveelheid of de omvang van een lening. Aan elk 
van deze grootheden zou men in beginsel een isin aantal waarden. 
kunnen toekennen. 
Een tweede reden om bij de planning van wiskendige technieken 
gebruik te maken, is het bestaan van interrelaties tussen de te be- 
schouwen grootheden. Indien de verschillende opbrengsten- en 
kostencomponenten wederzijds samenhangen, is men soms genood- 
zaakt deze samenhang in vergelijkingen weer te geven, teneinde een 
„ voldoend nauwkeurig inzicht in de consequenties van de overwogen 
maatregelen te verkrijgen. 





kr 
4 


4: Doorgaans kunnen alleen meetbare verschijnselen wiskundig 
worden geformuleerd. In het volgende zal de aaridacht dan ook . 
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voornamelijk worden geconcentreerd op de kwantitatieve aspecten 
van de planning. 

Hiermede wordt geenszins ontkend, dat andere omstandigheden 
van grote betekenis voor de te volgen politiek kunnen zijn. 


5. De juiste wiskundige formulering van planningsvraagstukken 
moet van geval tot geval worden vastgesteld. Door het algemene 
karakter van de wiskunde kunnen de meeste problemen echter in 
een beperkt aantal standaardvormen worden onderverdeeld. Een 
van de meest bekende wiskundige standaardvraagstukken uit de 
planning is het lineaire programmeringsvraagstuk. Hierbij is de 
opgave, het maximum of minimum van:een lineaire functie te 
bepalen onder de voorwaarde, dat de variabelen in deze functie 
niet negatief zijn en bovendien aan een stelsel gelijkheden of onge- 
lijkheden voldoen. 

In het volgende deel van deze paragraaf wordt i in een viertal voor- 
__ beelden toegelicht, hoe de planning tot een dergelijk vraagstuk kan 
leiden. Daarna wordt in $2 een algemene formulering hiervan 
gegeven. 


6. In de eerste plaats wordt een voorbeeld van aa en  pro- 
duktieplanning beschouwd. 

Stel, in een afdeling van een fabriek kunnen vier verschillende 
produkten worden gemaakt, aangegeven door de cijfers 1, 2, 3 en 4. 
Voor de fabricage beschikt men over drie machines A, B en C. 
Het aantal machine-uren, nodig om 1000 eenheden van produkt 1 
te laten. bewerken door machine A wordt voorgesteld door het 
symbool a11 (lees: a één één). Het aantal machine-uren, nodig om 
1000 eenheden van produkt 1 te laten bewerken door machine B 
wordt voorgesteld door aja (lees: a één twee). De van de verschil- 
lende apparaten vereiste uren voor produktie van 1000 eenheden 
der overige produkten worden op overeenkomstige wijze aangegeven 
door een letter a en twee indices. Men verkrijgt dan het volgende 
beeld. 
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Aantal machine-uren vereist per 1000 eenheden produkt 


\ produkt 

. 1 2 3 4 

machine \ 
A au a12 413 ara 
B aa1 423 423 424 
C 431 432 433 434 


In een zeker kwartaal worden de aantallen beschikbare draaiuren 
van de machines, na aftrek van de verwachte stilstand wegens 
onderhoud en reparaties, resp. geschat op 


Capaciteiten 

Machine Beschikbare uren 
A bi 
B ba 
C b3 


Er wordt aangenomen, dat er voldoende arbeiders en andere pro- 
duktiemiddelen aanwezig zijn, om de beschouwde machines te 
laten draaien. De capaciteiten bi, ba en bg vormen dus de ‘bottle 
necks’ voor de produktie van de betrokken afdeling. De variabele 
kosten van grond- en hulpstoffen, energie, arbeid enz. verbonden 
aan de produktie van 1000 eenheden van elk produkt worden 
geraamd op: 


Variabele kosten per 1000 eenheden 


Produkt 2 3 4 
Variabele kosten ki ka R3 ka 


De af fabriek opbrengsten van de vier goederen worden geschat op 


_ Opbrengsten per 1000 eenheden 


Produkt Î 2 3 4 
Opbrengst 01 02 03 Oa 
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Gevraagd wordt, het verkoop- en produktieprogramma zodanig 
vast te stellen, dat de winst maximaal is bij de gegeven capaci- 
teiten, vaste en variabele kosten en opbrengsten. 

Ziet men af van lange termijn overwegingen — later wordt hierop 
verder ingegaan — dan behoeft bij de oplossing van bovenge- 
noemd vraagstuk alleen rekening te worden gehouden met de feite- 
lijke opbrengsten en variabele kosten. De vaste kosten waarmede 
in de boekhouding wordt gerekend (afschrijvingen op gebouwen en 
machines, lonen van het vaste personeel enz.) kunnen buiten be- 
schouwing worden gelaten. Deze beïnvloeden wel de hoogte van 
de winst, maar niet de keuze van het meest winstgevende program- 
ma. Het probleem kan dan als volgt worden geformuleerd: 

De te verkopen en te produceren hoeveelheden van de vier goederen 
worden voorgesteld door «1, x2, x3 en xa (in 1000 eenheden). 

Op grond van de beschikbare capaciteiten moet dan gelden: 


Ari + A12%2 + A13%3 + arava S bi 
Q2iXi + 42242 + AazX3 + Azara S ba 
Qgix1 + agate + agg%3 + agava S b3 


Daar goederenhoeveelheden niet negatief kunnen zijn, moet voorts 
gelden 


%1, %2, X3, Ha =O 


Stelt men het verschil tussen de opbrengsten en de variabele kosten 
van een verkoop- en produktieprogramma voor door W, dan kan 
het gevraagde worden voorgesteld door 


W = (o1 — ki) + (o2 — ka)xa + (og — ka)%3 + (oa — Ra)va = MAX. 


Hiermede is het probleem wiskundig geformuleerd. 


7. Als volgend voorbeeld kan men een farmaceutische fabriek 
beschouwen waarin tabletten worden gefabriceerd welke o.a. 
twee stoffen A en B moeten bevatten. De produktie vindt 
plaats in twee etappen: eerst worden vier verschillende grond- 








6 Inleiding 


stoffen dooreen gemengd, daarna worden uit het mengsel tabletten 


geperst. 
De gehalten A en B per gram van de verschillende grondstoffen 


worden door de volgende symbolen voorgesteld. 


Gehalten der grondstoffen per gram 


\ Grondstof 1 2 3 4 
Kenmerk \ 
Hoeveelheid A in mg ar a12 413 414 
Hoeveelheid B in mg Q21 422 423 424 


Het is mogelijk, dat enkele van de a’s gelijk aan O zijn. Dit zou 
betekenen, dat een grondstof alleen A of alleen B of geen van beide 
bevat. In het laatste geval is sprake van een ‘vulstof’. 

Men wenst nu tabletten te maken, die minimaal bj mg van A en 
minimaal be mg van B bevatten. Aan de overige samenstelling van 
de tabletten worden geen bijzondere eisen gesteld. en ook aan het 
totale gewicht per tablet worden geen verdere pen opge- 
legd. De kwaliteitseisen van de tabletten zijn dus. 


Vereiste kenmerken per tablet 


Hoeveelheid A in mg b1 
Hoeveelheid B in mg be 


De prijzen van de grondstoffen worden voorgesteld door 


Prijzen der grondstoffen 


Grondstof 1 2 8 4 
‘Prijs P1 Pa Ps Pa 


Het probleem is hier, het grondstoffenpakket zodanig te kiezen 
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dat aan de kwaliteitseisen van de tabletten is voldaan en dat de 
grondstofkosten minimaal zijn. 

De wiskundige formulering van dit vraagstuk kan gemakkelijk 
worden gevonden door de per tablet gebruikte hoeveelheden der 
vier grondstoffen voor te stellen door resp. 41, %2, %3, %4. 

De kwaliteitseisen luiden dan 


A11%1 + Grete + Ai3%3 + A1aXa > bi 
A21%1 + a2aX2 + aagr3 + A2a%a > ba 


Voorts moet gelden 
X1, X2, %3, Xa =O 


De eis dat de grondstofkosten (K) minimaal moeten zijn, betekent 
dat voldaan moet zijn aan 


K == pix + Para + Pax3 H Paxa = min. 


8. Als derde voorbeeld wordt het geval beschouwd waarin een 
handelaar twee goederen 1 en 2 wil kopen teneinde deze na zes 
maanden te verkopen. 

In totaal.is f4000 voor de voorraadvorming beschikbaar. De 
inkoopprijzen van de twee artikelen bedragen franco resp. f 100 
en f40 per ton. Er wordt verwacht, dat de verkoopprijzen af 
fabriek na zes maanden zullen zijn opgelopen tot resp. f 120 en 
{65 per ton. 

De handelaar beschikt over een opslagruimte met een totaal vloer- 
oppervlak van 1800 m2. Een ton van artikel 1 neemt 20 m? in 
beslag; een ton van artikel 2 neemt 60 m? in beslag. 

Gevraagd wordt, hoeveel de handelaar van elk van de artikelen 
moet kopen en verkopen opdat zijn winst bij de gegeven prijs- 
verwachtingen, het gegeven vermogen en de gegeven opslagcapa- 
citeit maximaal is. 

Teneinde dit vraagstuk op te lossen stelt men de van de goederen 1 
en 2 gekochte hoeveelheden resp. gelijk aan #1 en x2. Het vraagstuk 








8 Inleiding 


luidt dan: Bepaal x1 en x2 zodanig, dat 


100x1 + 40x9 < 4000 
20x1 + 60x2 < 1800 
X1, X2 > 0 
W == (120 — 100)x1 + (65 — 40)x2 = 
== 20x1 + 25x32 — max. 


9. Tenslotte een eenvoudig transportvraagstuk. Stel, dat een zeker 
produkt vanuit twee havens aan de Amerikaanse Oostkust ver- 
scheept moet worden naar drie West-Europese havens. Door de 
Amerikaanse havens moeten resp. a1 en az ton worden afgeleverd. 
De Europese havens moeten resp. b1, bz en b3 ton ontvangen. Er 
geldt aj + a2 = b1 + be + bs. 

Beide Amerikaanse havens kunnen elke Europese haven beleveren. 
Er is dus een zestal verschillende transporten denkbaar. Deze 
kunnen door het volgende schema worden kb 


Amerika 


LAN 


Europa 


De transportkosten van de aen Amerikaanse naar de 
verschillende Europese havens worden voorgesteld door de vol- 
gende symbolen. 


Transportkosten per ton van Amerika naar Europa 


Amerikaanse Europese | 2 3 
haven haven 


2 kar ko kas 
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Van de indices geeft de eerste de afleveringshaven, de tweede de 
bestemmingshaven aan. 

De vraag is nu, op welke wijze het transport geregeld dient te 
worden opdat de totale transportkosten minimaal zijn. 

Om deze vraag te beantwoorden kan men de te verschepen hoeveel- 
heden voorstellen door 


Kil X12 X13 
X21 X22 X%23 


waarbij de eerste index weer de afleveringshaven, de tweede de 


bestemmingshaven aangeeft. 
Deze variabelen moeten dan in de eerste plaats voldoen aan de 


volgende gelijkheden 


Xi1 + K12 + %13 = 41 
Xa1 + 22  K23 — 42 

Xi + %21 = b1 
12 + 22 =— ba 

X13 + X23 =— bg 


Voorts moeten de uitkomsten positief zijn, dus 
Kil, X12, %13, X21, X22, 23 2 O 


Stelt men de totale transportkosten voor door K, dan moet ten- 
slotte gelden 


K = kuxu + kiaxia + kigxi3 + Ronxen + kooxaa + kagteg = min. 


$ 2. Probleemstelling 


1. De in $1, 6 en 8 beschreven vraagstukken zijn voorbeelden 
van het maximumvraagstuk van de lineaire programmering. De 
algemene vorm hiervan is: 
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Bepaal x1, 42, +" ** ‚Xn zodanig, dat 
Arixi + Arata Herre: + aintn < b1 
Aaii + Aaa Jerrrr + aanAn < ba 
Ami%1 + Amata Herre + AmnXn < bm 
K1, Katten ‚An 20 
C1k1 + Ca%g + eere + CnAn — max. 


De in $ 1, 7 en 9 beschreven vraagstukken zijn voorbeelden van 
het minimum-vraagstuk van de lineaire programmering, luidende: 


Bepaal x1, #2, «+++ ‚Xn zodanig, dat 
Aiki + Arata Herer + Ainkn Z b1 
Azika H Agora Herer + aantn > b2 
Ami + Ama + rr + AmnAn > bm 
Bj ennen ‚An 20 
C1%1 + Cara Jr HH CnXn — min. 


2. Er wordt gesproken van ‘programmeringsvraagstukken’, om- 
dat het doel van de berekeningen het vaststellen van een program- 
ma voor toekomstige activiteiten is. De programmering is ‘Lineair’, 
omdat zowel in de functie cixi + +--*- + entn als in het bijbe- 
horende vergelijkingenstelsel alleen eerste machten van de te 
bepalen variabelen voorkomen. 

Het laatste betekent, dat de lineaire programmering een onderdeel 
is van een meer algemene wiskundige problematiek. 

Er zijn reeds ingewikkelder, bijv. kwadratische, programmerings- 
vraagstukken bestudeerd. Tot nog toe wordt de hierbij ontwikkelde 
theorie echter zeer zelden bij de planning toegepast. Dit vloeit in de 
eerste plaats voort uit gebrek aan statistisch materiaal, waardoor 
men vaak genoodzaakt is van lineaire betrekkingen uit te gaan. 
Daarenboven wordt de toepassing van kwadratische programmering 
en programmering van hogere orde natuurlijk belemmerd door 
de grotere wiskundige gecompliceerdheid hiervan. 
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3. De wiskundige aspecten van de lineaire programmering zijn 
uitvoerig in de theoretische literatuur behandeld. Na de tweede 
wereldoorlog is een stroom van geschriften op dit gebied verschenen 
als resultaat van een intensief wetenschappelijk onderzoek. De 
toepassing van lineaire programmering op praktische vraagstukken 
is daarentegen nog niet zeer verbreid. Slechts enkele grote onder- 
nemingen en sommige overheidsdiensten maken bij hun planning 
van lineaire programmering gebruik. De plannen van vele bedrijven 
worden gebaseerd op een vergelijking van slechts enkele alternatie- 
ven uit een veelheid van mogelijkheden. Vaak wordt niet volledig 
rekening gehouden met de interdependentie van de te bepalen 
variabelen. 

Vanzelfsprekend zal ook een eenvoudige berekening in, de meeste 
gevallen niet ver van de juiste oplossing uitkomen. Door ervaring 
kan men immers vaak al intuitief aangeven in welke richting de 
oplossing moet worden gezocht. De optimale waarden van de be- 
schouwde grootheden kunnen echter in het algemeen alleen door 
een gecompliceerder berekening exact worden vastgesteld. Veelal 
zijn de consequenties van een kleine afwijking van het optimum 
belangrijk genoeg om een uitgebreidere berekening lonend te 
maken, te meer daar deze niet meer werk met zich behoeft te 
brengen als men eenmaal hierop is ingesteld. Waar mogelijk ver- 
dient het daarom meestal aanbeveling, lineaire programmering toe 
te passen. i 


4. In het volgende wordt getracht, een bijdrage te leveren tot 
het overbruggen van de bestaande kloof tussen de theorie en de 
toepassing van de lineaire programmering. Daar het in de praktijk 
veelal alleen op de fundamenten van de wiskundige theorie aan- 
komt, zijn deze zoveel mogelijk belicht. Voorts is relatief veel aan- 
dacht geschonken aan de wijze waarop praktische planningvraag- 
stukken als lineaire programmeringsproblemen kunnen worden 
geformuleerd en opgelost. 

Aan verschillende theoretische vraagstukken is voorbijgegaan. 
Daar het onderwerp zich evenwel moeilijk leent voor een opper- 
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vlakkige kennismaking, is wel gestreefd naar een exacte behandeling 
van de ter sprake gebrachte problematiek. 


5. Het plan van dit boek is eenvoudig: Voorafgaande aan het 
eigenlijke onderwerp worden in II enkele grondslagen van de 
lineaire algebra behandeld. In III worden de belangrijkste delen 
van de wiskundige theorie van de lineaire programmering uiteen- 
gezet. In IV wordt ingegaan op de toepassing van lineaire pro- 
grammering op verkoop- en produktieplanning en op transport- 
planning. 

De lezer die vertrouwd is met vectoren, matrices, determinanten 
e.d., kan II in eerste instantie overslaan en de daarin genoemde 
definities en stellingen bij latere verwijzing hiernaar nalezen. 








IL, Vectoren en lineaire vergelijkingen 


Evenals elk onderdeel van de wiskunde bestaat de lineaire algebra 
uit de definiëring van een aantal begrippen, het invoeren van een 
aantal regels voor het rekenen met deze begrippen en vervolgens 
uit het afleiden van stellingen uit de gegeven begrippen en reken- 
regels. Deze opzet is in dit hoofdstuk gevolgd. In $ 1, 2 en 3 worden 
o.m. definities gegeven van de begrippen ‘vector’, ‘matrix’, ‘deter- 
minant’ en ‘verzameling’, terwijl tevens enige rekenregels hiervoor 
worden geïntroduceerd. In $4 worden op grond hiervan enkele 
stellingen omtrent lineaire vergelijkingen afgeleid. Er wordt alleen 
ingegaan op die onderwerpen, waarvan bij de ontwikkeling van 
de theorie van de lineaire programmering gebruik moet worden 
gemaakt. 


$ 1. Vectoren en matrices 


1. Een aantal naast elkaar geplaatste getallen noemt men een 
rij getallen of wel een rijvector. 

Een aantal onder elkaar geplaatste getallen noemt men een kolom 
getallen of wel een kolomvector. 

Een rechthoekige opstelling van m Xx n getallen noemt men een 
matrix. Een matrix bestaat dus uit m onder elkaar geplaatste rijen 
van » getallen (of m.a.w. uit # naast elkaar geplaatste kolommen. 
van m getallen). 
De getallen waaruit een vector of een matrix bestaat, noemt men 
de elementen van deze vector, resp. van deze matrix. Het is gebrui- 
kelijk om rond vectoren en matrices haken te plaatsen. 
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Voorbeelden: De eerste twee van onderstaande voorbeelden zijn 
rijvectoren; de volgende twee zijn kolomvectoren. 


[B 25 2 —12 [4 O0 0 
3 — 02 

6 2 

3 

5 Se OW 4 0 18 
4 8 9 _—12 12 6 
ä B 1 Hi 


2. Als een matrix uit m rijen en n kolommen bestaat, zegt men, 
dat de matrix van de orde m Xx n is. Men noemt bij het opgeven 
van de orde van een matrix dus eerst het aantal rijen en daarna 
het aantal kolommen. Door de orde te noemen geeft men de vorm 
van een matrix aan. 

Vectoren kan men als bijzondere matrices beschouwen. Een 
rijvector van #» getallen is een matrix van de orde 1 X #; een 
kolomvector van m getallen is een matrix van de orde m X 1. 
Een enkel getal kan als een matrix van de orde 1 X 1 worden 
opgevat. 


Voorbeelden: De orden van de onder 1 genoemde matrices zijn 
achtereenvolgens 1 X 4, 1 Xx 3,3 X1,2xX 1,3 xX4, 2x 2. 


3. In de algebra stelt men getallen vaak voor door letters. Ook 
de elementen van een matrix kan men door letters voorstellen. 
Het is gebruikelijk de elementen van een bepaalde matrix alle 
door dezelfde letter met rechts onder twee verschillende zndices 
aan te duiden. De eerste van de indices geeft de rij aan waarin 
het betrokken getal staat; de tweede index geeft de kolom aan. 
Zo stelt b.v. xj het j-de element op de #-de rij of m.a.w. het d-de 
element in de j-de kolom voor. Men noemt # de rij-index en j de 
kolom-index. 
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Voorbeelden: Een matrix van de orde 3 X 4 kan men voorstellen door 
X11 %12 %13 X14 
X%21 X22 X%23 %24 
X31 X32 %33 34 
Een matrix van een niet nader te specificeren orde m Xx » kan 
men b.v. aanduiden door 


411 G1o; 2E ZE Ain 
421 aag err aan 
Amt Ama" Amn 


In plaats van de x en de a in deze voorbeelden kan men ook andere 
letters kiezen; de indices worden altijd op dezelfde wijze als 
hierboven gebruikt. 


4. Niet alleen de elementen van een matrix kan men door een 
letter aanduiden. Men kan ook de matrix als geheel door een letter 
voorstellen. 

Men pleegt matrices in het algemeen voor te stellen door een vet- 
gedrukte hoofdletter. Wil men aangeven dat de beschouwde matrix 
een kolomvector is, dan stelt men deze voor door een vetgedrukte 
kleine letter. Rijvectoren worden weergegeven door een kleine 
vetgedrukte letter met rechtsboven een accent. 


Voorbeelden: Men schrijft 


A voor Q11 412 °°° Ain 
421 422 *°°rr Aan 
Ami amar ch 
@ voor 411 
421 
Am1 


a' voor [am 412***" ain] 
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5. Onder de getransponeerde matrix van een zekere matrix A ver- 
staat men de matrix welker kolommen gelijk zijn aan de rijen van A 
(of m.a.w. de matrix welker rijen gelijk zijn aan de kolommen 
van A). 

Men kan de getransponeerde van een matrix dus altijd gemakkelijk 
vinden door rijen en kolommen van deze matrix te verwisselen. 
Is A van de orde m x n dan is de getransponeerde van A van de 
orde ” XxX m. 

Voor de zinsnede ‘de getransponeerde van A’ schrijft men A’. 
Het is duidelijk, dat de getransponeerde van A’ gelijk is aan A. 
Het is nu ook duidelijk waarom kolomvectoren en rijvectoren door 
middel van een accent van elkaar worden onderscheiden: een rij- 
vector kan als de getransponeerde van een kolomvector worden 
beschouwd. 


Voorbeelden: De getransponeerde van 
1 3 5 6| is |l 4 1 
4 3 8 2 3 3 3 
1 3 6 7 9.8 6 
6 2 7 


A’ kan worden voorgesteld door 


411 174 Am1 
412 aaa err aAm2 
din Aant" Amn | 


6. Onder de nulmatrix van de orde m X n verstaat men de matrix 
van de orde m X n waarvan alle elementen gelijk aan O zijn. Bij 
elke orde behoort dus een nulmatrix. Elk van de nulmatrices stelt 
men voor door 0. 


Voorbeeld: De nulmatrix van de orde 2 Xx 3 is 


0 0 0 
0 0 O0 
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7. Een matrix wordt vierkant genoemd, wanneer zij evenveel 
rijen als kolommen heeft. De orde van een vierkante matrix is 
dus m X m. 

Onder de hoofddiagonaal van een vierkante matrix verstaat men 
de elementen met gelijke indices. Dus: 


411 
422, 


en, 
ee, 


“ee 
“Amm- 


Voorbeeld: De hoofddiagonaal van de vierkante matrix 


3 —5 8 
9 —7 12 
14 17 25 


bestaat uit de elementen 3, —7 en 25, 


8. Onder de eenheidsmatrix van de orde m X m verstaat men de 
vierkante matrix van de orde mm Xx m waarvan alle elementen van 
de hoofddiagonaal gelijk aan 1 zijn, terwijl de overige elementen 
gelijk aan O zijn. Voor niet vierkante matrices wordt geen eenheids- 
matrix gedefinieerd. 

De eenheidsmatrix van elke orde wordt voorgesteld door 1. 


Voorbeelden: De eenheidsmatrices van de orden 1 X 1, 2 X 2 en 
3 X 3 zijn respectievelijk: 


[1] 


1 0 

0 1} 
1 0 0 
0 1 0 
o o 1} 


9, Onder de eenheidsvectoren van de orde m Xx 1 (of van de orde 
1 X #) verstaat men de m rijvectoren (resp. de „ kolomvectoren), 
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waarvan één element gelijk aan 1 is, terwijl de overige elementen 
gelijk aan O zijn. 


Voorbeelden: De eenheidsvectoren van de orde 6 X 1 zijn: 


EEE 


0 0 
De eenheidsvectoren van de orde m X 1 vormen de kolommen 
van de eenheidsmatrix van de orde m XxX m. De eenheidsvectoren 
van de orde 1 X n vormen de rijen van de eenheidsmatrix van 
de orde n Xx n. 


oo=oO 
oo=-oOo0o 

Ooo 
„-oOoooo 


o 
oo 


10. Een vierkante matrix A van de orde m Xx m wordt symmetrisch 
genoemd, wanneer voor alle 2= 1 +---« men allej=lt"…- m 
geldt 


dij = Uji 
Voorbeeld: Een symmetrische matrix is 


8 4 3 11 
4 7 9 —8 
3 —9 12 he) 
Il —8 lo) l 


Il. Een vierkante matrix A van de orde m x m wordt scheef- 


symmetrisch genoemd, wanneer voor alle #= 1 ++": m en alle 
jelte m geldt 
ay = —4jt 


Het is duidelijk, dat de elementen op de hoofddiagonaal van een 
scheef-symmetrische matrix gelijk aan O moeten zijn. Er geldt 
immers alleen a — — 44 wanneer a — 0. 
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Voorbeeld: Een scheef-symmetrische matrix is 


0 —4 —3 1 
4 0 9 8 
3 —9 0) —5 
Il —8 5 0 


12. Twee matrices A en B worden gelijk genoemd, wanneer zij 
van dezelfde orde zijn en elk element aj; van A gelijk is aan het 
overeenkomstige element 54; van B. 

Zijn A en B gelijk, dan schrijft men A = B; zijn zij miet gelijk, 
dan schrijft men A # B. 


Voorbeeld: Voor een symmetrische matrix A geldt: A’ = A. 


13. Men zegt, dat een matrix A van de orde m Xx n groter is dan 
een matrix B van dezelfde orde, wanneer elk element aj van A 
groter is dan het overeenkomstige element 54 van B is. Is hieraan 
voldaan (dus aj >> bij voor alle # en j) dan schrijft men A > B 
of B <A. 


Voorbeeld: 
8 —_12| > | —1 — 14 
5 27 3 25 
32 1 —6 0 


14. Men zegt, dat een matrix A van de orde mm x n miet kleiner 
is dan een matrix B van dezelfde orde, wanneer elk element a4j 
van A gelijk aan of groter dan het overeenkomstige element bij 
van B is. 

Is hieraan voldaan (dus ai > bij voor alle # en j) dan schrijft 
men A>=B of B <A. 

A >= B betekent dus, dat een van de volgende relaties geldt: 

— hetzij A = B 

— hetzij A > B 

— hetzij sommige elementen van A zijn gelijk aan de overeenkom- 
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stige elementen van B, terwijl de overige elementen van A groter 
zijn dan de overeenkomstige elementen van B. 


Voorbeeld: De eerste van onderstaande matrices is niet gelijk aan 
en niet groter dan de tweede 


7 3 2J2e{f5 38 1 
5 8 0 S 7 —2 


15. Wanneer a en b getallen zijn geldt a >b of a <5 of beide. 
Het is echter denkbaar, dat twee matrices A en B niet voldoen 
aan A > B en ook niet aan A < B. 

Dit is het geval, als sommige elementen van A groter zijn dan de 
overeenkomstige elementen van B terwijl andere elementen van A 
kleiner zijn dan de overeenkomstige van B. 

Is hieraan voldaan dan zijn A en B niet rangschikbaar. 


16. Men noemt A positief wanneer A>0 
niet negatief wanneer A > 0 
negatief wanneer A <0 


niet positief wanneer A <0 


17. Onder de som van twee matrices A en B van dezelfde orde 
verstaat men de matrix die men verkrijgt door bij elk element 
ai van A het overeenkomstige element bj van B op te tellen. 
Men duidt de som aan door A + B. 

Onder het verschil van twee matrices A en B van dezelfde orde 
verstaat men de matrix die men verkrijgt door van elk element 
ay van A het overeenkomstige element bij van B af te trekken. 
Men schrijft A — B. 

Optellen en aftrekken van matrices van verschillende orden wordt 
niet gedefinieerd. 

Het achtereenvolgens optellen en/of aftrekken van matrices (van 
gelijke orde) geschiedt door herhaald toepassen van de hierboven 
gegeven regels. 
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Voorbeelden: 

5 7 + {—-4 8] =| 1 15 

3 —2 7 3 10 1 

3 8 11] — f—8 4 12} = 11 4 —l1 
—l 5 10 18 2 7 —19 7 3 
A) +5) (6) +[-3|=| CO 

3 7 9 8 9 

— 34 2 1 5 — 28 


18. Stelling: Voor elke matrix A van de orde m Xx n geldt 
Ad0=A 


waarbij O0 de nulmatrix van de orde m X »n is. 

Stelling (A + B) + C= A + (B + C) 

Of in woorden: De som A + B + C kan men bepalen door eerst 
de som A + B te berekenen en daarbij C op te tellen, of wel door 
eerst de som B + C te berekenen en daarbij A op te tellen. Men 
noemt de optelling van matrices daarom ‘associatief’. 

De bewijzen van deze stellingen volgen rechtstreeks uit de definities. 


19. Onder het produkt van een getal A en een matrix A verstaat 
men de matrix die men verkrijgt door alle elementen van A met À 
te vermenigvuldigen. Men schrijft AA of Aà. 


Voorbeeld: 
3 25 48 = 75 144 
— 30 4 — 90 12 
8 — 1 24 —33 


Voor —1:A (of A-—1) schrijft men — A. 


Voorbeeld: Met behulp van deze notatie kan men de volgende 
stelling beknopt weergeven: Is A een scheef-symmetrische matrix, 
dan is A’ = —A 
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20. Laat A een matrix van de orde m x n en B een matrix van 
de orde n Xx p zijn. 
Het produkt AB van deze matrices wordt gedefinieerd als de 


matrix CG van de orde m Xx p waarvan de elementen cj; (2 = 1++-** m; 
jm leser p) gelijk zijn aan: 
Ci = aarbij + aiabag vreet + dinÒnj 


Volgens deze definitie vindt men dus het element c/ van het 
produkt C = AB door de elementen van de #-de rij van A te ver- 
menigvuldigen met de overeenkomstige elementen van de j-de 
kolom van B en deze produkten op te tellen. (Daar A van de orde 
m Xn én B van de orde n X p is, heeft een rij van A evenveel 
elementen als een kolom van B.) 

Men kan het produkt AB als volgt volledig uitschrijven: 


Aar Aaa erree Ain Buu Bija eres bip | = 
aai Aaa rrrer aan bai baar: bep 
Amt Ama Amn bni bna: bnp 
n n Nn 
Yi A1sbe1 B AisDs2 vrrrr pe Q1sbsp 
s=1 = s=1 
Nn n 
Xi 42s0s1 D aasbsa +: pd aasbsp 
s=1 s=1 s=1 
n 
A= 5, Amsbs1 3 Amsbsa-*…*+ p> AmsÖsp 
s=1 
Voorbeelden: 


5 8 5 3 —8 

3 —2 7 —2 1 

={5x548X7 5X3H8Xx--2 5X—-848XI 
3X5-2X7 ZX3-2X—-2 3X—-8—2XI1 


Ì 5 4) 


l 13 — 26 
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6 17 33 x) — | 6x +4 17y + 332 
12 —9 24 y 12x — 9y 4 242 
16 7 —3| [z l6x + 7y— 32 
3 —8 7 } [—!4 28 25 17) ={f—103 —68 190 —39 
s 8 8 12 33 4 20 —19 287 200 175 
[9 8 —7) [| 5 16 21 10 —65 404 110 243 


Op grond van de definitie van vermenigvuldiging van matrices 
kan alleen over het produkt AB gesproken worden, wanneer het 
aantal kolommen van A gelijk is aan het aantal rijen van B. Is A 
van de orde m x n en B van de orde n X p dan bestaat dus het 
produkt AB. Het produkt BA bestaat echter niet (aangenomen 
dat niet tevens » — #). 

Slechts wanneer A van de orde m X nen B van de orde n X m is, 
bestaat zowel AB als BA. Echter AB # BA, hetgeen reeds blijkt 
uit het feit, dat het produkt AB van de orde mm X m is, terwijl 
het produkt BA van de orde n X # is. 


21. Stelling: Voor elke matrix A van de orde m x n geldt 

IA = A 

Al = A 
In de eerste betrekking is I de eenheidsmatrix van de orde m Xx m; 
in de tweede betrekking is I de eenheidsmatrix van de orde n X #. 


Stelling: (AB)C — A(BC) 

of in woorden: Het produkt ABC kan men bepalen door eerst 
het produkt AB te berekenen en dit met C te vermenigvuldigen, 
of wel door eerst het produkt BC te berekenen en A hiermede te 
vermenigvuldigen. De volgorde van de matrices in de vermenig- 
vuldiging mag echter niet worden verwisseld. Dit betekent, dat 
de vermenigvuldiging van matrices wel ‘associatief’ maar niet 
‘commutatief’ is. Hierin verschilt de vermenigvuldiging van ma- 
trices van de vermenigvuldiging van getallen. 


Stelling: (A + B)C = AC + BCG. 
De bewijzen van deze stellingen volgen rechtstreeks uit de definities. 
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22. Deling van matrices wordt niet gedefinieerd. Wel behoort 
zoals in $ 2 zal blijken onder zekere voorwaarden bij een vierkante 
matrix A een inverse matrix, aangegeven door A-1, die in dien zin 
als omgekeerde van A kan worden beschouwd, dat 


A-1A = AAT =I. 


23. De invoering van de hiervoor besproken definities maakt 
het mogelijk, stelsels van vergelijkingen beknopt weer te geven. 
Beschouw het stelsel van m vergelijkingen met # onbekenden: 


Bl jp marre + Gin*n = b1 
Amika Herre + AmnXn = Ôm- 


Definieer A als hierboven en stel verder 


X1 = b 
an 


Op grond van 1) de definitie van veemerigvoldieing van matrices 
(in casu A van ne orde mm X n en # van de orde n X 1) en 2) de 
definitie van gelijkheid van matrices (in casu de vectoren Ax en b 
beide van de orde m x 1) geldt dan 


Ax =b 
Op analoge wijze kan men een stelsel van m ongelijkheden en 7 


onbekenden weergeven door Ax < b 
resp. Ax > b 


$2. Determinanten en inverse matrices 
1. De determinant van een vierkante matrix is een getal, dat op 


een hierna te bespreken wijze uit deze matrix kan worden afgeleid. 
Voor niet-vierkante matrices wordt geen determinant gedefinieerd. 
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In plaats van de zinsnede ‘de determinant van A’ schrijft men 
|A| of wel 


In plaats van haken plaatst men dus rechte strepen aan weerszijden 
van het schema van elementen. 


2. Beschouw A van de orde m x m. Uit alle rijen van A kan 
men een element kiezen, zodanig dat de m gekozen elementen alle 
in verschillende kolommen liggen (of m.a.w. uit alle kolommen 
van A kan men een element kiezen zodanig dat de m gekozen 
elementen alle in verschillende rijen liggen). Men verkrijgt dan m 
elementen die alle uit verschillende rijen en uit verschillende 
kolommen afkomstig zijn. Zowel de rij-indices als de kolom- 
indices van deze elementen hebben de waarden 1 --:-- m. 

Het produkt van m aldus bepaalde elementen noemt men de 
absolute waarde van een term van |Al. 


Voorbeeld: 


411 412 413 G14 415 
421 422 423 424 425 
431 432 433 434 435 
441 442 443 @44 445 
451 452 453 454 455 


De absolute waarden van enkele termen van de determinant zijn 
hieronder weergegeven: 


412. 424.431. 445.453 
441-412.423. 454.435 
433.422. 454.445. A11 
442.423. 455.414. 431 


Vanzelfsprekend blijft elk van deze produkten gelijk, wanneer 
men de volgorde van de elementen daarin verandert. De absolute 
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waarde van een term is dus onafhankelijk van de volgorde waarin 
de verschillende elementen worden gekozen. 


3. Is A van de orde m X m, dan kan men op m! wijzen * een 
term van |A| bepalen. 

Dit is als volgt in te zien: Uit de eerst beschouwde rij (kolom) 
kan men op m verschillende manieren een element kiezen. Heeft 
men een element bepaald, dan kan men op (m — 1) wijzen een 
element uit de als tweede beschouwde rij (kolom) kiezen. Er zijn 
dus m(m — 1) mogelijkheden om de eerste twee elementen te be- 
palen. Bij elke keuze van het eerste en het tweede element kan 
men het derde op (m — 2) manieren kiezen. De eerste drie ele- 
menten kan men dus langs m(m — 1)(m — 2) wegen verkrijgen. 
Op grond van deze redenering kan men aantonen, dat een term 
van |A|, bestaande uit m elementen, op m! verschillende manieren 
kan worden bepaald. 


4, Onder inversie verstaat men het verschijnsel, dat twee getallen 
in een reeks niet naar opklimmende grootte zijn gerangschikt. In 
een reeks komt dus een of een aantal inversies voor, wanneer een 
getal links van één of van een aantal kleinere staat. 


Voorbeelden: In de reeks 1, 2, 3, 4 is het aantal inversies nihil; in 
3, 2, 3, 4 is het aantal inversies één; in 5, 2, 3, 4 is het aantal in- 
versies drie; in 5, 4, 2, 3 is het aantal inversies vijf. 


5. Beschouw de absolute waarde van een term van |A|. Rang- 
schik de elementen in deze term in een willekeurige volgorde. 

Men kent nu aan de term Aet positieve teken toe, wanneer de som 
van het aantal inversies in de reeks rij-indices en het aantal in- 
versies in de reeks kolom-indices even is. Men kent aan de term. 
het negatieve teken toe, wanneer de som van het aantal inversies in 
de rij-indices en het aantal inversies in de kolom-indices oneven is. 


* ml (spreek uit m-faculteit) is een beknopte schrijfwijze voor 
° mlm — Ilm — 2) «--e 43.21, 
Bijvoorbeeld: 4! = 4.3.2.1 = 24; 71 = 7.6.5.4.3.2.1 —= 5040. 
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Voorbeelden: De sommen van het aantal inversies in de rij-indices 
en het aantal inversies in de kolom-indices is hieronder weergegeven 
voor de in 2 gegeven voorbeelden. 


OH4= 4 
4AHO0= 4 
oH6= ll 
645 = 11 


De eerste twee termen krijgen dus het even teken; de laatste twee 
termen zijn oneven. 


6. Stelling: De bepaling van het teken van een term van |A 
volgens 5 is onafhankelijk van de volgorde waarin de elementen 
zijn geplaatst. 


Bewijs: 

a) Eerst wordt de volgende hulpstelling bewezen: Als men twee 
willekeurige getallen in een reeks verwisselt, verandert het aantal 
inversies van even tot oneven, resp. van oneven tot even. 

De juistheid van deze stelling is evident, wanneer de twee verwis- 
selde getallen naast elkaar staan. Is immers het getal dat oor- 
spronkelijk voorop stond, kleiner dan het tweede, dan ontstaat 
door de verwisseling een inversie meer in de reeks; is het oorspron- 
kelijk voorop staande getal groter dan het tweede, dan wordt door 
de verwisseling een inversie teniet gedaan. In beide gevallen 
verandert het aantal inversies met 1. Dit geldt onafhankelijk 
van de waarden der getallen die voor en na de twee beschouwde 
staan. 

Stel nu, men verwisselt de getallen a en b waartussen k getallen 
staan; a staat links van b. De verwisseling van a en b kan plaats 
vinden door eerst a telkens een plaats te laten opschuiven. Bij elke 
verwisseling van a met een daarop volgend getal, verandert het 
aantal inversies met 1. Wanneer men tenslotte a met b verwisselt, 
is het aantal inversies k + 1 maal met 1 veranderd. Vervolgens 
verschuift men b naar de oorspronkelijke plaats van a. Wanneer 
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b op de plaats van a staat, hebben dan nog eens k verwisselingen 
plaats gehad en is het aantal inversies k maal met 1 veranderd. 
In totaal verandert het aantal inversies derhalve (2k + 1) maal 
met 1. Ongeacht het feit of k even dan wel oneven is verandert het 
aantal inversies in de beschouwde reeks dus een oneven aantal 
malen met 1. Dit betekent dat het aantal inversies van even tot 
oneven, resp. van oneven tot even wordt. 


b) Verwisselt men nu twee elementen in een term van |Al, dan 
veranderen 
1) het aantal inversies in de rij-indices 

— hetzij van even tot oneven 

— hetzij van oneven tot even 
2) het aantal inversies in de kolom-indices 

— hetzij van even tot oneven 

— hetzij van oneven tot even. 
Bij elke combinatie van deze mogelijkheden blijkt de som van 
het aantal inversies in de rij-indices en het aantal inversies in de 
kolom-indices even dan wel oneven te blijven. 


7. De determinant van een matrix A van de orde mm X m kan 
nu als volgt gedefinieerd worden: 

|Al is de som van alle m! volgens 2 gedefinieerde termen, waarvan 
het teken volgens 5 is bepaald. 

Met deze definitie ligt het begrip determinant vast. Het is echter 
mogelijk de determinant van een gegeven matrix op eenvoudiger 
wijze te berekenen dan met behulp van de besproken definitie. 
Dit zal in 11 worden behandeld. Hier kan reeds worden opgemerkt, 
dat de determinant van een matrix van de orde 1 X 1 volgens 
de definitie gelijk is aan het element waaruit deze matrix bestaat. 
Dus: lara = 41. 


8. Elke matrix, gevormd door een aantal rijen (kolommen) of 
delen van rijen (kolommen) van A, noemt men een deelmatrix 
van A. 
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Onder de minor van het element ais van A van de orde m Xx n 
verstaat men de deelmatrix van de orde (m — 1) X (n — 1) 
die uit A ontstaat door de j-de rij en de j-de kolom weg te 
laten. 

Voor de zinsnede ‘de minor van a’ schrijft men Aij. 


Voorbeeld: De minoren van a11 en a23 van A van de orde 3 Xx 4 zijn: 


A23 = |a11 a12 G14 
aa1 


432 434 


A10 = |a22 423 aza 
432 433 434 





9, Is A van de orde m X m, dan is Aj eveneens vierkant en wel 
van de orde (m — 1) X (m — 1). De determinant van de minor 


Aij wordt voorgesteld door det. Aj. 
Onder de cofactor van het element aj van de vierkante matrix A 


verstaat men 
(— 14 x det Aj 


Deze cofactor wordt voorgesteld door |Ayy|. 


10. De invoering van cofactoren maakt het mogelijk, de bereke- 
ning van een determinant op beknopte wijze te beschrijven. Er 
geldt namelijk: 


Stelling: Men vindt |A] door de elementen van een willekeurige 
rij of kolom van A te vermenigvuldigen met hun cofactoren en 
deze produkten op te tellen. 

In symbolen 


1) |A = anlAal + as |Asa) +°°°…: + aam |Aam| = D ass |Assl 
s=1 


voor elke j= 1 +…+": m. 
2) |A| = auylAijl + aaslAaj| +-:°"- + Ams |Amjl = Di aap |Asjl 
s=1 
voor elke j= | ++." Mm. 


Men kan deze berekening noemen ‘het berekenen van |A] volgens 
de elementen van de i-de rij (resp. j-de kolom) van A’. 
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Bewijs: Alleen het eerste deel van de stelling wordt bewezen; het 
bewijs van het tweede deel kan geheel analoog hieraan worden 


gegeven. 
De m! termen van |A] kunnen in m groepen worden verdeeld 
op grond van het criterium of a;1, az2, «+*** dan wel a;m erin 


voorkomt. De som van de termen die aj bevatten, kan worden 
geschreven als ajj X de som van (m — 1)! produkten van m — 1 
elementen. De absolute waarden van de laatste produkten zijn gelijk 
aan de absolute waarden van de termen van det Ajj. Het teken 
van een term van det Aj; kan volgens 5 worden bepaald. Plaatst 
men aj; voor één van de (m — 1) termen van det Ay, dan neemt 
het aantal inversies in de rij-indices met £ — 1, het aantal inversies 
in de kolom-indices met j — 1 toe. Volgens 6 geldt dit ook als 
men ai een andere plaats dan de eerste in de reeks geeft. Het totaal 
aantal inversies in de rij- en de kolom-indices neemt dus met 
1 +i— 2 toe. Als nu j +7 — 2 (dus j +1) even is, hebben de 
termen van |A] die a; bevatten, hetzelfde teken als de bijbehorende 
termen van det Ass; als j + 2 — 2 (dus j + 7) oneven is, hebben de 
termen van |A| die aj bevatten het tegengestelde teken van de 
bijbehorende termen van det Ajs. Daar |Ajy| gedefinieerd is als 
(—1) det Aj kan men het laatste weergeven in de vorm van 
het gestelde. 


11. Met behulp van voorgaande stelling kan men de berekening 
van de determinant van een matrix van de orde m x m her- 
leiden tot de berekening van m determinanten van matrices 
van de orde (m — 1) X (mm — 1). De berekening van de deter- 
minant van een matrix van de orde (m — 1) X (m — 1) kan 
men herleiden tot de berekening van (m — 1) determinanten 
van matrices van de orde (m — 2) X (mm — 2). Uiteindelijk kan 
de determinant van elke matrix dus worden geschreven als 
een veelterm van determinanten van matrices van de orde 
1 x 1, dus als een veelterm van determinanten van afzonder- 
lijke elementen. De laatste zijn, zoals in 7 werd opgemerkt, gelijk 
aan deze elementen zelf. 
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Berekent men de determinanten bijv. volgens de elementen van 
de eerste rij (zie 10), dan verkrijgt men de volgende berekeningen: 








| 411 | = d11 

411 A12 | S 411422 — 412421 

421 422 

411 412 G13l= Gril 422 423|— A12) A21 G23| A13) 421 A22 
421 @22 423 432 433 431 433 431 432 




















431 432 433 
== U11422433 — 11423432 — 4124210433 + 
+ 412423431 + 413421432 — 4134224310 


Met matrices van hogere orden gaat men op analoge wijze te werk. 


Voorbeelden: 
6 8 3|—= —136 8 O0 3 —10|== 13479 
Z2 1 | —7 5 4 
5 2 4 12 2 6 25 
2 7 O0 —9 


12. Stelling: De determinant van een matrix waarvan twee rijen 
of twee kolommen aan elkaar gelijk zijn, is gelijk aan 0. 


Bewijs: De determinant van de matrix die men uit een willekeurige 
vierkante matrix A verkrijgt door twee rijen: (kolommen) van A 
te verwisselen, is gelijk aan —|A|. De verwisseling betekent im- 
mers, dat in elke term van |A| twee rij-indices (kolom-indices) 
van plaats verwisselen, terwijl de kolom-indices (rij-indices) niet 
van plaats veranderen. Indien tussen de twee beschouwde rijen 
(kolommen) k andere liggen, verandert het aantal inversies in de 
indices bij stap voor stap verschuiven van de elementen 2k + 1 | 
maal met 1. Geeft men de nieuwe matrix aan door B, dan geldt | 
derhalve |B] — — [A]. | 
Als de twee verwisselde rijen (kolommen) gelijk zijn, is B — A, 

dus |B| = |A. Derhalve is —|A| = |A|. Dit is alleen mogelijk, 

indien |A| =O. 
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13. Men kan nog andere eigenschappen van determinanten dan 
de in 12 behandelde bewijzen. Zo kan men bijv. aantonen, dat 
de determinant van teken verandert als men twee rijen (kolommen) 
van A verwisselt, dat de determinant 4 x groter wordt als men 
een rij (kolom) van A met 4 vermenigvuldigt, en dat de determi- 
nant gelijk blijft als men A x een rij (kolom) van A bij een andere 
rij (kolom) optelt. Door deze eigenschappen toe te passen kan men 
de berekening van determinanten vaak bekorten. Daar hiervan 
echter in het volgende geen gebruik behoeft te worden gemaakt, 
wordt hierop niet verder ingegaan. 


14. De inverse van A van de orde m x m, waarvan de determi- 
nant ongelijk aan O is, wordt gedefinieerd als de matrix die men 
uit A verkrijgt door: 

— de elementen 44; te vervangen door hun cofactoren |A4y| 

— deze cofactoren te delen door |Al 

— de aldus verkregen matrix te transponeren. 

Stelt men de inverse van A voor door A-1, dan geldt dus: 











Al=(lAul |Anl Am] 
LAI |AI [Al 
lul lAal  |Ame 

|A| |Al |Al 


|A 1m! |Aam| ee (Am! 
|Al LAI [Al 








Is |A] = O dan heeft A geen inverse. 
De betekenis van bovenstaande definitie blijkt uit de in 15 behan- 
delde eigenschap van de inverse matrix. 


15. Stelling: 
AAI = AlA=lI. 











Determinanten en inverse matrices 33 


Bewijs: Volledig uitgeschreven is 














AAT = (ay «---- Aim lAul Amal 
EE U B HK Om wr |A] |A| 
Ami ns 
Aum] lAmml 
|A] [Al 
1 mm Me 
En lak} zen sl Ams 
|A| 2e | f | zn | | 
5 Ams |Ars! een Di Ams |Amsl 


Volgens 10 zijn alle elementen in de hoofddiagonaal van de laatste 
matrix gelijk aan |A. 

Elk van de overige elementen Deer Gisl Agel (jl m, 
i #j) kan men eveneens als een determinant beschouwen. De 
matrix, waarvan deze uitdrukking de determinant vormt, kan 
men direct uit A afleiden. De determinant van de betrokken matrix 
moet ook te bepalen zijn als Ym, ajslAzsl. Dit betekent, dat de 
i-de en de j-de rij aan elkaar gelijk zijn. Volgens 12 is dus 
Die-1 Giel Agsl =O voor d #j. 

Bovenstaande uitdrukking wordt derhalve: 


Hiermede is het eerste deel van het gestelde bewezen. Het bewijs 
van A-tA = I kan op analoge wijze worden geleverd. 


16. Door de definitie van inverse matrices is het mogelijk, de 
oplossing van een stelsel van m vergelijkingen en m onbekenden 
in het kort weer te geven. 
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Uit Ax =b 

volgt A-lAx = Alb mits althans |A| #0. 

Dus: Is = x = Alb. 


Als |A| = O is het vergelijkingenstelsel onoplosbaar. Hierop wordt 
in $4 nog teruggekomen. 


$3. Verzamelingen van vectoren 


1. In de wiskunde wordt elke groep objecten een verzameling ge- 
noemd. Hier worden verzamelingen van vectoren beschouwd. 
Een verzameling wordt in het algemeen door een hoofdletter voor- 
gesteld. Om aan te geven, dat de vector a tot de verzameling V 
behoort, schrijft men ac V. 


2. Het totaal van alle vectoren met m elementen stelt men voor 
door Rm. 

Volgens deze schrijfwijze geeft R de verzameling van alle reële 
getallen van — oo tot + oo aan. Deze verzameling kan op de 
bekende wijze worden voorgesteld door de punten van een rechte 
lijn. De verzameling van alle vectoren bestaande uit twee elemen- 
ten, R2, kan men voorstellen door de punten van een plat vlak, 
waarin twee loodrecht op elkaar staande assen zijn getrokken. 
R3 kan worden voorgesteld door de punten van de drie-dimensio- 
nale ruimte waarin drie loodrecht op elkaar staande assen zijn 
getrokken. Indien m > 3 kan R” niet ruimtelijk worden voor- 
gesteld. Toch blijft men in overdrachtelijke zin ook bij m > 3 
van Rm spreken als van de m-dimensionale ruimte. 


3. In overeenstemming met het onder 2 genoemde spraakgebruik 
noemt men een vector van de orde m X 1 of 1 X m veelal een 
punt van R” of wel een vector in R”. Om aan te geven dat a uit m 
elementen bestaat kan men schrijven ae R®, Een verzameling 
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vectoren van de orde m X l of 1 X m noemt men een verzame- 
ling in Rm, 


4. Stel, ad, a, ----- ‚a zijn » vectoren in R® en u1, ua, ***** 5 
Hn Zijn n reële getallen. In deze paragraaf wordt in het midden 
gelaten of van rijvectoren dan wel van kolomvectoren sprake is. 
Men noemt dan de vector b die gelijk is aan 


n 
b—= YX uia® 
i=1l 
een lineaire combinatie van ad -..-. a), 
De getallen gi: kn noemt men in dit verband gewichten. 


Men spreekt van een niet-negatieve lineaire combinatie wanneer 
Har: Un > 0. 


tweede 
element 


15 


10 





5 10 15 eerste element 





36 Vectoren en lineaire vergelijkingen 


Voorbeeld: In de tekening op de vorige bladzijde zijn 


218 


in het platte vlak voorgesteld door pijlen, lopend van de oorsprong 
naar resp. de punten (4, 1) en (2,3). De positieve lineaire com- 
binatie 


2f4l +3[2| = {4 

1 5 17 
is voorgesteld door de pijl van de oorsprong naar het punt 
(14, 17). Met behulp van eenvoudige analytische meetkunde 
kan men bewijzen, dat de derde pijl een diagonaal is van het 


parallelogram waarvan twee zijden worden gevormd door 2 x de 
eerste en 3 X de tweede pijl. 


5. Men noemt b= 5%, MA} een convexe lineaire combinatie 


van ad …---. a in Rm, indien de gewichten voldoen aan 
n 
zm! 
i=1 
Rd, kn > 0 


Bij een convexe lineaire combinatie van twee vectoren is 
ui + m2 = 1. Meestal stelt men u1 = u en u2 = | — u. Er geldt dan 


b=(l — wad + ua O<Sus<l 
Voorbeeld: Een convexe lineaire combinatie van 
É) * É) 
6 3 
is bijvoorbeeld 


4 Hi +$ [5] a K (zie bladzijde 37). 








Verzamelingen van vectoren 37 


Met behulp van analytische meetkunde kan men bewijzen, dat 
een convexe lineaire combinatie van twee vectoren in R? steeds 
ligt op de lijn door deze twee punten. 


tweede 
element 


10 





5 10 eerste element 


6. Onder het segment tussen twee punten ad en a@ in R” ver- 
staat men de verzameling van alle mogelijke convexe lineaire 
combinaties van ad en a@®. Het segment tussen ad en a@ is 
dus m.a.w. de verzameling van alle punten 


x= (l —u)ad + u4a® voor OSu Sl. 


Voorbeeld: Het segment tussen twee punten in R? kan worden 
voorgesteld door de punten op de verbindingslijn tussen de twee 
vectoren. Voor u =O verkrijgt men # =a@®d. Voor u = 1 ver- 
krijgt men « = a@. Hiertussen zijn alle combinaties O < u < | 
mogelijk. In de tekening op de volgende bladzijde is gesteld 


ad =f4| en a =[f6 
8 2 


De waarden van « zijn getekend voor u=t, u=tenumu=fÂ. 
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tweede 
element 


10 







alt (44 — 0) 


a@ (py. = 1) 


0 ed 10 eerste element 


7. Een verzameling van vectoren V wordt convex genoemd, wan- 
neer elke convexe lineaire combinatie van elk tweetal vectoren 
van V eveneens tot V behoort. 


Voorbeelden: Het eenvoudigste voorbeeld van een convexe ver- 
zameling is het segment tussen twee punten a® en a®%. Zijn 
immers bD en b® twee punten van het segment, en wel 
bD = (1 — mad + wa en b = (1 — ua)aD + gea, dan 
is een convexe lineaire combinatie c = (1 — A)bW + 4b® even- 
eens een punt van het segment, want: 


e= (1 —A{(L — 2D + 12) + AI — u) + 2} 
= (1 —A(L — u) HAL — JA + (1 — Aan + Aja 


Het is gemakkelijk in te zien, dat de coëfficiënten van ad en a@ 
in de laatste betrekking niet negatief zijn indien wi, ua en ÀÂ > 0, 
en dat de som ervan gelijk aan 1 is. 

In onderstaande tekening zijn vier verzamelingen in R?2 weer- 
gegeven. Hiervan zijn I en II convex; III en IV zijn niet convex. 
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Ì tweede T 
element 


eerste 
element 





8. Stelling: De verzameling van alle convexe lineaire combinaties 
van een eindig of oneindig aantal gegeven punten is convex. 
Bewijs: Stel de gegeven punten voor door a®. De index # kan 
een willekeurig aantal waarden aannemen. 

Beschouw nu twee convexe lineaire combinaties b® en b@® van a@: 


bD = 5 mja Smul, wi 20 
ï ï 


b®) — XE Aja Eh = 1, u 20 


KA 
Een willekeurige convexe lineaire combinatie van bb en b@ is 
c=(l —tjbD + tb O<St<I 
Hiervoor kan men schrijven: 


e= (ll —t) XY ma) Ht 9 ad 
ï ï 


=E(l — Bs + Aja. 
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Het blijkt dus, dat c als een lineaire combinatie van de vectoren 
a® kan worden geschreven. Bovendien blijkt deze combinatie 
echter een convexe te zijn. In de eerste plaats zijn de coëfficiënten 
positief, aangezien zowel (1 — £)u4, t als À; > O. In de tweede plaats 
is de som ervan gelijk aan 1: 


DU Du HAJ = (ll) Du HEA =l. 
i ï i 
Hiermede is bewezen, dat ook ce V. 


9, De punten van een convexe verzameling V kunnen nu worden 
onderscheiden in die welke als een convexe lineaire combinatie 
van andere punten van V zijn te schrijven, en die welke niet als 
een convexe lineaire combinatie van andere punten van V zijn 
te schrijven. De laatste noemt men de extreme punten van V. 


Voorbeelden: In tekening I van 7 zijn A, B, C, D en E extreme 
punten. In II is er een oneindig aantal extreme punten, namelijk 
alle punten op de omtrek. 


10. Uit 8 en 9 blijkt: Elk punt van een convexe verzameling V 
is te schrijven als een convexe lineaire combinatie van de extreme 
punten van deze verzameling. 


$ 4. Lineaire vergelijkingen 


1. Een lineaire vergelijking a's —=b wordt homogeen genoemd, 
indien 5 =O. Men spreekt van een niet-homogene vergelijking, 
indien b #0. 


2. Men zegt, dat een matrix van de orde m x n de vang r heeft, 
als de determinant van tenminste één van de vierkante deel- 
matrices van de orde 7 Xx r ongelijk aan Ois, terwijl de determinanten 
van alle vierkante matrices van hogere orden gelijk aan O zijn. 
Indien A van de orde m x n is, kan de rang van A hoogstens m 
zijn, wanneer mm < n; de rang van A kan hoogstens # zijn, wan- 
neer m > 1. 
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Een vierkante matrix van de orde m X m wordt singulier ge- 
noemd, indien de rang kleiner dan m is. Een singuliere matrix is 
dus een vierkante matrix waarvan de determinant gelijk aan O is. 
Een vierkante matrix van de orde m Xx m wordt niet-singubier 
genoemd, indien de rang gelijk aan m is. De determinant is dan 
ongelijk aan 0. 


3. Beschouw het stelsel van m homogene vergelijkingen met # 
onbekenden A« = 0 (A is van de orde m X 7). 

De vector « = 0 voldoet voor elke matrix A aan dit stelsel. 
Voorts geldt: Wanneer z een (eventueel van 0 verschillende) op- 
lossing is, is Az eveneens een oplossing voor elk reëel getal À. 


4. Stelling: Het stelsel van m homogene vergelijkingen met # 
onbekenden, Ax = 0, heeft een oplossing « # 0, wanneer en alleen 
wanneer de rang 7 van A kleiner dan # is. 


Bewijs: Per definitie is / <n. De rang van A kan immers niet 
groter zijn dan het aantal kolommen. Aan het gestelde is dus 
voldaan, wanneer het aantal onbekenden groter is dan het aantal 
vergelijkingen. Alleen het omgekeerde moet dus worden beschouwd. 
IL. Stel, 7 = n. Rangschik de rijen van A zodanig, dat de eerste 7 
rijen een niet singuliere deelmatrix van de orde » Xx » vormen. 
Dit kan volgens het gestelde. Deze deelmatrix wordt voorgesteld 
door A1; de deelmatrix, bestaande uit de overige m —n rijen 
van A wordt voorgesteld door As. 

Voor Ax = 0 kan men dan schrijven: 


_ le 


Aix = 0 
Asx = 0 


Omdat A; niet singulier is, bestaat Aj'. Uit Aix = 0 volgt dus 
Ar ‘Aix = Ar !0. 
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Dus %= 0. 
Uit r=n volgt dus #—= 0. Ma.w.: r <# is een noodzakelijke 
voorwaarde voor het bestaan van een oplossing « # 0. 


II. Stel r < n. Rangschik de rijen van A zodanig, dat de eerste 7 
elementen van de eerste 7 rijen een niet singuliere deelmatrix R 
van de orde 7 X 7 vormen. Daar de volgorde van de vergelijkingen 
in het stelsel Ax —= 0 niet ter zake doet, wordt de door de her- 
schikking van de rijen van A verkregen matrix gemakshalve even- 
eens voorgesteld door A. 

Zoals uit het gestelde volgt is elke deelmatrix van de orde 
(7 + 1) X (fr +1) die men uit R verkrijgt door hieraan r + 1 
elementen van een niet tot R behorende rij en 7 + 1 elementen 
van een niet tot R behorende kolom van A toe te voegen, singulier. 
Stel, men voegt (r + 1) elementen van de (7 + IJste rij van A 
en (r + 1) elementen van de (r + IJste kolom van A aan R toe. 
De aldus verkregen singuliere matrix van de orde (r + 1) X (rf + 1) 
wordt voorgesteld door C. De cofactoren van de laatste rij van G, 
die worden voorgesteld door Cosi Corio 00° ‚ Geir) 
zullen in het algemeen niet gelijk aan O zijn. In ieder geval is 
Cirie+i) #0, aangezien deze gelijk is aan |R| #0. 

Bepaal nu de vector z van de orde n X 1 zodanig, dat de eerste 
r + 1 elementen gelijk zijn aan bovengenoemde cofactoren van C, 
terwijl de laatste n — (r + 1) elementen gelijk zijn aan 0. 

Az is dan gelijk aan 


r+1 n 
Bauters + DX 440 ee en. 
jl j=r+2 


De tweede van deze sommen is vanzelfsprekend gelijk aan O voor 
alle z= 1 +: m. 

De eerste som is voor 2= 1 ++." r te beschouwen als de som 
van de elementen van een rij van C vermenigvuldigd met de co- 
factoren van de elementen van een andere, met name de (7 + iste, 
rij van CG. Deze uitdrukking is te beschouwen als de determinant 
van een matrix waarvan twee rijen aan elkaar gelijk zijn en is 
derhalve volgens $ 2.12 gelijk aan 0. 
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De eerste som is voor  —= (r + 1) gelijk aan |C] = 0. 

Voor £—= (ff 4e Aj er m is de eerste som te beschouwen als de 
determinant van een andere deelmatrix van de orde (7 + 1) X (r + 1) 
van A dan C. Deze determinant is per definitie eveneens gelijk aan 0. 
Derhalve is Az = 0. Voorts is z #0, omdat in ieder geval het 
(r + IJste element Cs ijgrri) — IRI #0. 

Wanneer 7 < n kan men dus een vector 2 #0 construeren, zo- 
danig dat Az — 0. M.a.w.: 7 < # is niet alleen een noodzakelijke 
doch tevens een voldoende voorwaarde voor het bestaan van een 
oplossing « #0 van Ax = 0. 


5. De vectoren ad ----- a in R® worden lineair afhankelijk 
genoemd, wanneer er » getallen xj «---- Xn, niet alle gelijk aan O, 
bestaan, zodanig dat 


xa) Jeen: + xam) —= 0. 


De getallen xj --*: Xn noemt men de gewichten in de lineaire 
onafhankelijkheidsbetrekking. 

Indien xj ……**: xn aan bovengenoemde eisen voldoen, voldoen 
Aje: Axn eveneens hieraan voor alle waarden van À. Naast 
deze oneindige verzameling van stellen gewichten A4 --*-- An, 
— co <Â << + vo, kunnen verschillende andere oneindige ver- 
zamelingen van stellen gewichten Ayj *:«"- Âym, — oo ZÂS + oo 
bestaan. Er bestaan dan tussen de vectoren ad ----- a ver- 
schillende lineaire afhankelijkheidsbetrekkingen. 

Indien vectoren ad ----- a niet lineair afhankelijk zijn, noemt 
men ze lineair onafhankelijk. 


6. Stelling: Wanneer en alleen wanneer de vierkante matrix A 
van de orde m x m singulier is, zijn de rijen en kolommen van Â 
lineair afhankelijk. 


Bewijs: 
IL. Als [Al = O0, bestaat er volgens 4 een oplossing # #0 van het 
stelsel Ax —= 0. Daar het vergelijkingenstelsel Ax —0 is te be- 
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schouwen als een lineaire combinatie van de kolommen van A 
met als gewichten xj ---*- Xn geldt: Als A singulier is, zijn de 
kolommen van A lineair afhankelijk. 

Hetgeen voor de kolommen van A geldt, geldt eveneens voor de 
rijen. Omdat de determinant van een matrix volgens $ 2.10 zowel 
volgens de elementen van een kolom als volgens de elementen van 
een rij kan worden ontwikkeld, is |A’| = |A. Derhalve is |A'| = 0, 
wanneer |A| =0. In dat geval bestaat er een oplossing y # 0 
van het stelsel A'y — 0. Als A singulier is, zijn dus eveneens de 
rijen van A lineair afhankelijk. 


II. Als |A| #0, bestaat er volgens 4 geen oplossing x #0 van 
Ax = 0. Ma.w.: de kolommen van A zijn in dit geval lineair 
onafhankelijk. 

Hetzelfde kan, op grond van het feit dat |A'| = |Al, worden be- 
wezen voor de rijen van À. 


7. Uit 6 volgt onmiddellijk het omgekeerde: Wanneer en alleen 
wanneer de kolommen van een vierkante matrix lineair afhankelijk 
zijn, is deze matrix singulier. 


8. Stelling: Wanneer de matrix A van de orde m x n de rang r 
heeft, zijn 7 en niet meer dan 7 rijen (kolommen) van A lineair 
onafhankelijk. 


Bewijs: 
TI. Wanneer A de rang r heeft, is een deelmatrix van de orde 
r Xr niet singulier. Er bestaan dus geen gewichten xj: Wp 


zodanig dat de gewogen som van de rijen of kolommen van deze 
deelmatrix gelijk aan O is. Voegt men aan elk van de rijen van 
de deelmatrix de overige (n — 7) elementen van de overeenkom- 
stige rijen van A toe, resp. voegt men aan elk van de kolommen 
van de deelmatrix de overige (mm — 7) elementen van de overeen- 
komstige kolommen van A toe, dan blijft het laatste uiteraard 
gelden. Er zijn dus # rijen en 7 kolommen van A lineair onafhan- 
kelijk. 
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II. Wanneer A de rang 7 heeft is elke deelmatrix van de orde 
(r +1) Xx (fr +1) gelijk aan O. Beschouw als eerste van deze 
deelmatrices de eerste (r + 1) elementen van de eerste (r + 1) 
kolommen van A. Stel de kolommen van deze deelmatrix voor 


door asin 0. Orgijgrsij Er bestaan dan gewichten 
Aj svenn Âir+1) zodanig, dat 
rtl 


DY Audis =O. 
i=1 


Beschouw vervolgens de eerste (r + 1) elementen van de tweede 
t/m de (7 + 2)-de kolom van A. Stel, de kolommen van deze 


deelmatrix voor door aja «+ +": Ar 1jr+ 2): Hiervoor geldt even- 
eens, dat er gewichten Âza «*«- Âogr+2) bestaan, zodanig dat 
+2 


Xi Alijn = 0. 
i=2 


Dit geldt voor elk stel opeenvolgende kolommen van (r + 1) ele- 
menten, tot en met: 
n 


NE  Am-erribern = 0. 
i=n—(r+1) 


Telt men deze 7 + 1 relaties op en stelt men de som van de bij de 
i-de kolom van (fr + 1) elementen behorende coëfficiënten A4 voor 
door p;, dan vindt men (r + 1) gewichten, niet alle gelijk aan O, 
zodanig dat 


r+1 
X wid; —= 0, 


i=1 

waarbij a1 «+*** Ari) de eerste (r + 1) rijen van A voorstellen. 
Hiermede is bewezen, dat de eerste (7 + 1) rijen van A lineair 
afhankelijk zijn. Hetzelfde kan op gelijke wijze worden bewezen 
voor elk ander stel van (r + 1) rijen. Voor een stel van (r + ) 
rijen kan men het gestelde op analoge wijze als hierboven aantonen 
door uit te gaan van deelmatrices van de orde (r + 4) X (r + {). 
Tenslotte levert men het bewijs voor de kolommen van A analoog 
aan het hierboven gegeven bewijs voor de rijen. 














46 Vectoren en lineaire vergelijkingen 


9. Uit 8 volgt rechtstreeks het omgekeerde: Wanneer 7 en niet 
meer dan 7 rijen en kolommen van A lineair onafhankelijk zijn, 
is de rang van A gelijk aan 7. 


10. Stelling: Elke verzameling van meer dan m vectoren in Rm 
is lineair afhankelijk. 


Bewijs: Men kan „ vectoren in R® beschouwen als de kolommen 
van een matrix A van de orde m x n. Daar de rang van A hoog- 
stens m is, is het stelsel A« —= 0 oplosbaar wanneer # > mm. 


11. De feitelijke oplossing van een stelsel homogene vergelijkingen 
kan worden herleid tot de oplossing van een stelsel niet-homogene 
vergelijkingen. Het stelsel Ax — 0 is oplosbaar, wanneer de rang 
van A van de orde m Xx 7 kleiner dan # is. Is de rang 7 dan bepaalt 
men eerst de 7 onafhankelijke kolommen van A (minstens één 
kolom valt hierbij af). Aan de variabelen, behorende bij de overige 
n — r kolommen, kent men willekeurige waarden, niet alle gelijk 
aan O, toe. Stelt men de eerste 7 kolommen voor door A1 en de 
gewogen som van de overige 7 kolommen door b, dan verkrijgt 
men een stelsel Aix — b. Door dit op te lossen verkrijgt men de 
waarden van de eerste 7 variabelen die tezamen met de willekeurig 
bepaalde waarden van de laatste „ — 7 variabelen een oplossing # 
van Ax — 0 opleveren. 

Op de oplossing van een stelsel niet-homogene vergelijkingen wordt 
in het volgende ingegaan. 


12. Stelling: Het stelsel van m niet-homogene vergelijkingen met 
n onbekenden Ax —=b, b #0 heeft een oplossing wanneer en 
alleen wanneer de rang van A gelijk is aan de rang van de matrix 
[Ab] bestaande uit de kolommen van A en uit b. 


Bewijs: Als de rang van A gelijk is aan 7, is een deelmatrix van 
de orde r X r niet singulier. Rangschik de rijen en kolommen van 
A zodanig, dat een dergelijke deelmatrix wordt gevormd door de 
eerste r elementen van de eerste 7kolommen. Stel de eerste 7 ko- 
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lommen van A voor door A1 van de orde m x 7. De rang van A1 
is dus gelijk aan 7. Volgens 4 betekent dit, dat er geen oplossing 
y #0 van het stelsel Ajy — 0 bestaat. 


I. Stel, de rang van [Ab] is 7. Dan is de rang van [A1b] van de 
orde m X (yr + 1) eveneens 7. Volgens 4 is er dan een oplossing 
z #0 van het stelsel [Ajb]z — 0. 

Het laatste element van deze oplossing z moet ongelijk aan O zijn. 
Zou dit niet het geval zijn dan zou er immers een oplossing, onge- 
lijk aan 0, van het stelsel Ajy = 0 bestaan, hetgeen in strijd 
met het gestelde is. 

Schrijft men nu voor [Aib]z = 0: 


n 
Xa + zb —=0, 


i=1 
dan volgt hieruit, omdat z # 0: 
1 r 
b= —— Xa 
R i=1 ° 
of wel 
n 
b= X xad, 
i=1 
waarbij 


1 5 
X1 FF — — 21 ane DE gp St 
of wel 
Ax — b. 


II. Stel de rang van [Ab] is (r + 1). De rang van [A1b] is dan 
eveneens (r + 1) en volgens 4 bestaat geen oplossing z #0 van 
het stelsel [Aib]z —= 0. 


Voegt men aan Aj 1,2, ----- of (# — 7) kolommen van A toe 
en stelt men de op deze wijze verkregen matrices voor door 
Ali =2"---: (n — 7 — 1)), dan bestaat wel een oplossing u: # 0 


van de stelsels [A;b]u; = 0. De orden van de matrices [A;b] zijn 








in 4 
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immers resp. m X (r + 2) +…+ m Xn, terwijl de rang van elk 
der matrices gelijk aan vr is. 

Daar volgens 8 de kolommen van A1 niet lineair afhankelijk zijn, 
maar de kolommen van elke deelmatrix bestaande uit (r + 1) 
kolommen van A wel, kan men elk der (n — 7) kolommen van A 
die niet tot A; behoren, als een lineaire combinatie van de kolom- 
men van A1 schrijven. Dit betekent, dat men b als een lineaire 
combinatie van de kolommen van Aj zou kunnen schrijven indien 
het laatste element van een der oplossingen u; gelijk aan O zou zijn. 
Het laatste is volgens het gestelde onmogelijk. Er bestaat dus geen 
oplossing u; van het stelsel [A;b]u; = 0 waarvan het laatste ele- 
ment ongelijk aan O is. Dus: het stelsel Ax = b is niet oplosbaar. 


13. Evenals in 5 e.v. na de bespreking van de voorwaarde voor de 
oplossing van een homogeen vergelijkingenstelsel enkele opmer- 
kingen werden gemaakt over de lineaire afhankelijkheid van vec- 
toren, kunnen hier na de behandeling van de voorwaarde voor de 
oplossing van een niet-homogeen vergelijkingenstelsel in het kort 
enkele eigenschappen van lineaire combinaties van vectoren worden 
genoemd. In $ 3,4 is een definitie van lineaire combinaties gegeven. 


14. In de eerste plaats geldt: Een vector b in Rm is te schrijven 
als een lineaire combinatie van ad ..... a in R® wanneer de 
rang van de matrix A, bestaande uit de vectoren ad --.-. an) 
als kolommen, dezelfde rang heeft als de matrix [Ab], bestaande 
uit de kolommen ad ----- a@® en b. Dit is niet meer dan een para- 
frase van de in 12 bewezen stelling. 


15. Het verband tussen lineaire afhankelijkheid en lineaire com- 
binaties is eveneens zonder meer duidelijk: 

Uit het feit, dat ad als een lineaire combinatie van a@ ----- a”) 
kan worden geschreven, volgt dat de vectoren ad... am) 
lineair afhankelijk zijn. Immers uit 


Nn 
ad = 5 ua 


i=3 
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volgt 


n 
De a® = 0, waarbij Âi == —l, Â2 = Haer 5 Àn = Un: 
i=1 
Het omgekeerde geldt niet: Uit het feit dat # vectoren lineair 
afhankelijk zijn volgt niet dat e/k daarvan als een lineaire combinatie 
van de (xn — 1) andere kan worden geschreven. Een bepaalde a@) 
kan alleen dan als een lineaire combinatie van de overige ad): « « -a®) 
worden geschreven, wanneer er minstens één lineaire afhankelijk- 
heidsbetrekking X% , wa® — 0 bestaat, waarin uj # 0. 


16. Een verzameling van m lineair onafhankelijke vectoren in Rm 
noemt men een basis in RP, 


Voorbeeld: De eenheidsvectoren [1 O--::- O) vaars [O---<: O0 1}, 
elk bestaande uit m elementen, vormen een basis in R. 


17. Stelling: Elke vector in Rm kan op één manier als een lineaire 
combinatie van een stel basisvectoren worden geschreven. 


Bewijs: Een stel basisvectoren ad -.--- am in Rm vormt een 
niet-singuliere matrix A; van de orde m x m. Voegt men hieraan 
een andere vector b van de orde m Xx 1 toe, dan verkrijgt men een 
matrix van de orde m X (m + 1). De rang hiervan ism. Volgens 12 
kan de laatste vector dus als een lineaire combinatie van de eerste 
m worden geschreven. De lineaire combinatie is eenduidig. De 
gewichten zijn namelijk gelijk aan « —= Aj 'b (zie $ 2, 15). 


18. Voor de feitelijke oplossing van het stelsel Ax — b elimineert 
men eerst ‘overbodige’ vergelijkingen en variabelen. Is m = n= 7 
dan kan men onmiddellijk één van de hierna te bespreken technie- 
ken toepassen. Is m > 7 dan beschouwt men alleen 7 vergelijkingen 
waarin lineair onafhankelijke rijen van A voorkomen. 

Stel bijv. een stel van 7 van dergelijke rijen en de bijbehorende 
elementen van b voor door [A1bi]. Deze matrix van de orde 
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r X (m1) heeft de rang 7, terwijl m > 7. Het stelsel is dus 
oplosbaar. 

Stel voorts een andere rij van A (niet tot A; behorend) en het 
bijbehorende element van b voor door [4,b;}. Daar de matrix 
A1b1 
a,b; 
[a4b;] als een lineaire combinatie van de eerste 7 rijen [A1bi] wor- 
den geschreven. Dus: 


(ab) = p[A1bi]. 
Stel, « is een oplossing van Aix = bi. Er geldt dan 


[Aibi][e — 1] = 0. 


van de orde (7 + 1) X (m +1) de rang 7 heeft, kan 








Volgens het bovenstaande is dan 

[asb;][e — 1] = u[Aibijl[x — 1] = 0. 
Dus: 

ajx — bj, 


m.a.w. % is eveneens een oplossing van de j-de vergelijking. 
Hieruit blijkt dat men voor het geval dat m > 7, de ‘afhankelijke’ 
(m — r) vergelijkingen buiten beschouwing kan laten. Het over- 
blijvende stelsel van 7 vergelijkingen met n onbekenden lost men op. 
Een oplossing hiervan voldoet eveneens aan de buiten beschouwing 
gelaten vergelijkingen. 

Is tenslotte „ > 7, dan beschouwt men alleen 7 lineair onafhanke- 
lijke kolommen van A en de daarbij behorende variabelen. Aan de 
overige variabelen kent men willekeurige waarden toe en vermenig- 
vuldigt deze met de bijbehorende kolommen. De aldus verkregen 
vectoren trekt men af van b. Het probleem is dan herleid tot het 
oplossen van een niet-homogeen vergelijkingenstelsel met 7 on- 
bekenden (waarbij de oplossing wordt bepaald door de gekozen 
waarden van » — r variabelen). 

In elk van de drie gevallen kan de oplossing van een niet-homogeen 
vergelijkingenstelsel dus gereduceerd worden tot de oplossing van 
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een stelsel waarin het aantal vergelijkingen gelijk is aan het aantal 
onbekenden. 


19. De meest gevolgde methode voor het numeriek oplossen van 
m niet-homogene vergelijkingen met m onbekenden is het succes- 
sievelijk elimineren van variabelen. Hierbij gaat men als volgt 
te werk: 

De twee eerste vergelijkingen vermenigvuldigt men met zodanige 
getallen, dat de twee coëfficiënten van een willekeurige variabele 
gelijk worden. Vervolgens trekt men de vergelijkingen van elkaar 
af. Hetzelfde doet men met de tweede en de derde vergelijking, 
de derde en de vierde, enz. Op deze wijze verkrijgt men tenslotte 
een stelsel van (m — 1) vergelijkingen met (m — 1) onbekenden. 
Op het aldus verkregen, beperktere stelsel past men deze techniek 
opnieuw toe. Men verkrijgt dan een stelsel van (m — 2) verge- 
lijkingen met (m — 2) onbekenden. Zo doorgaande vindt men uit- 
eindelijk één vergelijking met één onbekende. Uit deze vergelijking 
kan de waarde van de laatst overgebleven variabele worden bere- 
kend. Door deze waarde in één van de vergelijkingen van het 
‘voorgaande!’ stelsel van twee vergelijkingen met twee onbekenden 
te substitueren vindt men de waarde van de op één na laatst 
overgebleven variabele. Aldus terugkerende kan men de waarden 
van alle variabelen bepalen. 


Voorbeeld: 


2x + 4x2 =BI3| 6x1 + 1249 = 24 
9x + 342 == 7 |4| 20x1 + 1243 —= 28 


14x1 = 4 


X1 = 8, Xg = En 


20. De oplossing van m niet-homogene vergelijkingen met m 
onbekenden kan tenslotte ook worden bepaald met behulp van 
de z.g. Regel van Cramer. 
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Uit Ax = b, waarin A vierkant en niet-singulier is, volgt: 
i Ag veert Ami bi 


nag ss En bm 
De j-de variabele is dus gelijk aan 


_ biAu + beAa Herr: + bmAmt 
E |Al 





Xi 


De teller van deze breuk kan men beschrijven als de determinant 
van een matrix, waarvan alle kolommen behalve de d-de gelijk zijn 
aan de overeenkomstige kolommen van A, terwijl de #-de kolom 
gelijk is aan b. Stelt men de aldus gedefinieerde matrix voor door 
Ai, dan wordt de oplossing van het vergelijkingenstelsel gegeven 
door de volgende formule, die men aanduidt als de Regel van 
Cramer 


|Ail 


Ki Pe. 
|Al 


Men kan x; dus vinden door beide determinanten uit te rekenen. 








Il. Lineaire programmering als 
wiskundig probleem 


Na een uiteenzetting van de algemene kenmerken van het lineaire 
programmeringsvraagstuk wordt in dit hoofdstuk een algemene 
oplossingstechniek, de z.g. simplexmethode, behandeld. In $ 1 wordt 
nagegaan, op welke wijze een gegeven vraagstuk, waarin zowel 
gelijkheden als ongelijkheden voorkomen, kan worden herleid tot 
een vraagstuk met uitsluitend gelijkheden. Deze herleiding is 
noodzakelijk, omdat de te bespreken simplexmethode alleen op 
stelsels van gelijkheden toepasbaar is. Er wordt in $ 1 bewezen, 
dat de oplossing van het gegeven, oorspronkelijke vraagstuk on- 
middellijk uit de oplossing van het herleide volgt. In $ 2 worden 
enkele stellingen omtrent de oplossingen van het lineaire program- 
meringsvraagstuk besproken. De volgende paragraaf is gewijd aan 
een bepaalde techniek — de substitutie van vectoren in een basis — 
waarvan in $ 4 gebruik moet worden gemaakt bij de afleiding van 
de simplexmethode. Deze simplexmethode is de meest toegepaste 
en doorgaans ook de eenvoudigste methode ter oplossing van 
lineaire programmeringsproblemen. In bijzondere gevallen gaan 
andere methoden soms met minder rekenwerk gepaard; de in 
IV $3 te behandelen techniek is hiervan een voorbeeld. Op deze 
andere methoden wordt hier echter verder niet ingegaan. In 85 
wordt nog een laatste kenmerk van de simplexmethode besproken, 
teneinde de theorie volledig sluitend te maken. Bij eerste lezing 
doet men er wellicht goed aan, dit betrekkelijk lastige onderwerp 
over te slaan. 


$ 1. De oorspronkelijke en de herleide vorm 


1. De in 1 $2 geformuleerde lineaire programmeringsproblemen 
kunnen op twee manieren in matrix-notatie worden geschreven. 
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Men kan de vraagstukken in de eerste plaats als volgt omschrijven: 
Bepaal de vector « zodanig, dat 


Ax < b resp. Ax > b 
«>= 0 «>0 
C'£ —= max C'« —= min 


Men kan ook stellen: 


Bepaal de gewichten x1 «---- Xn zodanig, dat 
n Nn 
N40) <b resp. ad > b 
i=1l dal 
Kie Xn > 0 ra Xn br 0 
n n 
Si CIX —= Max Di Cx: = Min 


i=1 i=l 


2. De betrekkingen Ax <b en « >0 (resp. Ax >b en « > 0) 
noemt men de voorwaarden waaraan de oplossingen moeten vol- 
doen. De relatie c'x wordt aangeduid als de maximerings-, resp. de 


minimeringsfunctie. 
Aangenomen wordt, dat A van de orde m Xx # is. De kolommen 
ad... a en b zijn dan punten in Rm. De mgee oplos- 


singen % zijn punten in R». 


3. De betekenis van de tekens < en > is in II $ 1, 15 verklaard. 
Bij een bepaalde oplossing «* zullen sommige vergelijkingen vol- 
doen aan a;x* < bj (resp. > bj), andere aan a;x* — bj. Als voor- 
waarden voor nog te bepalen oplossingen #, zullen sommige ver- 
gelijkingen luiden a; < b; (resp. > bi), andere a;« — bj. In het 
eerste geval heeft men de keuze tussen gelijkheid en ongelijkheid; 
in het tweede geval bestaat deze keuze niet. 


4. De in $4 te bespreken oplossingstechniek kan alleen worden 
toegepast op stelsels van gelijkheden. In het volgende wordt aan- 
getoond, dat elk vergelijkingenstelsel, ongeacht of in sommige 
vergelijkingen het <, het > of het = teken geldt, als een stelsel 
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van gelijkheden kan worden geschreven en dat de oplossing van 
het oorspronkelijke stelsel direct uit die van het tweede volgt. 


5. Elke ongelijkheid kan als een gelijkheid worden geschreven door: 

— een nieuwe variabele in te voeren 

— deze bij de term voor het ongelijkteken op te tellen (resp. ervan 
af te trekken) 

— de eis te stellen, dat de ingevoerde variabele positief is. 

In plaats van a;x < b; kan men immers schrijven 


ax + Â= bi 
Â =0. 


In plaats van a;x > bj kan men schrijven 


aix — Â= bi 
A>=0. 


De nieuw ingevoerde variabele A kan men een verschilvariabele 
noemen (Engels ‘slack variable’). 


6. Van de in 5 weergegeven methode wordt bij het oplossen 
van lineaire programmeringsvraagstukken gebruik gemaakt: Be- 
schouw de voorwaarden van een maximum-vraagstuk Ax < b 
en « >=0. Stel, s van de voorwaarden zijn ongelijkheden; de 
overige (m — s) zijn gelijkheden. Aan de s ongelijkheden voegt 
men dan s verschilvariabelen Àj «---: Âs =O toe. Het stelsel 
wordt dan: 


Ax + Sed = b 


i=l 


Hierbij zijn ed ----- e® eenheidsvectoren in Rm, 
Geeft men vervolgens [& Â1--** As]’ weer door y, [A ed ----: e(»)] 
door B en [c’ 0:---- 0} door d’, dan kan bovenstaand stelsel 
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worden voorgesteld door 


By = b 
y=0 
d'y —= max. 


Het laatste is een lineair maximeringsvraagstuk waarvan alle 
voorwaarden gelijkheden zijn. Heeft men het tweede vraagstuk 
opgelost, dan kan de oplossing van het eerste, oorspronkelijke 
probleem hieruit zonder meer worden afgeleid. Immers: 

a) De eerste m elementen van een oplossing y* van By = b en 
y >0 voldoen aan Ax <b en x > 0. 

b) Daar de laatste s elementen van d’ (de coëfficiënten van de ver- 
schilvariabelen) alle gelijk aan O zijn, voldoen de eerste elementen 
van een vector y* waarvoor d'y — max, tevens aan C'4 — max. 


7. Wat voor een maximeringsvraagstuk is opgemerkt, geldt mu- 
tatis mutandis eveneens voor een minimeringsvraagstuk. In plaats 
van Ax > b, « > 0 en C'x = min, stelt men 


8 
Ax — ed =b 


i=l 


€'& + OÂr + -ee + OA; = min. 


Heeft men het tweede probleem opgelost, dan is de oplossing van 
het eerste gelijk aan de eerste n elementen van de oplossing van 
het tweede. De waarden van de verschilvariabelen in het tweede 
probleem zijn voor het eerste irrelevant. 


8. Het is thans zonder meer duidelijk dat ook een stelsel, bestaande 
uit vergelijkingen a;x < bi, uit vergelijkingen a; x > bj en uit ver- 
gelijkingen aix — bj, als een stelsel van louter gelijkheden kan 
worden geschreven. In de eerste soort vergelijkingen introduceert 
men verschilvariabelen met een coëfficiënt + 1, in de tweede soort 
vergelijkingen introduceert men verschilvariabelen met een coëffi- 
ciënt — 1 en de derde soort vergelijkingen laat men onveranderd. 
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De oplossing van het oorspronkelijke vraagstuk is gelijk aan de 
waarden van de oorspronkelijke variabelen in het tweede vraagstuk. 


9. Het gegeven lineaire programmeringsvraagstuk kan men aan- 
duiden als de oorspronkelijke vorm. Het hieruit afgeleide probleem, 
dat men verkrijgt door de ongelijkheden in gelijkheden om te zetten 
via de introductie van verschilvariabelen, wordt aangeduid als de 
herleide vorm. 


10. In het volgende zal voornamelijk aandacht worden geschonken 
aan de herleide vorm. Om de invoering van nieuwe letters te voor- 
komen, wordt de herleide vorm door dezelfde symbolen voorgesteld 
als de oorspronkelijke. Dus: 


Ax = b 
«> 0 
c'x = max, resp. min. 


Steeds wordt aangenomen, dat A van de orde m Xx # is. 


$ 2. Bruikbare en optimale oplossingen 


1. Een vector « die voldoet aan de voorwaarden Ax =b en 
x > 0 van de herleide vorm, noemt men eeh bruikbare oplossing 
(in het Engels ‘feasible solution’). 


2. In het algemeen voldoet aan een stelsel lineaire vergelijkingen 
Ax = b een oneindig aantal oplossingen «, als A van de orde 
m X nisen n > m. Hieromtrent kan verder niets worden gezegd. 
Er is in het algemeen ook een oneindig aantal bruikbare oplossingen, 
voldoende aan Ax=b en « > 0. Hieromtrent kan men echter 
een aantal stellingen afleiden. 


3. Stelling: De bruikbare oplossingen vormen een convexe ver- 
zameling V in Rm, 
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Bewijs: Voor Ax = b kan men schrijven X% , #4a® — b. 
Stel, u en wv zijn twee bruikbare oplossingen. Dus: 


n 
wad =b en wy =O 


iel 


n 

va) =b en vu >= 0. 

i=1 
Beschouw een convexe lineaire combinatie van deze vectoren 
z=(l — uu uv 0<u<l. Er geldt: 


n 
5 440 — 


i=l 


py (1 Ts u): + uv;ja® == 
i=1 


i=l 


n n 
(1 — u) E wad + u SY ua® — b, 
i=1 


terwijl tevens: 
z=(l — uu +uv >= 0. 


Als u en wv bruikbare oplossingen zijn is een convexe lineaire com- 
binatie hiervan dus ook een bruikbare oplossing. M.a.w.: V is een 
convexe verzameling. 


4. Stelling: Een bruikbare oplossing eV is een extreem punt 
van V, wanneer en alleen wanneer de positieve elementen van « 
(de overige zijn gelijk aan O) de gewichten zijn van lineair onaf- 
hankelijke vectoren a® in %%, #4a® = b. 


Bewijs: 

IL. Stel, « is een extreem punt van V. Stel voorts, de eerste k ele- 
menten van % zijn positief, de laatste n — k gelijk aan 0. Er geldt 
dan: 


n k 
A) = Y HA) — b. 
i=1 i=1 
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Stel nu, dat de vectoren ad ----- a®) lineair afhankelijk zijn. 
Er bestaan dan k getallen Aj +++" Ax, niet alle gelijk aan O, zo- 
danig dat 

k 

Ss ha® = 0. 

i=1 


Voor elke constante c geldt in dat geval: 


k k k 
Eua® Head =D (ús + cdj)a® =b 
i=1 i=1 i=1 


en 

k k k 

Sad — Cc 5 had = 5 (x; — chj)a® = b 

i=l i=1 i=1l 
Daar x4 > 0 kan men, ongeacht de waarden van Âj «-:*: Àr, de 
constante c altijd zodanig kiezen, dat 

x: + ch; > 0 
en voor #= 1 ++"" k. 

Ki — ch > 0 
Voeg nu aan elk van de vectoren [(%1 + cli) «+ (xx + cÂp)] en 
[(x1 — chi) «rr (ex — cÂi)] (n — K) nullen toe en stel de aldus 


verkregen vectoren met # elementen voor door u en v. Daar ook 
EL wad =ben XL, ua) =b, is ue V en ve V. Nu is 


Ju tw = 
Al(x1 + chi) eene (xx J cÂx)O … eee 0] + 
Hw — chi) «°° (er — cÂ)O «---- OJ = 4. 


Hiermede is « geschreven als een convexe lineaire combinatie van 
ueV en ve V. Het laatste is echter in strijd met de veronder- 
stelling, dat « een extreem punt van V is. Dus: Wanneer # een 


extreem punt van V is, kunnen a® «---- a% niet lineair afhan- 
kelijk zijn. 
IL. Stel, al neer a zijn lineair onafhankelijke punten in Rm. 


Neem onder deze veronderstelling aan, dat 
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waarvoor geldt X£, xa® — b, niet een extreem punt van V is. 
Het laatste betekent dat « als een convexe lineaire combinatie 
van twee andere punten ue V en ve V kan worden geschreven: 
x= (ll — uu tuv, 0OSpu<l. Op grond van de vorm van & 
moeten ook de laatste (n — k) elementen van u en w gelijk aan O 
zijn. Er geldt derhalve 


k k k 
NAO = Suad —= Sud =— b. 


t=l i=l i=1l 
Dus 

k k 

5 (% — W) a} = 5 (44 — vj) aÌ = 0. 

i=1 tel 
Daar ad ----- g@ lineair onafhankelijk zijn, zijn de gewichten 
in beide bovenstaande betrekkingen gelijk aan 0. 
Dus: 

KL Uijen Ak = Uk 
en 

KU Kk = Uk 


Daar de laatste (n — k) elementen van #, u en wv gelijk (aan 0) 
zijn, is: 

= U ==. 
Dus: # kan alleen geschreven worden als een convexe lineaire 


combinatie van zichzelf. M.a.w.: Als ad ----- a lineair onaf- 
hankelijk zijn, is # een extreem punt van V. 


5. Aangezien volgens II 84,10 niet meer dan m vectoren a® 
in R® lineair onafhankelijk kunnen zijn, volgt onmiddellijk uit 4: 
Een extreem punt van V kan hoogstens m positieve elementen 
_ hebben (de overige n — k zijn gelijk aan 0). 


6. Stelling: Het aantal extreme punten van V is eindig. 


Bewijs: Het aantal kolommen van A is eindig. Uit de kolommen 
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van A kan derhalve ook slechts op een eindig aantal manieren een 
stel lineair onafhankelijke kolommen worden gekozen. 


7. Beschouw nu de volledige herleide vorm 


Ax = 0 
x > 0 
C'& —= max, resp. min. 


Hiervoor is bewezen dat de oplossingen die aan de voorwaarden 
voldoen (de bruikbare oplossingen) een convexe verzameling V 
met een eindig aantal extreme punten vormen, en dat de positieve 
elementen van de extreme punten de gewichten van lineair onaf- 
hankelijke vectoren a® zijn. 

De lineaire funktie C'4 —= max, resp. min kan worden beschouwd 
als een funktie van de punten van V. Een bruikbare oplossing %, 
waarvoor C'x optimaal (max of min) is, noemt men een optimale 
oplossing van het lineaire programmeringsvraagstuk. 


8. Stelling: Het optimum van c'« wordt bereikt in een extreem 
punt van V. 


Bewijs: Volgens 6 is het aantal extreme punten van V eindig. 
Stel de extreme punten voor door e® ----- e® in R®. 

Stel, c’# is optimaal in het punt # —= u. Aangezien elk punt van 
een convexe verzameling volgens II, $ 3, 10 kan worden geschreven 
als een convexe lineaire combinatie van de extreme punten, is 


8 8 
u= YX mwed, waarbij X uy == l en uy > 0. 
j=l j=l 
Elk element w; van u is dus gelijk aan 
8 


== Eu, 
j=l 


waarbij == 1 +." n en ef) = het j-de element van e@). 
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Voor de maximum-, resp. minimum-waarde van de funktie kan 
men dus schrijven 


8 
= Y wyC'ed) = max, resp. min. 
j=1 
Beschouw nu alleen het maximum-vraagstuk. (Het verdere bewijs 
verloopt voor het minimum-vraagstuk mutatis mutandis op de- 


zelfde wijze). 
Stel, e® is van de extreme punten datgene waarvoor de funktie- 


waarde maximaal is. Dan is 
8 
2 ujc'ed < E uee) —=c et) De uy = ce). 
j=1 j=1 
Dus: 
c'u < c'ek). 


Daar volgens het gestelde c'u — max kan in de laatste betrekking 
alleen het gelijkteken gelden. Dus: u — e&). 


9. Stelling: Wanneer de lineaire funktie c'# haar optimum. in 
verschillende extreme punten ed .---- e® van V bereikt, dan 
wordt het optimum ook bereikt in elke convexe lineaire combinatie 
van ed ----. en). 


Bewijs: De funktiewaarde van een convexe lineaire combinatie 
ud van ed -…-.- em (Sl, en m0) is 
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gelijk aan: 


n 
c'u == Di Ci: 


i=l 
n h ‚ h 
i=l j=l j=1 


h 
— E uy (max €'#) 
=1 


10. Samenvattend kan nu worden opgemerkt: 

— Het totaal van de bruikbare oplossingen « is een convexe ver- 
zameling V met een eindig aantal extreme punten in R®. 

— De extreme punten zijn de gewichten in lineaire combinaties 
van onafhankelijke kolommen a® in Rm, 

— Het optimum van de maximerings-, resp. minimeringsfunktie 
wordt bereikt in een extreem punt of in een convexe verzameling | 
van een aantal extreme punten van V. | 

Op deze eigenschappen van het lineaire programmeringsvraagstuk | 

is de in $ 4 te behandelen simplexmethode gebaseerd. 


11. Als illustratie in R? van de in deze paragraaf behandelde 
begrippen wordt hier nog eens het in I, $ 1,8 genoemde maxi- 
meringsprobleem beschouwd: 


100x1 + 4042 < 4000 
20x1 + 60x2 < 1800 

X1, %2 > O0 
W =— 20x1 + 25x3 —= max. 








De oplossingen & van dit vraagstuk, bestaande uit twee elementen, 
kunnen worden voorgesteld door de punten van het platte 
vlak. Daar « > 0 moet zijn, liggen alle punten in het positieve 
kwadrant. 

Van de vectoren x in R?2, waarvan de elementen voldoen aan 
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100x1 + 40x3 =—= 4000 en xj, 42 > 0, is de ene uiterste 


4000 = f40|; 
100 | 0 | 
0) 
de andere uiterste is 
0) = ff 0 
4000 fnal 
40 


De verzameling van alle vectoren in R?2 waarvan de elementen 
aan bovenstaande gelijkheden voldoen is 


40f en | O E 
het segment tussen | el f il (zie II, 8 3, 6). 


De verzameling van de vectoren #, waarvan de elementen voldoen 
aan 100xj + 40xa < 40 000 en #1, x2 >= 0 kan worden voorgesteld 
door het oppervlak van de driehoek, gelegen tussen de beide assen 


en het segment tussen 6) ie jk 


{ al . Dit oppervlak is hieronder 


gearceerd voorgesteld. 


tweede 
element 


100 


50 


50 100 eerste element 
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De verzameling der vectoren «, voldoende aan 
20x1 + 60x2 < 1800 en x1,x2 > O 


kan worden voorgesteld door het oppervlak van de driehoek, gelegen 
tussen de beide assen en het segment tussen de punten. 


1800 Ì —= 5 en e) = Lal 
20 0 1800 30 
0) 60 


Dit oppervlak is hieronder, tezamen met het eerste, gearceerd 
voorgesteld. 


tweede 
element 






50 





2555 
ee 
eee, 


4, ", 





ee; 
25 
25 
25 
25 
25 
25 
25 
523 





OON 
AN ANKN 


0 A 50 100 eerste element 


Pen 





Het is duidelijk, dat de vierhoek OABC de convexe verzameling 
van bruikbare oplossingen aangeeft. De extreme punten hiervan 
zijn O, A, B en C. 

Welke van de extreme punten, of welke convexe lineaire combina- 
ties van extreme punten, een optimale oplossing is resp. optimale 
oplossingen zijn, kan in dit geval ook grafisch worden bepaald: 
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Een verzameling van bruikbare oplossingen welke alle een zelfde 
winst Wo opleveren, wordt voorgesteld door het segment tussen 


Wo en 0 
20 Wo 
0 25 


Zoals gemakkelijk is in te zien worden verzamelingen met verschil- 
lende winst Wo en Wi voorgesteld door evenwijdig lopende lijn- 
stukken. Naarmate de winst groter is ligt het betrokken lijnstuk 
verder van de oorsprong verwijderd. Teken nu de verzameling 
punten met een bepaalde winst, bijv. W == 1000 (zie onder- 
staande tekening). Het blijkt dat binnen de convexe verzameling, 
aangegeven door OABC nog punten rechts van deze lijn liggen. 
Er zijn dus oplossingen met een grotere winst dan f 1000. De 
meest naar rechts gelegen lijn, evenwijdig aan de eerste, gaat door 
het punt B. Dit punt geeft dus de optimale oplossing aan. Hierbij 
is 41 =ca32ton en x2 — ca 20 ton. 


tweede 
element 





A 50 409 eerste element 
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Het merkwaardige van dit resultaat is, dat men voor het bereiken 
van de maximale winst beide artikelen moet opslaan, hoewel in 
de gegevens geen beperkingen t.a.v. de afzetmogelijkheden zijn 
opgenomen. Zonder toepassing van de hiervoren besproken tech- 
niek zou men niet gemakkelijk tot deze conclusie zijn gekomen. 
In dit geval is er slechts één optimale oplossing. Men zou een on- 
eindig aantal convexe lineaire combinaties van twee extreme pun- 
ten hebben verkregen, wanneer de lijnen welke punten met gelijke 
winst aangeven, evenwijdig zouden hebben gelopen aan AB of 
aan CD. 


$ 3. Substitutie van vectoren in een basis 


1. Stela@® ----- a@) zijn n vectoren in de m-dimensionale ruimte. 
Neem aan, dat de eerste m lineair onafhankelijk zijn. Deze vormen 
dus volgens II, $ 4, 16 een basis in Rm. 

Men kan ad ----- a® volgens II, $ 4, 17 op één en niet meer dan 
één manier schrijven als lineaire combinaties van ad ..-.- a. 
Als men deze lineaire combinaties voorstelt door 


m 
a% == DN uya® ] == 1 NONE n, 
i=1 


dan kunnen de coëfficiënten wij worden weergegeven zoals in 
Tabel 1. ’ 


Tabel 1 
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2. Men kan een vector a®) uit de basis ad ----- a vervangen 
door een vector al) van de overgebleven a@mtD .…... a. Aan- 
genomen, dat ad. as... am) evenals ad - ---- a... am) 
lineair onafhankelijk zijn, kan men ad ----- a ook als lineaire 


combinaties van deze nieuwe basisvectoren uitdrukken. De coëffi- 
ciënten van de nieuwe combinaties worden in Tabel 2 weergegeven. 


Tabel 2 





3. Stelling: Vervangt men a@) in de oorspronkelijke basis door 
a® dan vormen ad ----. a... am) wederom een basis 
wanneer en alleen wanneer uns #0 in Tabel 1. se 


Bewijs: ad ……... a... am vormen een basis in R® als 
zij niet lineair afhankelijk zijn, d.w.z. wanneer er niet ge- 
wichten Aj ---*- Agere Àm, niet alle gelijk aan O, bestaan waar- 
voor geldt: : 


E HAD + A,4© — 0. 


ikh 


Volgens tabel 1 is 


m 
a) —= uisa®. 
i=l 
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Substitutie van de tweede uitdrukking in de eerste geeft 


m 
MAO + Ae S m000 — 
ikh i=l 
DN (A + Âgpis)a® + Asunsa) == 0. 
ih 
De laatste uitdrukking is een lineaire combinatie van de vectoren 
ad... ah) am). Daar deze vectoren lineair onafhanke- 
lijk zijn, moeten alle coëfficiënten in deze combinatie gelijk aan O 
zijn. Met name is dus ook 


Astins == 0. 
Is ups # 0, dan moet dus A; = O zijn. Wanneer Â; — O, moeten de 
coëfficiënten Ag «---- Am, i #h in de eerste vergelijking eveneens 
gelijk aan O zijn. Dit betekent, dat a «--- a atm) lineair 


onafhankelijk zijn. 

Is daarentegen uns — O, dan kan men voor As een waarde ongelijk 
aan O kiezen. Volgens de eerste vergelijking betekent dit, dat 
ad... a am) lineair afhankelijk zijn. 


4. Stelling: 

(1) De A-de rij van Tabel 2 wordt gevonden door de h-de rij van 
Tabel 1 te delen door gns- 

(2) De overige rijen van Tabel 2 worden gevonden door bij elk 
van de overeenkomstige rijen van Tabel 1 zoveel maal de h-de rij 
op te tellen, dat op de s-de plaats een O komt te staan. 


Bewijs: In symbolen uitgedrukt luidt het gestelde 





Hhj 
(Ĳ) Phi. = : 1 = | eee. n 
Uhs 
Hit b 
(2) py ag Fee n 
hs 
voor elke 4 = 1 +: m,i #h. 


Volgens (2) is het s-de element van een willekeurige rij 2 # h 
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van Tabel 2 inderdaad gelijk aan O, omdat voor j = s geldt: 





d 
Pis S Mis ed = 0. 
Uhs 


Er moet nu bewezen worden, dat de volgens (1) en (2) bepaalde 
elementen van Tabel 2 de coëfficiënten van de verschillende lineaire 
combinaties zijn, d.w.z. dat 

ad —= 5 pjad + psa) voor elke j= | +": n. 


ih 


Dit bewijst men door deze vergelijking als volgt te herleiden 





as 


ad =Y (vo iE nja + EA 
Hhs 


ikh hs 





d 
=E ( Hij — Le nao je ed Xi isa® + pna) 
Uhs ih 


ikh Uhs 


=E (ut en een Uni bn and Me) + nja) 


ih 


=S uG4Ò + urjad 
ikh 


i=l 


De laatste uitdrukking is de s-de lineaire combinatie volgens 
Tabel 1. Hiermede is het gestelde bewezen. 


5. Beschouw het vergelijkingenstelsel Ax = b, waarvan de coëffi- 
ciëntenmatrix A van de orde m Xx n is. Voeg aan A de eenheids- 
matrix van de orde m Xx m en tevens b toe. De aldus verkregen 
matrix kan worden beschouwd als een tabel zoals Tabel 1, wanneer 
men zich voorstelt, dat alle kolommen ervan als lineaire combina- 
ties van de eenheidsvectoren zijn uitgedrukt. Elke vector is immers 
een lineaire combinatie van de eenheidsvectoren met de elementen 
van de betrokken vector als coëfficiënten. 

_ Volgens de in deze paragraaf besproken techniek kan men nu de 
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eenheidsvectoren in de basis successievelijk door kolommen van A 
vervangen. Alle kolommen, incl. b, worden hierdoor als lineaire 
combinaties van m kolommen van A geschreven. Dit betekent, 
dat het stelsel Ax — b is opgelost; de laatste kolom van de tabel 
geeft de oplossing aan. 


Voorbeelden: In onderstaande voorbeelden is in de achtereenvol- 
gende tabellen het element uns, dat in de gesubstitueerde rij en 
kolom staat, omlijnd. Men noemt dit element wel het spil-element. 
Stel, het volgende vergelijkingenstelsel is gegeven 

4x A 6xa + 1243 — 22 

6x1 4 3xa + 15x3 —= 45 

2x H+ 7x — A3 == — 15 
Om dit vraagstuk op te lossen zijn in de eerste van onderstaande 
tabellen de coëfficiënten van dit stelsel en de eenheidsvectoren in 
R3 weergegeven. In de tweede tabel is a in de basis vervangen 
door a@®; in de derde tabel is a® vervangen door a® en in de 
vierde tabel tenslotte a® door a@. De oplossing van het verge- 
lijkingenstelsel blijkt dan te zijn xj = 4, 42 == —3 en %3 = 2. 


ad a) a) a a5) a(6 b 


a 
a(5) 
a(6) 


a 
a?) 
a{6) 


ad 
a) 
a(6) 
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an a) a) a a65) a) b 


at) 
a) 
a) 





Uit dit voorbeeld blijkt, dat het niet noodzakelijk is, de eenheids- 
vectoren expliciet te noemen en deze als lineaire combinaties van 
de basisvectoren uit te drukken. Dit is in het voorgaande uit- 
sluitend gedaan om de berekeningsmethode duidelijk voor te stel- 
len. In het volgende voorbeeld zijn de eenheidsvectoren weggelaten. 
Gegeven is het volgende stelsel 


24x1 + 122 — 7x3 3 


18x1 + 942 — 343 == 9 
2x1 + Mx + Bxzg —= 42 


In de eerste tabel zijn alleen de coëfficiënten van deze vergelijkin- 
gen weergegeven. Vervolgens worden resp. de vectoren a@, al) 
en a@ in de basis gebracht. De oplossing is dan xj = 0, x2 = 2 
en X3 = 3. 
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an 1 ) 0 0) 
a(2) (€) Ì 0 2 
a8 0 0) 1 3 


$4. De simplexmethode 


1. Beschouw de herleide vorm van een lineair programmerings- 
vraagstuk 


Ax = b 
«> 0 
C'x = optimaal (max of min). 


2. Aangenomen wordt, dat een tabel is gegeven, waarin elke 
kolom van A als een lineaire combinatie en b als een niet-negatieve 
lineaire combinatie van m lineair onafhankelijke kolommen van 
A is geschreven. 

Deze m kolommen van A noemt men de basisvectoren in de be- 
trokken tabel. De bijbehorende variabelen noemt men de basis- 
variabelen. 

Ter illustratie wordt Tabel 3 beschouwd. Voor het gemak is aan- 
genomen, dat de eerste m kolommen van A in dit geval basis- 
vectoren zijn. 





Tabel 3 
‚ a----a).---qln) amtD ..... as) ----- a b 
ag) l ar vet 0 vee ee 0 Bron 77 His aes. Hin % î 
8 8 : , ; : : de 
am 0 messe Ì EEE Or 0 Hnn+ 1) VET Hns vee Uhn %, 


an) Oo ermee Omen Ì Bman Hams Hann a 
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3. Een tabel als Tabel 3 is gemakkelijk te bepalen, wanneer het 
oorspronkelijke vraagstuk (zie $1) een maximum-probleem is, 
waarvan geen enkele der voorwaarden Ax < b een gelijkheid is 
en waarvan bovendien b >> 0. In dat geval zijn de m kolommen 
van de coëfficiëntenmatrix van de herleide vorm, welke bij de 
ingevoerde variabelen Aj ----- Am behoren, immers de eenheids- 
vectoren in R®. Zoals in $3,5 is opgemerkt, kunnen de kolommen 
van A en b onmiddellijk als lineaire combinaties hiervan worden 
geschreven. Omdat b > 0, zijn de coëfficiënten in de kolom onder b 
bovendien niet negatief. Als basisvectoren kan men dus de een- 
heidsvectoren kiezen; de verschilvariabelen (zie 81,5) zijn de 
basisvariabelen. 


Voorbeeld: Gegeven is het vraagstuk 


x1 + 2xa + 343 <6 
2x + 32 H+ 43 <6 
3x + Ka + 223 S6 
X1, %2, %3 > 0 
Xi + A2 + K3 = Max. 


De herleide vorm hiervan is 


%1 + 22 + 3x Ka =6 
21 + 32 d Kz + + % =6 
31 + Ka + 2%3 + xe = 6 
Xi, X3, X3, X4, X5, x6 =O 


Xi + Kad 43 + O.44 + O.x5 + 0O.x6 = Max. 


Het coëfficiëntenstelsel van de lineaire combinaties der eenheids- 
vectoren met als uitkomsten de verschillende kolom-vectoren in 
bovenstaande gelijkheden, is dan: 


ad) a) a) a® a5) a b 
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4. De in 3 beschreven methode ter bepaling van een tabel als 
Tabel 3 kan in de praktijk vaak worden toegepast bij de oplossing 
van een maximum-vraagstuk. Vrijwel alle in de economie be- 
schouwde grootheden — dus ook de coëfficiënten van b — nemen 
alleen positieve waarden aan. Voorts wordt veel vaker de eis 
gesteld, dat een bepaalde som kleiner dan of gelijk aan een zekere 
waarde is, dan dat gelijkheid wordt geëist. 

Bij de oplossing van een minimum-vraagstuk is de onder 3 be- 
schreven methode evenwel niet toepasbaar. De herleide vorm van 
Ax =b, x >0 is immers [A — IJ[xA] = b [A] >= 0. De kolom- 
men van — IF kunnen niet (zoals de kolommen van I) als basis- 
vectoren in Tabel 3 worden gebruikt. Daar de coëfficiënten in de 
lineaire combinatie van de kolommen van — I met als uitkomst b 
gelijk zijn aan de elementen van — b, zijn deze coëfficiënten name- 
lijk niet-negatief wanneer b > 0. 

In 18 zal worden nagegaan, hoe in andere gevallen dan het 
in 3 genoemde een tabel als Tabel 3 kan worden bepaald. 
Voorlopig wordt een dergelijke tabel, zoals gezegd, als gegeven 
beschouwd. 


5. De elementen in de laatste kolom van Tabel 3 vormen m ele- 
menten van een bruikbare oplossing x* van Ax = b, waarvan de 
overige (n — m) elementen gelijk aan O zijn. Daar hoogstens de 
m weergegeven elementen xj positief zijn, is «* een extreem punt 
van de convexe verzameling V van bruikbare oplossingen van het 
beschouwde vraagstuk (zie 82,4). De bij deze oplossing 


behorende waarde van de funktie c'« is gelijk aan 


m 
2 Ot. 
i=l 


6. Vervangt men een vector a® in de basis van Tabel 3 door een 
ander, bijv. a®, volgens de in $3 besproken techniek, dan ver- 
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krijgt men een nieuwe tabel met als waarden van de basisvaria- 
belen in de laatste kolom: 





De elementen van deze kolom vormen m elementen van een nieuwe 
basisoplossing «°, waarvan de overige (n — m) eveneens gelijk aan 
O zijn. Rangschikt men de elementen van deze nieuwe basisop- 
lossing naar volgorde van de bijbehorende kolommen van A, dan 
verkrijgt men 





18 ms 
x= Pl GER 0 x5 RE Eine ant 0 


7. Stelling: De nieuwe basisoplossing 40 is bruikbaar (d.w.z. 
* 











x0 > 0) wanneer —- het minimum is van de quotiënten. uit de 
Uhs 
“ * 
x 
reeks waarvan mis > 0. 
Mis Mms 


ak 
Bewijs: Wanneer uns > 0 is Ì_ >0 want xj is als het h-de 
â Khs 
element van x* per definitie niet-negatief. Dus: het s-de element 
van #0 is niet-negatief. 
: *® 
Voor de quotiënten —_- met mis > 0, die dus eveneens niet- 
pC His 
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* * 
Xx X 

h En ï 
Uhs Mis 


‚ hetgeen na vermenigvuldiging 








negatief zijn, geldt 











met mis oplevert: x} — En xj} =0. De uitdrukking voor het 
Uhs 
ongelijkteken is het z-de element van «°. 
xr xj Rd 
Voor de quotiënten —- met mjs <0 geldt —&- >-Z, Dus 
Hijs Hhs His _ 
eveneens is x} — Lie pe > 0. 
Hhs 


Ed 





8. Het is denkbaar, dat het minimum van de quotiënten 
‘ Mis 


bij een gegeven s voor verschillende waarden van # optreedt. 
Zoals in $5 blijkt, kan dit belangrijke consequenties hebben. 
Voorlopig wordt dit echter buiten beschouwing gelaten. Er wordt 


aangenomen, dat een willekeurige van de vectoren, waarvoor 
pn 


min 





optreedt, wordt gekozen om door a@® uit de basis te 
Mis 
worden verdreven. 


9. De bij de nieuwe oplossing «0 behorende waarde van de maxi- 
merings-, resp. minimeringsfunktie is gelijk aan 














se 
m 
* His 4 h 
= 5 ef — zi) + Cs 
i=1 Mhs Mhs 
m m se m 
= NC% + Cs — Xi Ciis 
i=1 Mhs i=1 
* 
% m 
=C'x* + L (e — DN cuz) 
Uhs i=l 


In deze uitdrukking is X%, ciuis de funktiewaarde van de ge- 
wichten uit de s-de kolom van Tabel 3. 


« 
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Kortheidshalve kan men schrijven 





mm 
XY Cikis — Ze 
i=1 
Dan is 
r , Xh 
c'x0 = c's* + (cs — Ze). 
Hhs 


10. Uit 9 blijkt onmiddellijk: 





* 
x 
Wanneer men a% door a@© vervangt en uns > o(aus Sp o) 
Hhs 


dan is €'x° > c'x* wanneer en alleen wanneer cs — 2; > 0. 
8 


11. De grootheden 2; = 2%, ciuij en dus de verschillen c; — 24 
kunnen voor alle kolommen j= 1": n van Tabel 3 worden 
berekend. 


Voorbeeld: Hieronder zijn de verschillen c‚ — zj berekend voor de 
in 3 gegeven tabel. Op de eerste regel van het onderstaande voor- 
beeld zijn de coëfficiënten c‚ van de lineaire funktie c'x weer- 
gegeven. Hieronder is de tabel uit 3 nog eens afgedrukt. Op de 
laatste twee regels zijn resp. Zi **"** ze en (c1 — 21)" ""** (ce — 26) 
vermeld. 


Cy 1 1 1 0 0) 0 





| a a) a) a a(5) a) | b 


0 AA 
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12. Stelling: Indien alle c; — 24, <0 (j= 1 --: n) in het maxi- 
mumvraagstuk (resp. indien alle c;— 24 > 0 in het minimum- 
vraagstuk) dan is de in Tabel 3 weergegeven oplossing «* een 
optimale oplossing. 

Bewijs: Deze stelling wordt alleen voor het maximum-vraagstuk 
bewezen. Het bewijs voor het minimum-vraagstuk volgt direct 
uit het onderstaande. 

Laat « in R” een willekeurige bruikbare oplossing van A« = b 
en % > 0 zijn. De funktiewaarde hiervan is Cc’. 

Bewezen moet nu worden, dat uit c; — 24, > 0 voorj=1l"--:-- ” 
in Tabel 3 volgt c'x < c'x*. 


Nn 
Ax = 5 xjAd — 
j=l 


En Ea) = 


j=1 i=1 


m gn 
= X ( zg Ja =— b. 
j=l 


i=l \j= 





Aangezien ook 


m 
Ed 
As*= 5 x;40) = b 


i=1 
en aangezien elke vector in R”, dus ook b, slechts op één manier 
als een lineaire combinatie van de basisvectoren ad... alm) 
is te schrijven (zie II 8 4, 17) geldt 
n 
Byng, Et Madi, m. 
j=1 
Uit 24 <0 (jz: n) volgt 
n n 
cx = Y Cj% SD 24% 
j=1 j=l 
n m 
ap > (2 page) 
=1 \s=1 
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Substitueert men in de laatste uitdrukking x} voor 2%, wij%, 
dan vindt men 


m 
« SN Mit 
i=1 


* 


IA 
a, 
& 


13. Stelling: Indien er in het maximum-vraagstuk een verschil 
Cs — 28 > O is (resp. indien er in het minimum-vraagstuk een 
verschil cs — ze < Ois), terwijl alle gewichten ujs <O(4=1-:--" m), 
dan heeft c'x geen eindig maximum (resp. minimum). Een optimale 
oplossing is 


Bewijs: Evenals in 9 kan er hier mede worden volstaan, het bewijs 
voor het maximum-vraagstuk te geven. Het bewijs voor het 
minimum-vraagstuk volgt rechtstreeks uit het onderstaande. 
Uit de in Tabel 3 gegeven oplossing £* kan men op (n — m) ma- 
nieren een oneindig aantal nieuwe oplossingen met (m + 1) posi- 
tieve coördinaten afleiden, en wel oplossingen van de vorm 


x= [(x{ — Aue) veen (xn — Apame) Oes Aere 0), 
waarbij O < A < oo. 
Deze vectoren zijn niet alle extreme punten van V. Wanneer alle 
gewichten mis onder a@ in Tabel 3 niet-positief zijn, zijn de op- 
lossingen x wel alle bruikbaar. In de eerste plaats is immers % > 0. 
In de tweede plaats is 


Ax = 5 (#7 — Auys)a® + Aa) 


m m 
=S xj a® —_ÂX uisa® + Aas) 
Û i=l 


b— Aa® + Ja) — b. 
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De funktiewaarde, behorende bij #, is 


C's = YS cilx; — Apis) + CsÂ 


i=l 


m m 
= XY cx) — Â XY Cipsis + CeÂ 
i=l 


i=l 


= C'X* + Alcs — Zo). 
Daar (cs — zs) > O geldt 
€'& — oo, indien Â — oo. 


In het gestelde geval kan dus steeds een oplossing worden gevonden, 
waarvan de funktiewaarde tot oneindig nadert. 


14. Uitgaande van Tabel 3 kan men nu de volgende bewerkingen 
uitvoeren: 

a) Men bepaalt voor elke kolom j=1----- n de verschillen 
(c; — 24). Is elke (c; — 24) <O in het maximum-vraagstuk (resp. 
is elke (c; — z2,) > O in het minimim-vraagstuk) dan is de gegeven 
oplossing «* volgens 12 optimaal. Het vraagstuk is hiermede 
opgelost. 

b) Wanneer niet elke (c; — 2,) <0 (resp. > 0), bepaalt men het 
maximale (resp. minimale) van deze verschillen. Zijn verschillende 
(c; — 2) maximaal (resp. minimaal), dan kiest men hiervan een 
willekeurige. Stel, (cs — 29) is het maximale (minimale) verschil 
of één van de maximale (minimale) verschillen. 

c) Vervolgens beschouwt men de gewichten wis in de kolom onder 
a®. Wanneer elk gewicht mis <0, bestaat er volgens 13 geen 
eindig optimum van het betrokken vraagstuk. De optimale 
oplossing is gelijk aan Ax*‚A —oo. Hiermede is het vraagstuk 
opgelost. 


d) Zijn niet alle gewichten wis <0, dan berekent men de posi- 
* 
tieve quotiënten — 





. Hiervan bepaalt men het minimum. Stel, 


is 
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ak 





5 

dit is &-. Voorlopig wordt aangenomen, dat het minimum slechts 
Uhs 

voor één waarde A van 2= 1 +………* m optreedt. 


e) Men vervangt a@) door a@® in de basis volgens de in $3 be- 
sproken techniek. Men verkrijgt dan een nieuwe tabel met ge- 
wijzigde coëfficiënten piy en een nieuwe bruikbare oplossing «0. 


Voorbeeld: In de onder 11 weergegeven tabel is nog geen optimale 
oplossing gevonden, want (ci — zi), (ca — 22) en (ca — 23) > 0. 
Daar (ci — zi) = (ca — 22) — (cg — 23) = 1 — max, kiest men hier- 
van een willekeurige. Stel, dat het eerste verschil gekozen wordt, 


d.w.z. dat besloten wordt ad in de basis te brengen. De quotiën- 
Ed 





ten zijn gelijk aan: 
is 
s=4 £=6 
s=5 $=3 
s=6 $=2. 


Het quotient is minimaal voor s —= 6. Derhalve moet a® uit de 
basis worden verwijderd. Als nieuwe tabel verkrijgt men dan: 


a a) a a® a(5) a(6) 


b 
a® 4 
a) 2 
an 2 





15. Volgens de in 14 beschreven methode kan men dus, uitgaande 
van een gegeven tabel, hetzij bewijzen dat de gegeven oplossing 
x* optimaal is (met een eindig of een oneindig optimum), hetzij 
een nieuwe tabel berekenen met een nieuwe oplossing «® zodanig 
dat c'x0 > c's* (resp. c'x? < c's* in het minimum vraagstuk). 
Op de tweede tabel kan Ímen de in 14 beschreven methode 
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opnieuw toepassen. Blijkt x® niet optimaal te zijn, dan kan men 
wederom een nieuwe tabel berekenen met een funktiewaarde 
xt zodanig dat c'4* > c'sl (resp. c'x+ < c'x?). Vervolgens past 
men de in 14 genoemde techniek opnieuw toe, enz. De bere- 
keningen worden beëindigd, wanneer een optimale oplossing is 
gevonden. 

Het herhaaldelijk toepassen van de in 14 genoemde techniek 
totdat een optimale oplossing is gevonden, noemt men de simplex- 
methode ter oplossing van lineaire programmeringsvraagstukken. 
Er geldt: 


16. Stelling: De simplex-methode leidt in een eindig aantal her- 
halingen tot een optimale oplossing (d.w.z. in een eindig aantal 
stappen wordt de in 14a of de in 14b genoemde toestand bereikt). 
Hierbij wordt aangenomen, dat het minimum, van de quotienten 


18 
van ? optreedt (zie 8). In $5 zal deze veronderstelling worden 
losgelaten. 


oe ‚ His > O in geen enkele herhaling voor verschillende waarden 
u 


Bewijs: Volgens $ 2,8 wordt het optimum van c'x bereikt in een 

extreem punt van de verzameling V van bruikbare oplossingen 

van Ax =b en # > 0. Volgens $ 2,6 heeft V een eindig aantal 

extreme punten. Zoals in 5 van deze paragraaf is opgemerkt, zijn 

de bij toepassing van de simplex-methode gevonden punten alle 
Xi 





extreem. Daar is aangenomen dat steeds voor één waarde 


His 
van ? optreedt, verkrijgt men bij elke herhaling een andere op- 
lossing x met een grotere (resp. kleinere) waarde van c'x dan de 
voorgaande (zie 85). Zou men de extreme punten van V naar 
grootte van de bijbehorende funktiewaarden rangschikken, dan 
komt men bij elke herhaling dus minstens één plaats verder in 
de reeks. De ervaring leert, dat men zelfs in het algemeen bij toe- 
passing van de simplex-methode niet alle extreme punten achter- 
eenvolgens behoeft te doorlopen; er worden doorgaans verschil- 
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lende punten overgeslagen. In ieder geval wordt een optimale 
oplossing echter in een eindig aantal stappen verkregen. 


17. Wanneer verschillende optimale oplossingen bestaan, kunnen 
deze alle met behulp van de simplex-methode worden bepaald. 
Stel, er is een situatie bereikt, waarin alle c; — 2, < 0. (Voor het 
minimum-vraagstuk geldt mutatis mutandis hetzelfde als voor het 
maximum-vraagstuk). Volgens 12 is de gevonden oplossing x* 
dan optimaal. De verschillen (c;— z,), behorende bij de basis- 
vectoren, zullen gelijk aan O zijn. Vervanging van een a@) in de 
basis door zichzelf verandert immers niets aan de oplossing £* en 
dus niets aan de functiewaarde. Wanneer alle andere verschillen 
(c; — zj) < O, dan is «* blijkens 10 de enige optimale oplossing. 
Wanneer daarentegen enkele van de andere verschillen gelijk aan O 
zijn, dan bestaan er nog andere optimale oplossingen van het 
beschouwde vraagstuk. Deze kunnen worden gevonden door de 
vectoren a@ waarvoor (c; — zj) = O achtereenvolgens in de basis 
te brengen. Volgens $2,9 is bovendien elke convexe lineaire 
combinatie van de op deze wijze bepaalde punten een optimale 
oplossing. 


18. Het is nu mogelijk om aan te geven, op welke wijze in het 
algemeen een oorspronkelijke tabel zoals Tabel 3 kan worden 
bepaald. Daarbij wordt nog steeds aangenomen, dat zich niet het 
in 3 genoemde bijzondere geval voordoet. 

Stel, men heeft een gegeven vraagstuk herleid tot Ax = b, « > 0 
en C'x — optimaal (max of min). Wanneer een eerste tabel niet op 
de in 3 beschreven wijze kan worden gevonden, gaat men als volgt 
te werk: 

In de eerste plaats vermenigvuldigt men alle vergelijkingen in 
Ax —= b, waarin de constante bj negatief is, met — 1. Men ver- 
krijgt dan een nieuw vergelijkingenstelsel Ax =b met b > 0. 
Aan A voegt men een eenheidsmatrix I van de orde m X m toe 
en aan &% voegt men m- nieuwe variabelen 1 «*-*" ym == y' toe. 
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Vervolgens beschouwt men het volgende vraagstuk: 


(A1) js 


lik 
y 
O Ay Fe vere HO Bij f Fr tr reen + Ym = min. 


Een eerste tabel van het laatste vraagstuk kan volgens 3 worden 
gevonden door de kolommen van I als eerste basisvectoren te 
kiezen. Lost men het vraagstuk, uitgaande van deze tabel, met 
behulp van de simplex-methode op, dan bestaan twee mogelijk- 
heden. 


a) Men vindt een oplossing ; waarin y — 0. In dat geval is een 


bruikbare oplossing # en tevens een tabel met coëfficiënten wij, 
behorende bij Ax — b, x > 0 en c'« —= optimaal, gevonden. Hier- 
mede is het doel bereikt en kan men aan de oplossing van het 
laatstgenoemde vraagstuk beginnen. 


b) Er kan geen oplossing de waarin y = 0, worden gevonden. In 


dat geval is bewezen, dat Ax — b geen (bruikbare) oplossing heeft. 


19. Er rest nu nog één probleem: Er moet nog worden onderzocht 
welke consequenties kunnen voortvloeien uit de mogelijkheid dat 
% 
His 
Dit probleem wordt in $ 5 besproken. 

Ter afsluiting van deze paragraaf worden enkele voorbeelden gegeven. 


> 0 voor verschillende waarden van 4 optreedt (zie 16). 





20." In de eerste plaats wordt op blz. 86 het in 3 genoemde (en voor 
een deel reeds in 3, 11 en 14 opgeloste) vraagstuk volgens de simplex- 
methode opgelost. In elke fase van de berekeningen is een vierkantje 
% 





geplaatst rond max (c; — 24) en min 


18 
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De optimale oplossing is #1 = 42 = x3 = |. Het maximum -van 
de lineaire funktie is 1 + 1 +1 =3. 


0 0 0) 





a a) a) a a® a 
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21. Als tweede voorbeeld wordt gekozen: 


Sx1 + 1242 — 12x3 > 24 
3x1 — Ova + 1843 — 944 = —9 
4x — 6x3 + 1243 + Axa = 6 


X1, %2, X3, X4 > 0) 
7Tx1 — 18x2 + 12x3 + 8x4 —= min. 
Dit probleem wordt in de eerste plaats door invoering van een 
verschilvariabele herleid tot 


Sx1 + 1242 — 1243 — — 45 == 24 
—3x1 + 943 — 1843 + Ova = 9 
4x — 6xa + 1243 — 4x4 


I 
on 


%1, %2, X3, Xa, X5 > O 
7x1 — 18x2 + 1243 + 8x4 + 0.45 —= min. 


Teneinde een bruikbare oplossing van dit tweede vraagstuk te 
vinden, beschouwt men volgens 18 het volgende probleem: 


Sx1 + 1242 — 1243 — — _X5 + %6 — 24 
—3x1 + Ma — 1843 Ha + + X7 = 9 
4x — 6x2 + 1243 — Axa + J-xg= 6 


X1, %2, %3, KA, X5, XK, A7, A8 2 O 


O.x1 HO. + 0.43 + Ova 4 O.x5 + K6 H- 47 + Kg —= Min. 





Het laatste vraagstuk is hierna op blz. 88 volgens de simplex- 
methode opgelost. 
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e | 0 0 0 0 0 1 1 1 
ed a a a) a) a a” a® bh 

a® WW -t Ad 1 e= al 

an | —3 [9] B 4e U B 1 ale 


a®) 


4 


—Ë 


12 —4 O 0 0) 1 6 





_ O vo Ofte ofte 
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Daar alle laatstberekende (c‚— 2) >O is aldus een bruikbare 
oplossing van het tweede vraagstuk gevonden. De eerste vijf 
kolommen en de basisvariabelen onder b van de laatst gevonden 
tabel, vormen een eerste tabel voor de oplossing van het tweede 
vraagstuk. Deze is hieronder weergegeven. 


0 


an a) a) a a) 





Het blijkt nu, dat cs — 2; < O terwijl alle mis <0. Volgens 13 
is er dus geen eindig minimum, terwijl de optimale oplossing gelijk 
is aan: 

[X1 #2 %3 A4 Al = 

[0 (9+H) (74+ EA) 6 A], A — oo. 
De optimale oplossing van het allereerste, oorspronkelijke pro- 
bleem is derhalve 


[x1 A2 A3 Kal = 
Hi Be PAD B, Ane 


Het is gemakkelijk te controleren dat deze oplossing inderdaad 
aan de voorwaarden voldoet. Immers 


of5) + (OH) 12 ]+ (742) —12} +6( O } (24+ 
3 8 18 zi eci 
4 == 12 4 6 


terwijl Â — + oo. 
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Het minimum van de lineaire funktie is: 


7.0 — 18(9 + $A) + 12(7 + A) +8.6= 
== — 30 — 2Â — — oo, daar À — + oo. 


$ 5. Degeneratie 


1. In $4, 16 is de fundamentele stelling, dat de simplexmethode 
in een eindig aantal herhalingen tot een oplossing leidt, bewezen 
Xi 





onder de veronderstelling, dat min ‚ His > O voor slechts één 


i8 : 

waarde van # optreedt. Laat men deze veronderstelling los, dan 
wordt het in $ 4, 16 gegeven bewijs aangetast. 
Teneinde dit te onderzoeken, wordt in het onderstaande in de 

5 … % 
eerste plaats nagegaan, wat er kan gebeuren, indien min kat 

Mis 

His > O voor verschillende waarden van # optreedt. 


2. Beschouw de volgende tabel. 


aD..ah) adam.  b 
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Stel, (cs — ze) = optimaal (max of min). Neem voorts aan, dat 
Xh Xk 





… @ 
ss =—= min een His > 0. Stel, men vervangt a) door a). 
Hhs Uks Mis 

Men krijgt dan een nieuwe oplossing x in de tweede tabel, waar- 


van een element, met name het /-de, gelijk aan O is. 


aD.-ah)...a)..a).. ag. Db 





3. Wanneer x® een optimale oplossing blijkt te zijn, rijzen uiteraard 
verder geen problemen. 

Wanneer dit niet het geval is, bepaalt men het optimale verschil 
(c; — zj). Stel, dit is (cr — zr). 


Xi 





Bij de bepaling van min Hir > O zijn er twee mogelijkheden: 


Mir 
a) upr < 0. In dat geval wordt een andere vector dan a® in de 
basis vervangen. De oplossing «+ in de derde tabel zal dan (aan- 
genomen dat niet weer twee of meer quotienten aan elkaar gelijk 
en minimaal zijn) geen element gelijk aan O bevatten. 


x 
b) unr > O0. Dan is min - —=0 en 





‚ Mir >O gelijk aan 





Hir Hhr 
wordt a@%) vervangen. In de nieuwe oplossing «+ is weer een ele- 
ment en wel x,; gelijk aan O. De nieuwe tabel kan in dat geval als 
volgt worden voorgesteld. 
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ad...) ab 





4. Wanneer «+ een optimale oplossing is, is het vraagstuk opgelost. 
Stel, dit is niet het geval. Het optimale verschil is (ct — zi). Men 
kan dan drie mogelijkheden onderscheiden: 

a) ure <0. Er wordt dan een andere vector dan a”) vervangen 
en men verkrijgt een vierde oplossing waarvan geen element 
gelijk aan O is. 


+ 
kp 





b) urt > O en t # h. Dan is = 0 het minimum van de quo- 


Urt 
+ 


R %i 
tienten — -, mit > O, zodat a” wordt vervangen door a® en 
Hit ' 


een nieuwe oplossing wordt verkregen waarvan een element gelijk 
aan 0 is. 

Cc) ur > O en t = h. In dit geval wordt a©) vervangen door a). 
Het laatste betekent, dat men als volgende tabel opnieuw de 
tweede (in 2 genoemde) tabel verkrijgt. Op deze wijze wordt de 
optimale oplossing dus niet bereikt. Door de simplexmethode te 
blijven toepassen voert men voortdurend dezelfde berekeningen uit; 
men blijft in een kringetje ronddraaien. Men spreekt hier van 
‘degeneratie’ van de simplexmethode. 


5. In het onder 4b genoemde geval wordt, zoals gezegd, een 
nieuwe tabel verkregen met een oplossing waarvan een element 
gelijk aan O is. Is deze oplossing niet optimaal en past men op 
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de tabel de simplexmethode toe, dan bestaan weer drie mogelijk- 
heden: 

a) er wordt een nieuwe oplossing gevonden waarvan geen enkel 
element gelijk aan 0 is. 

b) er wordt een nieuwe oplossing gevonden waarvan een element 
gelijk aan O is. 

c) er treedt degeneratie op doordat een oplossing gelijk aan &@ 
wordt verkregen. 

Wanneer b het geval is, bestaan in de volgende etappe opnieuw 
de drie bovengenoemde mogelijkheden. Het is dus denkbaar, dat 
de degeneratie zich niet onmiddellijk doch eerst na een aantal 
herhalingen voordoet. 


6. Alvorens verder te gaan moet erop worden gewezen, dat dege- 
neratie in de praktijk uiterst zelden voorkomt. Het verschijnsel, 


Xi 





dat twee of meer quotienten gelijk en minimaal zijn, komt 


His 

vrij veel voor. In het algemeen ontstaan echter geen moeilijkheden 
wanneer men hieruit een willekeurige kiest en verder de in $ 4 be- 
sproken simplexmethode zonder meer toepast. In één van de vol- 
gende etappen verkrijgt men doorgaans, hetzij direct een optimale 
oplossing met een element of meer elementen gelijk aan O, hetzij 
een oplossing > 0 (situatie a) en daarna een optimale oplossing > 0. 
Theoretisch is de degeneratie niettemin van belang. In het volgende 
deel van deze paragraaf zal daarom worden onderzocht, op welke 
wijze deze kan worden voorkomen. 


7. Gegeven is het vraagstuk I, luidende: 


n 
j=l 
xj == 0 


n 
2 C4%j = MAX. 
j=1l 


Voor het gemak wordt hier uitsluitend het maximum-vraagstuk 
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bestudeerd. De analoge resultaten voor het minimum-vraagstuk 
kunnen gemakkelijk uit het onderstaande worden afgeleid. 

Stel, A is een positief getal. Beschouw een vraagstuk II, dat als 
volgt wordt gedefinieerd: 


p) Had — p+ +E 5 wao). 


y=0 


n 


2 ejyj = max. 
j=1 


Hierin betekent Ai: ‘A tot de macht 4’. 


8. Beschouw vervolgens opnieuw Tabel 3 in $4. Uit de hierin 
gegeven oplossing #* = [xj «……*- Km Oers 0] van vraagstuk I 
kan een oplossing van vraagstuk II worden afgeleid. Een oplos- 
sing van II is namelijk: 


y= (er se 5 Hu) eneen (25 + E Hum) On 0. 
j=l j=1 
Immers: 
m 
Za 
i= 
Mm 
E (zr +E > Hug a =b+EHE ad 
i=l j=l í=l 


j=1 


Voorts zijn de coëfficiënten ui, behorende bij vraagstuk I en «*, 
uiteraard gelijk aan de coëfficiënten u4j, behorende bij vraagstuk 
HI en y*. 


9. Heeft men een optimale oplossing y® van II gevonden, dan 
__ kan hieruit een optimale oplossing «® van I worden afgeleid. De 
optimale oplossing van I is gelijk aan: 


x= (5 — E a) eaten ( tE zn) 
j=1 j=l 
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Volgens 10 voldoet deze oplossing x® aan Ax0 = b. Daar y® opti- 
maal is, is c'y® > c'y voor elke y die voldoet aan 


Ay = + a> 5 Hao) 


Derhalve is c 0 > €'& voor elke # die voldoet aan Ax — b. 


10. Stel, in Tabel 3 van 84 is cs — 2; = max en 








Xh Xk …_ Xt 
zn = min ‚ His > 0. 
Hhs Hks Mis 


Er bestaat dan een (zeer geringe) kans, dat in de etappen na de 
volgende degeneratie optreedt. Teneinde deze kans uit te sluiten 
zou men als volgt te werk kunnen gaan. 

Vraagstuk 1 vervangt men door een vraagstuk van de vorm van ÏÌ. 
De waarde van Â moet nog worden bepaald. 

In het nieuwe vraagstuk blijft c‚— z; = max. De quotienten, 
waardoor wordt bepaald, welke vector in de basis moet worden 
vervangen door a(®, zijn resp. 


n 
“+ XE Huy 
j=1 


Mls 











5 
Xn + DS Hung 
j=l 


Khs 


xe + X Mur 
j=l 
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In de eerste plaats kan de nog te bepalen constante A nu zo klein 
worden gekozen, dat het A-de en het k-de quotient in de laatste 
reeks kleiner zijn dan de andere positieve quotienten hieruit. Dit 
is in elk geval mogelijk, omdat 5, Muij — 0, indien A —0. 
Van het A-de en het k-de quotient moet vervolgens het kleinste 
had * 
worden bepaald. Daar EEEN 
Khs MKks 
vergelijking buiten beschouwing worden gelaten. Er moet dus het 
minimum worden bepaald van de volgende twee sommen: 





‚ kunnen deze termen bij de 


Tr zn En 
Mhs Khs Uhs 
en 
Rasen) + A2 ed dn anne An men 
Uks Mks Hks 
Van deze sommen kunnen de respectievelijke j-de termen niet 
aan elkaar gelijk zijn. Daar ad ----- am) basisvectoren zijn, 
kan elke vector ad ----. a immers slechts op één wijze als 


een lineaire combinatie hiervan worden geschreven. 
Stel, de eerste verschillende termen zijn 


RER on A 
Hhs Uks 
en neem aan, dat de eerste kleiner is dan de tweede. De sommen 
van de eerste 7 — 1 termen van de twee reeksen zijn dan aan 
elkaar gelijk, ongeacht de waarde van A. Men kan verder de waarde 
van À altijd zo klein kiezen, dat het verschil tussen de twee i-de 
termen 


1 (ee _ Ht ) 
Uhs Uks 


groter is dan het verschil tussen de sommen van de laatste n — í 
termen 


5 zn ee, Be 


jei+1 


hs Mks 
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Alle machten van À in het laatste verschil zijn immers groter 
dan 7. 
Kiest men voor À een dergelijke lage waarde, dan is tevens 


n Nn 
at SMumg Ht ZE Nur 
j=1 j=1 
nd en 
Uhs Mks 
Dit betekent, dat in vraagstuk II geen degeneratie optreedt, wan- 
neer men 4% uit de basis verwijdert. Aangezien echter de coëffi- 
ciënten uij in beide problemen I en II steeds dezelfde zijn, zal ook 
in probleem I geen degeneratie optreden wanneer men a uit 
de basis verwijdert. 


Il. Uit 10 volgt, dat men steeds de uit de basis te ver- 
wijderen vector zodanig kan kiezen, dat in de volgende etappen 
geen degeneratie optreedt. Het blijkt, dat degeneratie vermeden 
kan worden door de in $4, 12d gegeven regel als volgt uit te 
breiden: 


… % 8 
Wanneer min —- voor twee of meer waarden van 7 (stel 4 en k) 
Ksi 


optreedt, dan deelt men de A-de rij van de coëfficiënten ui door 
Uns en de k-de rij door prs. Men vergelijkt dan de quotienten 


DN mn EE ‚jl n twee aan twee. Is EL < LEL dan 


Hhs Hks Mhs Uks 
wordt a@® uit de basis verwijdert; in het tegengestelde geval ver- 


wijdert men a@%) uit de basis. Is nl ‚ maar Ze < Es 


Uhs Uks Hhs Uks 
dan verwijdert men eveneens a®) uit de basis, enz. 





12. Met het bovenstaande is bewezen dat elk willekeurig lineair 
programmeringsprobleem met behulp van de simplexmethode kan 
worden opgelost. In het algemeen kan met de in $ 4 behandelde 
techniek worden volstaan. Zijn verschillende van de quotienten 


t 
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TE His > O minimaal dan kiest men gewoonlijk een willekeurige 
Mis 


van de hierbij behorende vectoren ter verwijdering uit de’ basis. 
Treedt degeneratie op, dan kan men een andere kiezen. Er is altijd 
één waarvoor geen degeneratie ontstaat. Ter vermijding van de- 
generatie kan men eveneens de in 11 beschreven selectiemethode 


toepassen. 














IV. De toepassing in het bedrijf 


In dit hoofdstuk wordt getracht een beeld te geven van een aantal 
mogelijkheden om de in het voorgaande behandelde lineaire pro- 
grammeringstechniek in een bedrijf toe te passen. Hiertoe wordt 
in de eerste plaats in $ len $ 2 nagegaan op welke wijze productie 
en verkoop redelijk nauwkeurig kunnen worden beschreven door 
een vergelijkingenstelsel (een ‘model’ in de economische termino- 
logie). Vervolgens worden in $ 3 enkele opmerkingen gemaakt over 
de wijze waarop het verkoop- en productieprogramma dat de 
grootste winst oplevert, met behulp van een dergelijk vergelijkingen- 
stelsel kan worden bepaald. Hierbij blijkt men gebruik te kunnen 
maken van de in III $ 4 besproken simplexmethode. In $ 4 wordt 
in het bijzonder ingegaan op de betekenis van de lineaire program- 
mering voor transportplanning. De vergelijkingen die het transport- 
vraagstuk beschrijven hebben een bijzondere vorm waardoor een 
eenvoudiger oplossingstechniek dan de simplexmethode kan wor- 
den toegepast. Deze techniek wordt behandeld. Tenslotte worden 
in $5 nog enkele opmerkingen gemaakt over een bijzondere eigen- 
schap van de lineaire Programmering, waarvan bij de planning 
gebruik kan worden gemaakt. 


$ 1. Het produktiemodel 


1. Elk bedrijf van enige omvang kan worden verdeeld in een aantal 
afdelingen, die elk een speciale groep van min of meer homogene 
prestaties kunnen leveren. Een afdeling bestaat uit een aantal 
mensen en een aantal productiemiddelen. Ook binnen één afdeling 
kunnen vaak groepen van werknemers of series machines worden 
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onderscheiden, die specifieke prestaties verrichten. Dit wordt hier 
evenwel buiten beschouwing gelaten; de totaal-prestatie van de 
onderneming wordt niet verder gesplitst dan in de prestaties der 
afdelingen. 


2. Zowel tussen de afdelingen onderling als tussen de personen 
binnen een afdeling bestaat een zekere hiërarchie: er is een geheel 
van voorschriften omtrent de bevoegdheden der verschillende 
onderdelen en omtrent de afdeling, resp. de persoon aan wie men 
verantwoording moet afleggen. De uiteindelijke verantwoordelijk- 
heid voor de gesties van het bedrijf draagt de directie tegenover 
de eigenaars. 

De opbouw van een onderneming in afdelingen en de regeling van 
de bevoegdheden en verantwoordelijkheden worden in het alge- 
meen neergelegd in een organisatieschema. Dit schema kan als 
uitgangspunt voor de bepaling van een planningtechniek worden 
gebruikt. 


3. De afdelingen kunnen in het algemeen in de volgende categorieën 
worden verdeeld: 

— Beheerafdelingen. 

— Productieafdelingen. 

— Commerciële afdelingen. 

— Hulpafdelingen. 

In deze paragraaf zal voornamelijk de tweede groep worden be- 
schouwd. 


Voorbeeld: Een machinefabriek zou uit de volgende afdelingen 
kunnen zijn opgebouwd: 


Beheerafdelingen: 
Directie. 
Bedrijfsbureau. 
Personeelszaken. 

_ Financiële zaken 
Administratie. 
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Productieafdelingen: 
Grondstoffenopslag. 
Draaierij. 
Halffabrikatenafdeling I 
Halffabrikatenafdeling U. 
Montageafdeling 1. 
Montageafdeling II. 
Spuiterij. 
Eindproductenopslag. 
Verpakking en verlading van eindproducten. 


Commerciële afdelingen: 
Inkoop. 
Verkoop. 
Expeditie. 


Hulpafdelingen: 
Onderhoudswerkplaats. 
Laboratorium/kwaliteitsdienst. 
Huisvesting. 


4. Onder de capaciteit van een productieafdeling wordt verstaan 
de hoeveelheid prestaties die maximaal per tijdseenheid kan wor- 
den geleverd onder normale verhoudingen. De capaciteit is dus 
niet gedefinieerd als een topprestatie, doch als het maximum, dat 
gedurende lange tijd is vol te houden. 


9. Vanzelfsprekend is de capaciteit van een afdeling en dus ook 
de capaciteit van het gehele bedrijf niet onbeperkt een gegeven. 
Het is mogelijk de capaciteit uit te breiden (en ook in te krimpen, 
maar dit wordt hier buiten beschouwing gelaten). Voor het uit- 
breiden van de productiecapaciteit is echter een zekere tijd vereist. 
Men moet een uitbreidingsplan maken, geldmiddelen aantrekken, 
overleg plegen en onderhandelen met leveranciers van machines, 
personeel werven, een nieuwe organisatie creëren, enz. 

Vergroting van de capaciteit van een afdeling kan dus alleen schoks- 
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gewijze plaats vinden. Tussen elk tweetal opeenvolgende uitbrei- 
dingen ligt een zekere periode waarin de capaciteit als constant 
kan worden beschouwd; de lengte van deze periode varieert uiter- 
aard met de grootte van de uitbreiding. 

Hetzelfde geldt voor het bedrijf als geheel. Ook de capaciteiten van 
alle afdelingen tezamen zullen-gedurende zekere perioden constant 
zijn. De bestaansduur van een bedrijf kan derhalve worden verdeeld 
in perioden met gegeven capaciteiten van alle afdelingen. Ziet men af 
van kleine uitbreidingen, dan is de lengte van deze perioden in een 
normaal industrieel bedrijf doorgaans minstens gelijk aan een jaar. 


6. In deze paragraaf wordt nu de aandacht beperkt tot één van die 
perioden waarbinnen de capaciteiten der afdelingen gegeven en 
onveranderlijk zijn. Nagegaan zal worden, hoe de technische moge- 
lijkheden en productiekosten binnen een dergelijke periode kunnen 
worden beschreven. 


7. Aan het ter beschikking hebben en houden van productie- 
capaciteit zijn kosten verbonden die moeten worden gemaakt wan- 
neer wel en ook wanneer niet wordt gefabriceerd. Deze kosten die 
niet variëren met de productietijd en het tempo van de productie 
(dus bij een gegeven efficiëntie met de productieomvang), noemt 
men vaste kosten. Enerzijds omvatten de vaste kosten bepaalde 
uitgaven (‘out of pocket expenses’) zoals de huur van gebouwen, 
te betalen interest en het loon van personeelsleden die in dienst 
blijven ongeacht de omvang van het te verrichten werk. Anderzijds 
rekent men tot de vaste kosten ook de afschrijvingen op gebouwen 
en installaties welker gebruiksduur niet door slijtage maar door 
veroudering wordt bepaald. In feite zijn deze afschrijvingen reser- 
veringen voor toekomstige aanschaf van nieuwe productiemiddelen 
_ (waarvan de hoogte uit de aanschaffingsprijs of de vervangings- 

waarde van de bestaande installaties wordt afgeleid). Het heeft 
echter voor verschillende doeleinden zin deze bedragen als ‘kosten’ 
te beschouwen; omdat het bedrijf op den duur niet zou kunnen 
worden voortgezet wanneer het verschil tussen de werkelijke ont- 
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vangsten en de ‘out of pocket expenses’ op lange termijn niet 
minstens gelijk aan de afschrijvingen zou zijn. 


8. De vaste kosten zijn gegeven als een totaal bedrag per beschouwde 
periode. Hiermede wil niet gezegd zijn, dat de vaste kosten per tijds- 
eenheid binnen de beschouwde periode niet kunnen variëren. Men kan 
huur opzeggen, arbeiders ontslaan en machines verkopen. Voorts is 
het mogelijk dat een wijziging in de lonen, de huur, de waarde der 
machines, enz. optreedt. Deze veranderingen worden evenwel buiten 
beschouwing gelaten. Aangenomen wordt, dat de ‘vaste’ kosten 
binnen de in 6 bedoelde periode in het geheel niet veranderen. 


9. De technische mogelijkheden van een bedrijf worden in de eerste 
plaats bepaald door de capaciteiten der productieafdelingen. In 
sommige gevallen kunnen ook de opslagruimten een beperkende 
factor voor de productie zijn. Eenvoudigheidshalve wordt hier 
evenwel van het laatste afgezien. Zijn de beschikbare opslag- 
loodsen een ‘bottle neck’ voor de productie, dan geldt ten aanzien 
hiervan mutatis mutandis hetzelfde als hierna ten aanzien van de 
capaciteiten der eigenlijke fabrieksafdelingen wordt opgemerkt. 


Voorbeeld: Van de onder 3 genoemde machinefabriek zijn de vol- 
gende afdelingen bepalend voor de productiemogelijkheden. 


Halffabr. afd. | Halffabr. afd. II 
Montageafd. | Montageafd. Il Montageafd. II 
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10. Om een eenheid van een zeker product te fabriceren moeten 
in het algemeen verschillende productieafdelingen prestaties leve- 
ren. Het totaal der opeenvolgende bewerkingen, vereist om een 
eenheid van een zeker product voort te brengen, wordt een pro- 
ductieactiviteit genoemd. 


11. In sommige gevallen gaat de productie van een zeker goed om 
technische redenen altijd gepaard met de productie van een ander 
goed of een aantal andere goederen (bijv. bij de graanbouw ver- 
krijgt men graan en stro, bij een hoogovenproces verkrijgt men ruw 
ijzer, hoogovengas en hoogovenslak). Men spreekt dan van gemeen- 
schappelijke productie. Vaak kan men hierbij één van de pro- 
ducten als hoofdproduct en de andere(n) als bijproduct(en) be- 
stempelen; dit behoeft echter niet het geval te zijn. 

Het is eenvoudig het begrip productieactiviteit ook in het geval van 
gemeenschappelijke productie toe te passen. Daarbij kan onder- 
scheid worden gemaakt tussen het geval waarin de verschillende 
producten slechts in één vaste verhouding kunnen worden gefabri- 
ceerd, en het geval waarin de verhouding binnen zekere grenzen 
varieerbaar is. 

Als er slechts één vaste productieverhouding bestaat, wordt een 
productieactiviteit hiervoor ingevoerd. Hiertoe definieert men een 
‘eenheidspakket’ van de betrokken producten. In principe kan elke 
combinatie van hoeveelheden, die voldoen aan de vaste productie- 
verhouding, als ‘eenheidspakket’ worden beschouwd. Doorgaans 
zal echter die combinatie als ‘eenheidspakket’ worden gekozen, 
waarin de hoeveelheid van het belangrijkste goed gelijk aan één 
eenheid is. De productieactiviteit kan dan worden gedefinieerd als 
het totaal der opeenvolgende bewerkingen, vereist om het ‘een- 
heidspakket’ voort te brengen. (In de terminologie van 10 wordt 
dus ‘eenheid product’ vervangen door ‘eenheidspakket’). 

Als de productieverhouding binnen zekere grenzen kan worden 
gevarieerd, kiest men uit het gehele scala van mogelijkheden een 
aantal verhoudingen die het totaal goed representeren (d.w.z. die 
regelmatig over het variatie-interval zijn verdeeld). De gekozen 















Het produktiemodel 


verhoudingen worden beschouwd als beschrijving van de gehele 
verzameling van mogelijkheden. Door het aantal beschouwde ver- 
houdingen te vergroten verkrijgt men uiteraard een beter beeld 
van de werkelijkheid. Voor elke beschouwde verhouding definieert 
men nu een ‘eenheidspakket’. Per verhouding verkrijgt men dan 
een productieactiviteit als totaal van de opeenvolgende bewerkin- 
gen, vereist om het bijbehorende ‘eenheidspakket’ te fabriceren. 


12. Als een bedrijf uit » productieafdelingen bestaat, kan een 
productieactiviteit worden voorgesteld door een kolomvector van 
m productietijden. De j-de van de onderscheiden productieactivi- 
teiten kan bijv. worden aangegeven door 


an be Gij 


Ami 
Elk element aj; van deze vector is gelijk aan het tijdsbeslag, dat 
voor de uitvoering van de j-de productieactiviteit wordt gelegd op 
de ú-de afdeling. Aangenomen wordt, dat de productieafdelingen 


in een vaste volgorde zijn gerangschikt. Uiteraard kunnen verschil- 
lende elementen a; gelijk aan O zijn. 


13. De tijdsbeslagen 44; worden aan de administratie ontleend. 
Meestal worden de grootheden aj; zelf niet geadministreerd. Men 
moet ze dan uit andere gegevens berekenen. Op welke wijze dit 
kan geschieden, kan het gemakkelijkst aan de hand van een voor- 
beeld worden toegelicht. 


Voorbeeld: Stel, dat een eenheid van een zeker product A op de 
volgende wijze in de onder 3 en 4 genoemde machinefabriek kan 
worden gefabriceerd: 

In de draaierij worden r assen type R en s assen type S gedraaid. 
In halffabrikatenafdeling I worden u eenheden van goed U en in 
halffabrikatenafdeling II worden v eenheden van goed V voort- 
gebracht. De assen en de halffabrikaten worden vervolgens in 
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montageafdeling II tot een eenheid van het beschouwde eindpro- 
duct samengesteld. Tenslotte wordt het product A in de spuiterij 
gespoten. 

Aan de administratie kunnen in het algemeen de volgende gegevens 
worden ontleend. 





draaitijd per as type R in de draaierij br 
draaitijd per as type S in de draaierij bs 
productietijd per eenheid van goed U in halffabrikatenafd. I | tu 
productietijd per eenheid van goed V in halffabrikatenafd. II | 4 
productietijd per eenheid van goed A in montageafd. IT ta 
spuittijd per eenheid van goed A in de spuiterij te 





Op grond hiervan kan de productieactiviteit a®, behorende bij 
de beschouwde productiewijze van goed A als volgt worden be- 
rekend: 


Draaierij ay = [7 Xr Hs X ts 
Halffabrikatenafd. I 4aj ux ty 
Halffabrikatenafd. II | a3; 

Montageafdeling 1 aas 

Montageafdeling II a5j 

Montageafdeling III aes 

Spuiterij ay 


14. Het aantal te onderscheiden productieactiviteiten is minstens 
even groot als het aantal producten (of vaste combinaties van 
producten), dat het bedrijf kan voortbrengen. (Teneinde omslach- 
tige formuleringen te voorkomen, wordt hierna niet gesproken van 
‘de vector ad die de j-de productieactiviteit aangeeft’, doch van 
‘de productieactiviteit a@)’). 

Wanneer de fabricages van twee verschillende producten (of vaste 
combinaties van producten) dezelfde bewerkingstijden van alle 
afdelingen zouden vereisen, zou hiervoor slechts één productie- 
activiteit behoeven te worden onderscheiden. Indien daarentegen 





tar wen GE EEE Sia 
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een zelfde product (of combinatie van producten) op uiteenlopende 
manieren kan worden voortgebracht, wordt voor elke productie- 
wijze een productieactiviteit ingevoerd. 


15. Boven de in 7 genoemde vaste kosten, die onafhankelijk van 
de productieomvang zijn binnen een in 5 gedefinieerde periode, 
heeft de uitvoering van een productieactiviteit bepaalde kosten tot 
gevolg, die in twee groepen zijn te onderscheiden: 

a) de variabele tijdkosten. Onder de variabele tijdkosten van een 
afdeling worden verstaan de kosten, die variëren met de productie- 
tijd, ongeacht welk product wordt voortgebracht. Tot deze kosten- 
groep behoren o.m. de kosten van het energieverbruik (stroom, 
kolen, olie), smeermiddelen, schroeven, enz. Voorts moeten hiertoe 
worden gerekend de afschrijvingen op machines waarvan de bruik- 
baarheidsduur afhankelijk is van de slijtage. 

De variabele tijdkosten van een afdeling zijn in het algemeen per 
uur gegeven. Gegevens hieromtrent moeten aan de administratie 
worden ontleend. De berekening van deze kosten vereist in de 
eerste plaats een zeer nauwkeurige analyse van de budgettarieven 
per afdeling en vervolgens een voortdurend aanpassen hiervan aan 
wijzigingen welke gedurende het boekjaar kunnen optreden. De 
methode van kostencalculatie wordt in hoge mate bepaald door 
de toegepaste administratieve techniek. Hierop wordt te dezer 
plaatse niet ingegaan. 7 

Om de variabele tijdkosten per productieactiviteit te bepalen dient 
men de tijdsbeslagen in uren op de verschillende afdelingen te ver- 
menigvuldigen met de variabele tijdkosten per uur en deze pro- 
ducten vervolgens te sommeren. Wanneer de kosten per uur een- 
maal zijn berekend, zijn hieraan geen problemen verbonden. 

b) de grondstofkosten. De grondstofkosten worden in de administra- 
tie per half- of eindproduct berekend. Hiertoe gaat men in het alge- 
meen uit van standaardverbruiken per eenheid product en vaste 
verrekenprijzen per eenheid grondstof. Aangenomen dat de stan- 
daardverbruiken overeenstemmen met de geschatte werkelijke 
verbruiken in de beschouwde, toekomstige periode kunnen deze 
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ook als uitgangspunt voor de planning worden gebruikt. De stan- 
daardverrekenprijzen zullen zelden overeenstemmen met de door 
de inkoopafdeling verwachte prijzen in de planperiode. De grond- 
stoftarieven moeten dus doorgaans opnieuw voor de planning wor- 
den berekend op basis van de laatste schattingen van de prijs- 
ontwikkeling. 

Wanneer de grondstofkosten per eenheid half- en eindproduct zijn 
gegeven, vindt men de grondstofkosten van een productieactiviteit 
door de bij uitvoering van deze activiteit verbruikte hoeveelheden 
met de eenheidsprijzen te vermenigvuldigen en de producten op 
te tellen. Bij het bepalen van de grondstofverbruiken per productie- 
activiteit hoede men zich voor dubbeltellingen. In de achtereen- 
volgende afdelingen worden grondstoffen verbruikt. Naarmate een 
product verder zijn voltooiing nadert, zijn steeds meer grondstoffen 
erin opgegaan. De waarde van de grondstoffen die voor de uit- 
voering van een basisactiviteit moeten worden verbruikt, is gelijk 
aan de toegevoegde grondstofwaarden in de achtereenvolgende 
afdelingen die het product doorloopt. 


16. De som van de variabele tijdkosten en de grondstofkosten 
per eenheid product noemt men de variabele kosten zonder meer. 
In het volgende wordt de som van de variabele tijdkosten en grond- 
stofkosten per productieactiviteit aangeduid als de variabele kos- 
ten van die activiteit. 

De variabele kosten van productieactiviteit a® worden voorgesteld. 
door Zj. 


17. Er wordt verondersteld, dat er # productieactiviteiten 
ad)... a bestaan. 
Blijkens de definitie stelt elke positieve lineaire combinatie hiervan 


Ajan L eerie sé + Aram 
AF ennen MED 0 


een productieprogramma voor. Niet elk productieprogramma zal 
evenwel binnen het raam van de beschikbare capaciteiten kunnen 
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worden gerealiseerd. Aan de waarden van de variabelen 4; zijn 
dus zekere beperkingen gesteld. Hierop wordt in het volgende 
ingegaan. Er kan nu reeds worden opgemerkt, dat de variabele 
kosten van een realiseerbaar productieprogramma gelijk zijn aan 


18. Of een gegeven productieprogramma binnen een zekere periode 
kan worden gerealiseerd is, behalve van de vereiste productietijden 
a in de basisactiviteiten, mede afhankelijk van de tijden gedurende 
welke de verschillende afdelingen voor de productie ter beschik- 
king staan. De beschikbaarheidsduur van elke afdeling vindt men 
door de kalendertijd te verminderen met de noodzakelijke stil- 
standstijden wegens schoonmaken, reparatie, onderhoud, vacantie 
van het personeel, enz. Ten dele worden de stilstandstijden bepaald 
door institutionele factoren, zoals bijv. de voorgeschreven werk- 
tijden van het personeel. Voorts beïnvloeden toevallige, en dus tot 
op zekere hoogte voorspelbare, factoren de stilstandstijden. Ten- 
slotte heeft het bedrijf zelf een zekere invloed op de stilstandstijden 
door verschuiving en/of verlenging of verkorting van de werktijden, 
de onderhoudsperioden, e.d. 


19. Hier worden de stilstandsperioden als gegeven beschouwd. 
De tijden, gedurende welke de m afdelingen binnen een zekere 
periode voor productie ter beschikking staan, kunnen dan resp. 
worden voorgesteld door de elementen van een vector b van de 
orde m X 1. 


b =| b1 


20. Het totaal van de technische mogelijkheden van een bedrijf 
in een periode waarbinnen de capaciteiten der verschillende af- 
delingen constant zijn, kan nu worden gedefinieerd als: de verza- 
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meling der vectoren [À1 *--*- An], welke voldoen aan 
had —- VERES + Ma <S 
Â1 ER Â en > le) 


De bij de verschillende productieprogramma’s behorende variabele 
kosten worden weergegeven door 


21. Zijn bepaalde waarden Aj «:**: AË gegeven, dan kunnen deze 
altijd worden ‘terugvertaald’ in hoeveelheden van de verschillende 
halffabrikaten en eindproducten. Hiertoe dient men, uitgaande 
van de definities der basisactiviteiten, een berekening uit te voeren, 
omgekeerd aan de onder 13 genoemde. 


Voorbeeld: Wanneer gegeven is, dat onder andere de in het voor- 
beeld van 12 genoemde basisactiviteit a@ A} maal moet worden 
uitgevoerd, dan moeten worden gefabriceerd: 


Af X r eenheden van as type R 
Af Xs eenheden van as type S 
Af X u eenheden van goed U 
Aj X v eenheden van goed V 
Aj eenheden van goed A, 


terwijl eveneens Aj eenheden van goed A moeten worden gespoten. 
Voert men soortgelijke berekeningen ook uit voor andere produc- 
ten Ajan, jl sense n, in een lineaire combinatie der basis- 
activiteiten en telt men de aldus gevonden hoeveelheden half- 
fabrikaten en eindproducten op, dan wordt het productieprogram- 
ma van elk der afdelingen gevonden. 


22. Tot slot kan worden opgemerkt, dat de hierboven besproken 
lineaire beschrijving van een productieproces en bijbehorende pro- 
ductiekosten op enige veronderstellingen berust. Men kan de wer- 
kelijkheid nooit volledig in een aantal formules vangen. Hoogstens 
kan men haar zodanig beschrijven, dat aan de formules zekere 
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richtlijnen voor het handelen kunnen worden ontleend. Een derge- 
lijke beschrijving kan een model worden genoemd; in het voorgaande 
is dus een productiemodel opgesteld. 

Hierbij is in de eerste plaats aangenomen, dat het mogelijk is de 
tijdsbeslagen aj, de variabele tijdkosten per uur en de variabele 
grondstofkosten per eenheid product vast te stellen. Vanzelfspre- 
kend zullen deze grootheden in achtereenvolgende herhalingen van 
een productieactiviteit niet steeds dezelfde waarden aannemen. 
Men kan nu eenmaal een bepaald werk niet tweemaal op precies 
dezelfde wijze uitvoeren. In het bovenstaande is ervan uitgegaan, 
dat men hetzij op grond van vroegere prestaties hetzij op grond 
van een analyse der te verrichten werkzaamheden ‘gemiddelde’ of 
‘normale’ tijdsbeslagen, variabele kosten per uur en grondstof- 
verbruiken kan bepalen. 

In de tweede plaats is verondersteld, dat zowel de basisactiviteiten 
als de variabele kosten optelbaar zijn. Dit houdt in, dat de productie 
van Â eenheden van alle goederen A maal zo veel tijd vergt en Â maal 
zo veel variabele kosten met zich brengt als de productie van 1 een- 
heid van alle goederen. Hiermede wordt voorbijgegaan aan het 
veelvuldig geconstateerde feit, dat grote hoeveelheden van een 
zelfde goed per eenheid sneller en met minder variabele (bijv. 
grondstof-) kosten kunnen worden gefabriceerd dan enkele stuks. 
De lineaire beschrijving van een productieproces is dus slechts een 
benadering van de werkelijkheid. Men kan zich voorstellen, dat 
een betere benadering zou worden verkregen door ingewikkelder 
functies in te voeren. In de praktijk blijkt evenwel, dat het gebruik 
van lineaire betrekkingen in het algemeen voldoende nauwkeurige 
resultaten oplevert. De gehele moderne kostenbudgettering en 
-administratie berust op de twee bovengenoemde veronderstel- 
lingen. 


$2. Het verkoopmodel 


Ll. Voor elk van de producten welke een bedrijf kan fabriceren, 
bestaan in het algemeen vele potentiële afnemers. Deze mogelijke 
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afnemers kunnen veelal in een aantal categorieën worden ingedeeld. 
De meest fundamentele indeling is wellicht die naar de gebieden 
waarin zij zijn gevestigd. Deze indeling kan meer of minder ver- 
fijnd zijn. De meest eenvoudige is die in twee gebieden, bijv.: het 
binnenland en het buitenland, de ‘eigen wijk’ en het gebied daar- 
buiten. Een meer verfijnde indeling zou die zijn waarbij de binnen- 
landse markt bijv. in een aantal rayons wordt verdeeld en het 
buitenland wordt onderscheiden in Europa, N. Amerika, Z. Ame- 
rika, N. Afrika en het Midden Oosten, de rest van Afrika, de rest 
van Azië en Australië. 

De verdeling van het totale potentiële afzetgebied in een aantal 
deelgebieden is doorgaans de basis van een analyse van de afzet- 
mogelijkheden. In de eerste plaats, omdat de vraag van de af- 
nemers naar de beschouwde artikelen van gebied tot gebied zullen 
verschillen. Voorts omdat de concurrentie in de verschillende ge- 
bieden ongelijk zal zijn en omdat de eigen verkooporganisatie 
niet overal in dezelfde mate is ingevoerd en op dezelfde wijze is 
georganiseerd. Tenslotte omdat de transportkosten van de plaats 
van productie naar de verschillende afzetgebieden zullen uiteen- 
lopen en de overheidsheffingen (invoerrechten, omzetbelasting e.d.) 
niet overal gelijk zijn. 


2. Binnen elk van de onderscheiden deelgebieden kunnen de af- 
nemers verder worden onderverdeeld. Een belangrijke onder- 
scheiding is bijv. die naar funktie in de maatschappij. Een grove 
indeling naar dit kenmerk is bijv. die in consumenten, handelaren 
en producenten. Deze zou kunnen worden verfijnd door de consu- 
menten te onderscheiden naar het beroep en de handelaren en 
producenten naar de soort verhandelde, resp. gefabriceerde goe- 
deren. 


3. Vrijwel elk kenmerk van de potentiële afnemers van een zeker 
artikel kan als indelingscriterium zinvol zijn. Het hangt af van de 
betekenis van de verschillende kenmerken voor de vraag der af- 
nemers, de beschikbare gegevens en de omvang van het onderzoek, 
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welke maatstaven in een gegeven geval zullen worden gekozen. 
Elk van de onderscheiden klassen van potentiële afnemers wordt 
aangeduid als een deelmarkt. 


4. Om te kunnen verkopen moet men beschikken over een acqui- 
sitie- en distributieapparaat. Hieronder wordt verstaan het totaal 
van de mensen, gebouwen, machines en transportmiddelen door 
middel waarvan de verkoopovereenkomsten worden afgesloten, de 
goederen van de fabriek naar de afnemers worden gebracht en de 
uit de verkoopovereenkomsten voortvloeiende verplichtingen wor- 
den uitgevoerd. Acquisitie (het aantrekken van gegadigden, de 
overreding tot de koop en het afsluiten van de overeenkomst) en. 
distributie (het afwikkelen van de afgesloten overeenkomst) zijn 
duidelijk gescheiden handelingen. In grote verkooporganisaties 
houden sommige personen zich zelfs alleen met de eerste en anderen 
alleen met de tweede funktie bezig. In het algemeen kan echter niet 
van twee afzonderlijke organisaties worden gesproken. Het acqui- 
sitie- en distributieapparaat wordt daarom in het volgende be- 
schouwd als een geheel, aangeduid met het woord verkooporga- 
nisatie. 

De verkooporganisaties van verschillende bedrijven lopen zeer 
uiteen. Een kruidenier heeft veelal als verkooporganisatie een 
winkel met een aantal personeelsleden. De verkooporganisatie van 
een grote industriële onderneming kan bijv. bestaan uit een cen- 
traal hoofdkantoor en verschillende filialen, die elk een aantal 
reizigers in dienst hebben en over verschillende transportmiddelen 
beschikken. Elke verkooporganisatie vertoont echter de overeen- 
komst met een productieapparaat, dat zij alleen met veel kosten 
kan worden opgebouwd en niet binnen korte tijd kan worden 
gewijzigd. 


5. Zoals in $ 1.5 is opgemerkt kan de levensduur van een bedrijf 
worden verdeeld in perioden waarbinnen de capaciteiten der ver- 
schillende productieafdelingen onveranderlijk zijn. Veelal zal ook 
de verkooporganisatie, zowel het acquisitie- als het distributie- 
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apparaat, binnen elk van deze perioden vast zijn, d.w.z. niet 
variëren met de omvang van de verkoop. Een verschil tussen de 
verkooporganisatie en het productieapparaat is, dat men niet 
gemakkelijk een capaciteit van de eerste kan aangeven. Weliswaar 
zal een verkoopafdeling niet meer dan een zekere acquisitieactiviteit 
kunnen ontplooien en niet meer dan een zeker aantal orders kunnen 
afwikkelen, doch de organisatorische en administratieve mogelijk- 
heden zijn in de regel groot genoeg om elk te fabriceren assortiment 
te kunnen behandelen. Hiervan zal dan ook worden uitgegaan. 


6. Uit hetgeen in 5 is opgemerkt volgt, dat verschillende acquisitie- 
en distributiekosten — zoals de lonen van een deel van het verkoop- 
personeel en de afschrijvingen op gebouwen en transportmiddelen — 
binnen een in $ 1.5 beschreven periode niet variëren met de ver- 
kochte hoeveelheden. In die zin zijn deze kosten ‘vast’. Naast de 
vaste productiekosten bestaan dus vaste verkoopkosten. 


7. Het resultaat van de ontplooide verkoopactiviteit op een deel- 
markt is, dat een zekere hoeveelheid aan de cliënten wordt ge- 
leverd en een zekere opbrengst wordt ontvangen. De verkochte 
hoeveelheid wordt de afzet op de betrokken deelmarkt genoemd. 


8. Er kan onderscheid worden gemaakt tussen de algemene prijs- 
notering en het gemiddelde bedrag dat de cliënten per eenheid 
product betalen. Het gemiddelde door de cliënten per eenheid 
product betaalde bedrag is gelijk aan de algemene prijsnotering 
minus de volgens de secundaire verkoopvoorwaarden (kortingen 
voor contante betaling, rabat, omzetpremies, commissies, e.d.) te 
verkrijgen aftrekposten. In het volgende zal onder de prijs van een 
artikel steeds worden verstaan het in feite door de cliënten betaalde 
bedrag per eenheid product. De prijs hangt dus af van de mate 
waarin de afnemers gebruik maken van de in de secundaire ver- 
koapvoorwaarden geboden mogelijkheden. 


9. Om te bepalen, wat de verkoop op een zekere deelmarkt uit- 
eindelijk oplevert, moet de bruto opbrengst (—= som van de ver- 
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kochte hoeveelheden x de door de cliënten in feite betaalde 
prijzen) in de eerste plaats verminderd worden met de kosten, ver- 
bonden aan de stimulering van de verkoop zoals reizigersbezoek, 
het plaatsen van advertenties, het uitgeven van folders en het 
deelnemen aan een jaarbeurs, voorzover deze kosten niet reeds zijn 
begrepen onder de in 6 genoemde vaste verkoopkosten. Het gaat 
hier dus om de extra kosten voor acquisitie op een bepaalde deel- 
markt boven de vaste verkoopkosten. 

Omtrent deze kosten kan een tweetal opmerkingen worden ge- 
maakt. In de eerste plaats kan men stellen dat de uitgaven voor 
acquisitie in een zekere periode niet alleen kosten van deze periode 
zijn, indien er voorzover de acquisitie ook de verkopen in latere 
perioden positief beïnvloedt. Kan het laatste met een grote mate 
van waarschijnlijkheid worden vastgesteld, dan kan worden besloten 
niet het totale bedrag als kosten van de verkopen in de beschouwde 
periode te berekenen. In de tweede plaats kan men opmerken dat 
met het noemen van de acquisitiekosten de aard van de acquisitie 
nog niet vast ligt. Bij een gegeven bedrag zijn verschillende ver- 
koopcampagnes denkbaar. In het volgende wordt er steeds van 
uitgegaan, dat dan die zal worden gekozen welke de potentiële 
kopers het sterkst beïnvloedt. 

De extra kosten van stimulering van de verkoop kan men de 
variabele acquisitiekosten noemen. Men dient dan evenwel in het oog 
te houden dat deze kosten niet variabel zijn inde zin van de varia- 
bele fabricagekosten: er bestaat geen duidelijk lineair verband 
tussen de variabele acquisitiekosten en de verkochte hoeveelheid. 
De variabele acquisitiekosten gaan aan de verkoop vooraf en het 
hangt van de markt af, hoe de afzet hierop zal reageren. Deze kosten 
zijn dan ook niet gegeven als een bedrag per eenheid product, doch 
als een totaal bedrag in de beschouwde planperiode. 


10. Een tweede categorie kosten die van de opbrengst moet wor- 
den afgetrokken zijn de variabele distributiekosten. Hieronder wordt 
verstaan die kosten verbonden aan het transport en het afleveren 
van de goederen welke direct met de verkochte hoeveelheid samen- 
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hangen. Hiertoe behoren in de eerste plaats de kosten van transport, 
overslag, opslag in loodsen van derden, kosten van formulieren, 
makelaarscourtage, enz. Deze kosten zijn dus variabel in afhan- 
kelijkheid van de verkochte hoeveelheid op dezelfde wijze als 
waarop de variabele fabricagekosten afhankelijk zijn van de gepro- 
duceerde hoeveelheid. De variabele distributiekosten zijn te bere- 
kenen als een bedrag per eenheid product. 


11. De laatste groep verkoopkosten zijn de heffingen en vestituties 
door de overheid. Bij verkoop in het binnenland behoeft in het 
algemeen alleen rekening te worden gehouden met de omzetbelas- 
ting. Bij export krijgt men doorgaans te maken met een veelheid. 
van heffingen, zoals de restitutie van de omzetbelasting op de 
verbruikte grondstoffen, invoerrechten en compenserende heffin- 
gen aan de grens. Het totaal van restituties en heffingen is variabel 
in dezelfde zin als de distributiekosten. 


12. Onder de bruto af fabriek opbrengst van de verkopen op een 
zekere deelmarkt kan worden verstaan: het in feite door de af- 
nemers betaalde bedrag minus de variabele acquisitiekosten minus 
de afzet X de variabele distributiekosten per eenheid minus de 
afzet X de som van heffingen en restituties per eenheid. De netto 
af fabriek opbrengst van de verkopen op een zekere deelmarkt kan 
worden gedefinieerd als het verschil tussen de bruto af fabriek 
opbrengst en de vaste verkoopkosten voor die deelmarkt. 


13. Op elke deelmarkt zal binnen een in $ 1,5 en in 6 van deze 
paragraaf genoemde periode een zeker verband bestaan tussen de 
te verkopen hoeveelheden enerzijds en de prijzen en variabele 
acquisitiekosten anderzijds. (De variabele acquisitiekosten zijn 
volgens 9 verbonden aan een en niet meer dan een acquisitie- 
activiteit.) Dit verband wordt aangeduid als een afzetfunktie. 

Met deze terminologie wordt niet bedoeld, dat de betrekking tussen 
de afzetten, prijzen en variabele acquisitiekosten in één vergelijking 
kan worden uitgedrukt. In haar algemene gedaante zal de ‘afzet- 
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funktie’ bestaan uit een vergelijkingenstelsel. Stel bijv, er zijn 
r producten. De te behalen afzet van product # op deelmarkt 4 
wordt voorgesteld door Vij. De prijs van product d op deelmarkt 
wordt voorgesteld door Pi. De variabele acquisitiekosten voor 
product # op deelmarkt j worden voorgesteld door Aij. In dat 
geval kan de afzetfunktie van deelmarkt 4 door de volgende relaties 
worden beschreven: 


Vij= hiy(Pig «0e Pap dag == A44) 


Voy = f(Paig «0e Poj, Ang oere A45) 


Bij deze formulering is ervan uitgegaan, dat voor elk product een 
afzonderlijke verkoopactiviteit wordt ontwikkeld. Zou een alge- 
mene reclame voor alle producten tezamen worden gemaakt, dan 
zou in plaats van of naast Aij «…-*- Ars één variabele A; moeten 
worden ingevoerd. 


14. De vorm van de afzetfunktie van een bepaalde deelmarkt 
wordt door een groot aantal factoren bepaald. In de eerste plaats 
wordt deze bepaald door de aard en de intensiteit van de behoefte 
der potentiële afnemers. De vraag van een consument is primair 
afhankelijk door fysieke, psychische en sociale factoren en daar- 
binnen door financiële, zoals het inkomen en de prijzen van andere 
goederen dan de beschouwde. De vraag vaneen bedrijf zal o.m. 
worden bepaald door de aard van zijn productieproces, de prijzen 
van andere grondstoffen dan de beschouwde en de opbrengsten 
van de eindproducten. Naarmate de prijs van een artikel lager is 
zal elke potentiële afnemer doorgaans een grotere hoeveelheid 
ervan willen aanschaffen; enerzijds omdat men bij een gegeven te 
besteden bedrag in staat is meer te kopen, anderzijds omdat men 
een goed waarvan de prijs wordt verlaagd, ter vervanging van 
andere zal gaan gebruiken (ballpoints en vulpennen, kolen en olie, 
enz.). Hetzelfde geldt t.a.v. de verkoopactiviteiten: naarmate de 
acquisitie intensiever wordt bedreven, zal elke potentiële afnemer 
meer willen gebruiken. 
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iS. Een tweede bepalende factor voor de afzetfunktie van een 
deelmarkt wordt gevormd door de concurrentieverhoudingen. Of 
een onderneming meer zal gaan verkopen door prijsverlaging en 
verhoging van de verkoopactiviteit, hangt niet alleen af van de 
reacties der verbruikers, doch ook van die der concurrenten. Als 
de concurrenten op een prijsverlaging van één van hen reageren 
door hun prijzen nog verder te verlagen, is het zelfs denkbaar, 
dat de onderneming die de prijs als eerste verlaagt, niet méér doch 
minder gaat verkopen. Hoe de concurrenten in feite zullen reageren, 
hangt ook weer van een groot aantal feiten af; o.a. de kwaliteit en de 
omvang van hun verkooporganisaties en hun beschikbare productie- 
capaciteiten kunnen hiervoor van betekenis zijn. - 


16. Als derde bepalende factor voor de afzetfunktie kan de in. het 
verleden opgebouwde marktpositie worden genoemd. Onder de 
marktpositie van een onderneming kan worden verstaan het totaal 
van de denkbeelden die de potentiële afnemers zich over de onder- 
neming hebben gevormd. In de eerste plaats behoort hiertoe de 
algemene bekendheid van de handelsnaam en de zakelijke reputatie. 
In de tweede plaats bestaat de marktpositie van een onderneming 
uit de voorkeur van een deel der afnemers om bij haar te kopen. 
Deze voorkeur van de z.g. vaste klanten kan ten dele zijn gebaseerd 
op een door geregeld contact gegroeide persoonlijke band tussen 
koper en verkoper. Er zijn echter ook zuiver financiële redenen die 
maken dat men bij voorkeur bij een vroegere leverancier koopt. 
Zoals een efficiënte fabricage in vele gevallen alleen denkbaar is, 
wanneer voortdurend handelingen moeten worden herhaald, zo is 
ook een vlotte afwikkeling van de ruiltransacties tussen verschil- 
lende huishoudingen alleen mogelijk, wanneer steeds met dezelfde 
personen of firma’s zaken worden gedaan. Als koper en verkoper 
elkaar kennen, kunnen bijv. inspectie der goederen vooraf, onder- 
handelingen over leverings- en betalingscondities, navrage omtrent 
de credietwaardigheid e.d. achterwege blijven. 

In het algemeen zal men meer op een bepaalde deelmarkt kunnen 
verkopen, naarmate men beter is ingevoerd. De feitelijke betekenis 
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van de marktpositie is echter in hoge mate afhankelijk van de 
marktorganisatie. Wanneer de koop en verkoop centraal wordt 
georganiseerd door een beurs of een veiling, is het contact tussen 
koper en verkoper doorgaans onpersoonlijk en niet-duurzaam. In 
dat geval kan soms nauwelijks van een marktpositie worden ge- 
sproken. Bij bilaterale verkoop daarentegen zijn de verkoop- 
mogelijkheden vaak sterk afhankelijk van de mate waarin het 
bedrijf is ingevoerd. Er wordt dan wel beweerd, dat ‘it is more 
important to own a market than a factory’. 


17. Er kan onderscheid worden gemaakt tussen 

a) de afzetfunktie op lange termijn, d.w.z. na een aantal, bijv. vijf 
perioden (jaren); 

b) de afzetfunktie op korte termijn, d.w.z. in de eerstvolgende 
periode (jaar). 

Beide afzetfunkties kunnen worden beschreven door een stelsel als 
in 13 is weergegeven. Tussen de waarden van de coëfficiënten in de 
afzetfunktie op lange termijn en die in de afzetfunktie op korte 
termijn kunnen echter belangrijke verschillen bestaan. 

Deze verschillen kunnen in de eerste plaats voortvloeien uit ver- 
schuivingen die de vraag van de potentiële afnemers en de concur- 
rentieperiode bepalen. Vrijwel altijd zal de afzetpositie van een 
onderneming geleidelijk door deze verschuivingen wijzigen. Ook 
wanneer tussen elk tweetal opeenvolgende perioden slechts geringe 
verschillen bestaan, zal veelal op lange termijn een belangrijke toe- 
of afneming van de verkoopmogelijkheden zijn te constateren. 
Doorgaans zijn de veranderingen in de marktbepalende factoren 
daarom reeds reden genoeg om onderscheid te maken tussen de 
afzetfunktie op lange termijn en die op korte termijn. 

Een tweede oorzaak van verschillen tussen beide afzetfunkties kan 
een verandering van de marktpositie zijn. Deze verandering heeft 
de onderneming tot op zekere hoogte zelf in de hand. Een algemene 
bekendheid en reputatie en een kring van vaste afnemers kan alleen 
in de loop van jaren worden verworven. Het hangt van de verkoop- 
politiek af, of deze worden nagestreefd. Een bedrijf, dat doelbewust 
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naar duurzaam contact met de afnemers streeft, zal op den duur 
een verlaging van de verkoopkosten bereiken, terwijl de afzet- 
funktie ten gunste wijzigt. Hier staat evenwel tegenover, dat men 
in sommige perioden aan de vaste klanten moet leveren, terwijl 
tijdelijk op andere deelmarkten een hogere opbrengst kan worden 
verkregen. Het is daarom ook mogelijk, dat men geen vaste relaties 
aanknoopt, maar in elke periode daar verkoopt, waar dan de afzet- 
mogelijkheden het gunstigst zijn. Zowel de eerste als de tweede 
verkooppolitiek kan gericht zijn op het maximaliseren van de 
winst in de loop der jaren; er bestaat dan alleen verschil tussen 
de wijzen waarop men dit doel tracht te bereiken. 

Uit het voorgaande volgt, dat eigenlijk onderscheid moet worden 
gemaakt tussen: 

1) de afzetfunktie op lange termijn onder de veronderstelling, dat 
in de tijd tot de beschouwde periode regelmatig op de betrokken 
deelmarkt wordt verkocht; 

2) de afzetfunktie op lange termijn onder de veronderstelling, dat 
in de tijd tot de beschouwde periode niet of slechts incidenteel 
op de betrokken deelmarkt wordt verkocht. 

Voor de hierna te maken opmerkingen is het onverschillig, welke 
van beide afzetfunkties wordt beschouwd. De meeste onderne- 
mingen, zeker die van enige omvang, zullen naar regelmatige ver- 
koop streven en dus hun politiek op lange termijn baseren op de 
in 1 genoemde funktie. Ter bepaling van de gedachten kan ook 
in het volgende hiervan worden uitgegaan. 


18. Het is duidelijk, dat het schatten van een afzetfunktie zeer 
moeilijk is. Vrijwel nooit zal een afzetprognose kunnen worden op- 
gesteld in de vorm van het in 13 weergegeven stelsel. In de meeste 
gevallen zal de verkoopprognose bestaan uit de opgave van een 
aantal mogelijk te verkopen hoeveelheden met bijbehorende prijzen 
en variabele acquisitiekosten. Een enkele maal kan statistisch een 
verband tussen enige afzetten, prijzen en acquisitiekosten worden 
bepaald. 

Over de methoden ter schatting van de verkoopmogelijkheden in een 








Het verkoopmodel 121 


bepaalde periode (op korte of lange termijn) wordt hier niet ingegaan. 
De lineaire programmering is alleen van betekenis voor de selectie 
van een optimaal verkoopprogramma uit gegeven mogelijkheden. 
Wel moeten nog enkele opmerkingen worden gemaakt over de 
manier waarop een afzetprognose zodanig kan worden ‘vertaald’, dat 
de lineaire programmeringstechniek hierop kan worden toegepast. 


19. Zoals hierna zal blijken is het gewenst de schatting der ver- 
koopmogelijkheden op een deelmarkt uit te drukken in een aantal 
combinaties van afzetten en bruto af fabriek opbrengsten. De wijze 
waarop deze kunnen worden berekend is afhankelijk van de aard 
der verkoopraming. In dit verband kan een zestal categorieën van 
verkoopramingen worden onderscheiden: 


I. De eenvoudigste raming is die, waarbij voor elk product af- 
zonderlijk één mogelijke afzet, één bijbehorende prijs en één bedrag 
aan variabele acquisitiekosten is gegeven. In dit geval kan men 
uit deze gegevens de bruto af fabriek opbrengst per product direct 
volgens 12 berekenen. De mogelijke afzet van goed 7 op deelmarkt 7 
wordt voorgesteld door V4j; de bijbehorende bruto af fabriek op- 
brengst wordt aangegeven door Oj. Als er 7 producten zijn, kan 
de afzetfunktie op deelmarkt j dus worden voorgesteld door 7 
combinaties 

Vy Oy En ERE: y. 
II. Een omvangrijker afzetprognose is die, waarbij van elk goed 
verschillende alternatieve afzetten en bijbehorende prijzen en 
variabele acquisttiekosten zijn bepaald. Men kan dan verschil- 
lende combinaties van afzet en bruto af fabriek opbrengst bereke- 
nen. Deze worden door accenten van elkaar onderscheiden. Zijn 
bijv. drie alternatieven voor de verkoop van artikel 7 op deelmarkt 
gegeven, dan kan de schatting der verkoopmogelijkheden als volgt 
worden voorgesteld. 

Vy 0; 


: Fa Os, 


] ü 


m mm 
i Vi Oi; 
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Soortgelijke combinaties kunnen ook voor de andere producten 
Dam Dern r worden bepaald. 


III. Het is voorts denkbaar, dat per artikel niet een beperkt aantal 
combinaties, doch een functionele relatie tussen afzet en prijs bij 
enkele niveaus van de acquisitiekosten gegeven zijn. Dit verband 
zou door een grafiek als onderstaande kunnen worden weergegeven. 
Door de curve dalend te tekenen wordt verondersteld, dat een 
grotere afzet kan worden gehaald naarmate de prijs lager is. Door 
de curve, die het dichtst bij de oorsprong ligt, met de laagste 
acquisitiekosten Aj te laten corresponderen neemt men aan, dat 
een vergroting van de acquisitiekosten eveneens de afzetmogelijk- 
heid vergroot (zie echter 15). 


afzet 


Ali < Allj < All, 


prijs 





Uit het door de gegeven funktie beschreven interval kiest men nu 
een aantal afzetten, bijv. drie V;,, V‚; en Vj. Bij elke afzet bepaalt 


men de mogelijke bruto af fabriek opbrengsten voor de verschil- 
lende combinaties van prijzen en variabele acquisitiekosten. Van 
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de bij elk der afzetten behorende bruto af fabriek opbrengsten kiest 
men de maximale. Stel deze maximale af fabriek opbrengsten resp. 
voor door O;;, O;; en O;,. De afzetrelatie kan dan worden beschreven 
door de drie combinaties 


4 , 
Vi Os; 
Nn ” 
Vi Oi; 
HIT Hr 
Vi fi 


Door deze beschrijvingswijze wordt uiteraard niet volledig van de 
beschikbare informatie gebruik gemaakt. Dit is nodig om de 
afzetraming voor de toepassing van lineaire programmering ge- 
schikt te maken. Men zou evenwel de curve willekeurig dicht kunnen 
benaderen door een groter aantal punten te kiezen. 


IV. Bij de in I, II en III genoemde verkoopramingen per product 
wordt de interdependentie tussen de verkoopmogelijkheden van de 
verschillende producten — zoals beschreven in 13 — buiten beschou- 
wing gelaten. Wil men hiermede rekening houden, dan dient men 
niet afzetten en opbrengsten per product afzonderlijk te bepalen, 
doch vaste combinaties van afzetten van verschillende producten 
met bijbehorende prijzen en acquisitiekosten. De verkoopprognose 
moet in deze vorm worden opgesteld, bijv. wanneer het noodzakelijk 
is een volledig assortiment aan te bieden en er een vaste verhouding 
tussen de vraag naar de verschillende artikelen bestaat. In deze 
gevallen zou naar analogie van $ 1, 11 gesproken kunnen worden 
van gemeenschappelijke verkoop. 

De verkoopmogelijkheden op deelmarkt j kunnen in dit geval wor- 
den beschreven door een vector V; van bij elkaar behorende hoeveel- 
heden vij, d= 1 +"……: 7. 


V;= | vy 
Urj 


Bij een dergelijke vector zijn r prijzen P14 +:*-: P‚; en hetzij r 
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variabele acquisitiekosten A1j *:*"- A+j, hetzij een totaal bedrag 
aan acquisitiekosten A; gegeven. Hieruit kan een totale bruto af 
fabriek opbrengst voor deelmarkt j worden berekend. Deze wordt 
voorgesteld door Oj. De afzetfunktie kan dan worden weergegeven 


door de combinatie 
V; 0; 

V. In II is als categorie van verkoopschattingen die genoemd, 
waarbij van elk goed een aantal alternatieve afzetten is bepaald. 
Hiernaast kunnen nu die verkoopschattingen worden onderschei- 
den, waarbij een aantal alternatieve vaste combinaties van af- 
zetten van verschillende goederen is gegeven. De laatste ramingen 
kunnen worden voorgesteld door een aantal vectoren Vj; en bij- 
behorende bruto af fabriek opbrengsten 0. Zijn bijv. drie alter- 
natieve verkoopmogelijkheden voor deelmarkt j gegeven, dan kun- 
nen de volgende combinaties worden bepaald: 


Vj 0j 
V; 9 
V; o 7 


VI. De meest omvangrijke verkoopraming is die, waarbij voor een 
bepaalde deelmarkt een vergelijkingenstelsel is bepaald van de- 
zelfde vorm als dat in 13. Vanzelfsprekend zal dit uiterst zelden 
voorkomen. Zou men voor deelmarkt j over een dergelijk stelsel 
beschikken, dan zou men hieruit voor een aantal willekeurige 
combinaties van afzetten de bijbehorende maximale bruto af fa- 
briek opbrengsten kunnen berekenen. Geeft men de combinaties 
der afzetten weer door een vector Vj; en stelt men de bij een der- 
gelijke vector behorende totale bruto af fabriek opbrengst voor 
door O;, dan verkrijgt men als benadering van de continue funkties 
eenzelfde beschrijving van de verkoopmogelijkheden als in V. 


20. Uit 19 blijkt, dat elke raming der verkoopmogelijkheden op 
deelmarkt j kan worden beschreven, hetzij door één of een aantal 
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afzetten en af fabriek opbrengsten voor elk product afzonderlijk 


Vi Ojj 
Vi Ög pim lenses r 
Vv” o” 


hetzij door één of een aantal vectoren, bestaande uit afzetten van 
verschillende producten als elementen en totale bruto af fabriek 
opbrengsten 


vj 04 
Vv! or 
Vv d 0 zi 


21. Om de verkoopmogelijkheden door een stelsel van lineaire 
vergelijkingen te kunnen beschrijven, moet nu nog één stel begrip- 
pen worden ingevoerd. Het is vereist, de verkoopprognose niet te 
formuleren in termen van afzetten en totale bruto af fabriek op- 
brengsten, maar in termen van additionele afzetten en gemiddelde 
(additionele) bruto af fabriek opbrengsten. Om na te gaan op welke 
wijze dit kan geschieden, worden voorlopig alleen de in 191, II 
en III genoemde categorieën van verkoopramingen beschouwd; 
de andere komen in 25 e.v. ter sprake. 

De gemiddelde bruto af fabriek opbrengst, behorende bij een be- 
paalde afzet, is eenvoudig gelijk aan de totale bruto af fabriek 
opbrengst gedeeld door de afzet. Een gemiddelde additionele op- 
brengst wordt gedefinieerd als het verschil tussen een gegeven 
totale opbrengst en de vorige gegeven totale opbrengst gedeeld 
door de bijbehorende additionele afzet; onder een additionele afzet 
wordt verstaan het verschil tussen een gegeven afzet en de vorige 
gegeven. afzet. 

Stel, de volgende combinaties zijn gegeven 


Vij Oi’ 

Vi O5 pim lee: r 

LA 05 
7 
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Hieruit kunnen de volgende additionele afzetten, vij, v;; en v#, en 
de volgende gemiddelde additionele bruto af fabriek opbrengsten 
oij, oj en 0, worden afgeleid. 





n Os n 
Vij = Vij Fi == Oij 
Á 


Vi == Vij == Vig De Kk = Oij 
m „ 
s-i BR 
Kr] 2 
22. Tenslotte wordt nog een veronderstelling ten aanzien van de 
verkoopramingen ingevoerd. Er wordt verondersteld, dat de ge- 
middelde (additionele) opbrengst tussen O en de eerste gegeven 
afzet vij en tussen elk een tweetal opeenvolgende afzetten vij en 
vi, enz. constant is. Er wordt m.a.w. aangenomen, dat wanneer 
men een hoeveelheid vj tegen een gemiddelde (additionele) op- 
brengst 04j kan verkopen, ook een hoeveelheid kleiner dan v4 
tegen deze opbrengst kan verkopen. De hoeveelheid vij wordt dus 
beschouwd als het maximum, dat tegen de bruto af fabriek op- 
brengst 04 kan worden afgezet. 
Dit betekent, dat het verband tussen de bruto af fabriek opbrengst 
en de afzet van een goed op een bepaalde deelmarkt tussen elk 
tweetal opeenvolgende combinaties 04 en vij, 04; en vij, enz. lineair 
wordt benaderd. Voor twee componenten van de bruto af fabriek 
opbrengst is dit exact juist: de variabele distributiekosten en de 
som van overheidsheffingen en -restituties zijn inderdaad lineair 
van de afzet afhankelijk. Er zal daarentegen in het algemeen geen 
lineair verband bestaan tussen de prijs en de afzet en tussen de 
variabele acquisitiekosten en de afzet. Ten aanzien van deze com- 
ponenten van de bruto af fabriek opbrengst betekent de hier ge- 
maakte veronderstelling een vereenvoudiging van de werkelijkheid ; 
zij komt erop neer, dat een verandering van de variabele acquisitie- 
kosten bij een gegeven prijs een evenredige verandering van de 


afzet tot gevolg heeft. 
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23. Stel, dat er r producten en s deelmarkten zijn. Zou per product 
en per deelmarkt slechts één combinatie vij en oy zijn bepaald 
(ramingen volgens 19 I) dan zouden dus 7 X s van deze combina- 
ties zijn gegeven. Zouden per product en per deelmarkt verschillende 
combinaties vij en 04, v‚j en 0;;, enz. zijn bepaald (ramingen volgens 
19II en 19 II), dan zouden meer dan fr X s combinaties zijn 
gegeven. Om niet steeds accenten bij de symbolen te moeten 
schrijven wordt hieronder alleen het eerste geval beschouwd. Het 
zal hieruit onmiddellijk duidelijk zijn, welke uitbreiding in het 
tweede geval aan de formules kan worden gegeven. 


24. De kolomvector met als elementen de r X s waarden vij wordt 
dan voorgesteld door wv. Dus 


De van artikel # op deelmarkt j verkochte hoeveelheid wordt 
voorgesteld door 44. 

De respectievelijke eenheidsvectoren van dezelfde orde als v worden 
voorgesteld door 


Het totaal van de verkoopmogelijkheden gedurende een periode 
als bedoeld in $ 1, 5 kan dan worden omschreven als de verzameling 
van vectoren # = [411 ***** Xrs] die voldoen aan 
r 8 

u Yrpeid Sv. 

i=l j=l 
Een bepaalde # > 0 wordt een verkoopprogramma genoemd. De 
bruto af fabriek opbrengst van elk der mogelijke verkoopprogram- 
ma's # is gelijk aan 


r 8 


E Dot. 


i=l j=l 
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Met dit stelsel is een lineaire beschrijving van de verkoopmogelijk- 
heden en af fabriek opbrengsten verkregen. Deze beschrijving is 
vergelijkbaar met die van de productiemogelijkheden en productie- 
kosten in $ 1, 20, zij het ook, dat het verkoopmodel een minder ge- 
detailleerd beeld van de werkelijkheid geeft dan het productiemodel. 


25. In de in I9IV, V en VI genoemde gevallen, waarin vaste 
afzetcombinaties V; van verschillende producten zijn gegeven, zijn 
de berekeningen wat ingewikkelder dan de voorgaande. Het model 
verandert echter niet wezenlijk. 
Als voor een deelmarkt slechts één vector Vj is gegeven (zie 19 III), 
behoeft niet het begrip additionele afzet te worden gedefinieerd. 
Overeenkomstig de definities van 21 is V; = vj. 
Om een gemiddelde bruto af fabriek opbrengst te kunnen bepalen, 
moet in dit geval het begrip ‘eenheidspakket’ (zie $ 1, 11) worden 
ingevoerd. Hieronder wordt verstaan een vector 
En V; = En U; 

€ 8 
waarbij e; een getal is. In principe kan het getal e; willekeurig 
worden gekozen. Doorgaans zal men het echter zo kiezen, dat één 


Ì ‚ 
element van — vj (en wel het element, corresponderend met het 
€ 


belangrijkste product) gelijk aan 1 is. De overige elementen geven 
dan de afzetten der overige producten in verhouding tot de afzet 
van dit product aan. 

De gemiddelde bruto af fabriek opbrengst, behorende bij de hier- 


boven gedefinieerde ‘eenheidspakket’ wordt vervolgens gedefi- 
0 
nieerd als + = 0. 
8 
Tenslotte kunnen de verkoopmogelijkheden worden beschreven 
door de combinatie 
1 
Es. ke 
€ 
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Als verschillende vectoren Vj, V‚ ‚ Vs enz. voor één deelmarkt 
zijn gegeven, worden de additionele afzetten gedefinieerd als de 
verschillen hiertussen, dus resp. 


Vj= vj, V; — Vj= vj, Vj — V; = vj, enz. 


Verder worden verschillende ‘eenheidspakketten’ ingevoerd, en wel 
u Vi vj Va Vi vj Vi; — V; 
rt == Dj 2 Om 17 > 

ej 


€; €; ë; 





De getallen ef, e;, e;, enz. zullen in de regel zo worden gekozen, 
dat in elk van de ‘eenheidspakketten’ hetzelfde element gelijk aan 
1 is. 
Tenslotte kan men de gemiddelde additionele bruto af fabriek 
opbrengsten definiëren als 

== 0 gp == 0, 


04 Be Ek 
heer » 3 î mn mek 0, 
Ee €; €; 


! 


es OU 
b] 

De afzetmogelijkheden kunnen tenslotte worden beschreven door 
de paren additionele afzetten en gemiddelde additionele op- 
brengsten 


vj 0j 
n „ 
v 0 5 


26. Evenals in 22 is verondersteld, dat de gemiddelde bruto af 
fabriek opbrengst tussen een tweetal opeenvolgende afzetten 
Vij en Vi, constant is, wordt thans verondersteld, dat deze op- 
brengst voor elke convexe lineaire combinatie van een tweetal 
opeenvolgende combinaties van afzetten Vjen V‚; constant is. 


27. Voeg nu weer alle w; resp. vj enz. aaneen tot een vector v. 
Voer vectoren ui in van dezelfde orde als v en bestaande uit één 


8 1 l : 
van de ‘eenheidspakketten’ — vj, resp. — vj enz. en verder uit 
€ € 
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: 1 l Á 
nullen, zodanig dat de elementen van — © resp. — vj enz. in 
€ Es 


dezelfde rijen staan als de overeenkomstige elementen van ®; resp. 
vj enz. in U. 
Stel tenslotte het aantal malen dat ‘eenheidspakket’ vj (algemeen 


1 1 ' 
voor — V, —v' enz.) wordt uitgevoerd, voor door u}. 

€ 8 
De verkoopmogelijkheden kunnen dan worden beschreven door 
het volgende lineaire vergelijkingen. stelsel 


Ee Zuu sv 
hj 
E Eojuj 
hj 

ui 20. 


$3. Verkoop- en produktieplanning 


1. Onder een verkoop- en productieprogramma wordt verstaan 
een opsomming van de op de verschillende deelmarkten te ver- 
kopen hoeveelheden en de door de verschillende productieafdelin- 
gen te fabriceren hoeveelheden. Een optimaal verkoop- en productie- 
programma is een programma dat onder de gegeven omstandigheden 
de maximale winst oplevert. Verkoop- en productieplanning is het 
vaststellen van een dergelijk programma. 


2. Er kan onderscheid worden gemaakt tussen 

a) de planning op lange termijn, d.w.z. de vaststelling van de na 
een aantal, bijv. vijf jaren te verkopen en te produceren hoeveel- 
heden; 

b) de planning op korte termijn, d.w.z. de vaststelling van de in 
het eerstvolgende jaar te verkopen en te produceren hoeveelheden. 
In feite kan niet bij deze indeling worden stilgestaan. Na de planning 
‘van de in het volgende jaar te verkopen en te produceren hoeveel- 
heden volgt de verdeling van deze hoeveelheden over de maanden 
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rekening houdend met noodzakelijke stilstandstijden van de instal- 
laties wegens revisies, vacanties enz. en met een eventuele seizoen- 
beweging in de afleveringen. Op de maandplanning volgt de detail- 
planning als directe voorbereiding van de uitvoering door de diverse 
afdelingen; deze detailplanning zou bijv. eens per week kunnen 
plaats vinden. 

In het volgende wordt uitsluitend enige aandacht geschonken aan 
de onder a en b genoemde onderdelen van de planning. Hiermede 
wil evenwel niet worden gesuggereerd, dat lineaire programmering 
niet met succes zou kunnen worden toegepast op de maandplanning 
en de hiermee samenhangende onderhoudsplanning en op de detail- 
planning. 


3. Bij de planning op lange termijn wordt beslist over investeringen 
in duurzame activa en over wijzigingen in de verkooporganisatie. 
Hiertoe dient men in de eerste plaats een inzicht te verkrijgen in 
het met het bestaande produktie- en verkoopapparaat te behalen 
toekomstige resultaat. Uit de gegeven produktiemogelijkheden en 
de raming van de afzetfunctie op lange termijn berekent men het 
optimale verkoop- en produktieprogramma door afweging van de 
variabele kosten en bruto af fabriek opbrengsten van de verschil- 
lende alternatieven. Bij deze berekening kunnen de vaste produktie- 
en verkoopkosten buiten beschouwing blijven, daar deze niet door 
de keuze van het programma worden beïnvloed. 

Nadat het optimale programma bij de bestaande installaties en de 
bestaande verkooporganisatie is bepaald, kan worden overwogen 
of het voordelen zou opleveren wijzigingen in de installaties en/of 
in de verkooporganisatie aan te brengen. Daarbij dient in de eerste 
plaats te worden berekend, welke veranderingen in het optimale 
programma uit deze wijzigingen zouden voortvloeien en hoe het 
verschil tussen bruto af fabriek opbrengsten en variabele kosten 
hierdoor zou veranderen. Voorts dient te worden nagegaan, welk 
bedrag moet worden geïnvesteerd om de wijzigingen in de capaciteit 
te bewerken en welke verandering van de vaste jaarkosten (lonen, 
huur, enz.) hiermede gepaard gaat. Tenslotte moet de investering 
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worden afgewogen tegen de verandering in de ‘cash flow’ per jaar. 
Om deze afweging te vergemakkelijken kunnen verschillende ken- 
getallen worden berekend, zoals de terugverdientijd, de ‘discounted 
cash flow rate’, de ‘present worth value’, etc. 


4. De planning op korte termijn bepaalt welke acties door een 
gegeven produktieapparaat en verkooporganisatie zullen worden 
ondernomen. 

Gegevens voor de planning op korte termijn zijn in de eerste plaats 
de capaciteiten der produktiemiddelen en de afzetfunctie op korte 
termijn. Streeft men in elk jaar naar de maximale winst op korte 
termijn zonder rekening te houden met eventuele gevolgen van 
de verkoop in een bepaald jaar voor latere verkoopmogelijkheden, 
dan zijn dit de enige gegevens Voor de planning op korte termijn. 
Het investeringsprogramma en het verkoop- en produktiepro- 
gramma voor de eerstvolgende periode worden dan geheel los 
van elkaar opgesteld. De planning op lange termijn wordt geba- 
seerd op de in $ 2, 17, a2 genoemde afzetfunctie. 

Heeft men daarentegen als verkooppolitiek de opbouw van duur- 
zame marktposities, dan zullen naast de produktiecapaciteiten en 
de verkoopmogelijkheden op korte termijn ook bepaalde uitkomsten 
van de planning op lange termijn gegevens vormen voor de plan- 
ning op korte termijn. Het is bijv. denkbaar, dat men met het oog 
op verkoopmogelijkheden in latere jaren een hoeveelheid in het 
eerstvolgende jaar tegen verlies moet afzetten terwijl elders winst- 
gevende verkoopmogelijkheden bestaan. De uit de verkooppolitiek 
op lange termijn voortvloeiende afzetten in het eerstvolgende jaar 
vormen dan een gegeven voor de planning op korte termijn. De 
planning op lange termijn wordt gebaseerd op de in $2, 17, al 
bedoelde afzetfunctie. 


5. Het uitgangspunt zowel voor de planning op lange termijn als 
voor de planning op korte termijn is de bepaling van een optimaal 
verkoop- en produktieprogramma bij gegeven produktiecapaciteiten en 
verkooporganisatie. In beide gevallen worden dezelfde produktie- 
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mogelijkheden beschouwd. Deze kunnen worden beschreven door 
een stelsel als het in $ 1, 20 genoemde 


EHaD <b 
ì 


5 ki 
1 
A >= 0. 


Bij de planning op lange termijn wordt de afzetfunctie op lange 
termijn beschouwd en bij de planning op korte termijn de afzet- 
functie op korte termijn. Beide afzetfuncties kunnen worden be- 
schreven door eenzelfde stelsel, hetzij dat van $ 2, 24 hetzij dat 
van $ 2, 27. Deze stelsels zijn hieronder nog eens weergegeven: 


EE med <w z 5 ujuj < 
dd 
5 DL 0yky Xx 5 of u} 
jd h j 
xy =O 74 > 0. 


De w in korte termijn modellen zal afwijken van de v in lange ter- 
mijn modellen om de in $ 2, 17 genoemde redenen. Wanneer naar 
duurzame marktposities wordt gestreefd zullen bovendien om de 
in 3 genoemde reden verschillende vergelijkingen in korte termijn 
modellen gelijkheden zijn; in lange termijn modellen zullen daaren- 
tegen meestal alleen < tekens voor de vij voorkomen. 


6. In 1 is de verkoop- en produktieplanning gedefinieerd als het 
vaststellen van de te produceren en de te verkopen hoeveelheden. 
Tussen de som van de te verkopen hoeveelheden van een bepaald 
artikel en de te produceren hoeveelheid kan een verschil bestaan. 
Dit kan worden veroorzaakt door aankoop van het betrokken 
artikel, door voorraadintering en/of door voorraadvorming. Met 
deze drie factoren kan gemakkelijk rekening worden gehouden. 

Als een produkt kan worden ingekocht, voegt men aan b een ele- 
ment toe gelijk aan de te kopen hoeveelheid. Aan elk van de ko- 
lommen a® voegt men een nul toe. Voorts voert men een nieuwe 
kolom a@® in, gelijk aan een eenheidsvector van dezelfde orde als b, 
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waarvan het element 1 in dezelfde rij staat als de te kopen hoeveel- 
heid in b. De variabele kosten A4, behorende bij deze nieuwe a®, 
zijn gelijk aan de inkoopprijs cif plus de variabele inladings-, 
opslag- en afleveringskosten per eenheid produkt. 

Met voorraadintering kan op dezelfde wijze rekening worden ge- 
houden door aan b een element toe te voegen, gelijk aan de maxi- 
maal mogelijke voorraadintering, en als nieuwe kolom a® een een- 
heidsvector van dezelfde orde als b in te voeren. De variabele kosten 
kj zijn in dit geval gelijk aan de geschatte constante waarde van de 
opbrengst in een latere periode minus de variabele voorraadkosten. 
In vele gevallen betekent huidige verkoop, dat later meer moet 
worden geproduceerd. De variabele kosten kunnen dan worden 
gelijkgesteld aan de (toekomstige) produktiekosten minus de varia- 
bele voorraadkosten. Is het produkt bederfelijk, dan kan A; veelal 
gelijk aan O worden gesteld. 

Aankoop en voorraadintering brengen dus een verandering in het 
produktiemodel. Voorraadvorming daarentegen brengt een wijzi- 
ging in het verkoopmodel. Er kunnen hierbij twee gevallen worden 
onderscheiden. 

Speculatieve voorraadvorming met het oog op verkoop in latere 
perioden kan in het model worden geïntroduceerd. door aan v 
een element toe te voegen gelijk aan de te vormen voorraad. Voorts 
voert men een nieuwe eenheidsvector ed), resp. een nieuwe vector 
u} van dezelfde orde als v in. De bruto af fabriek opbrengst 04; 
resp. o} is gelijk aan de contante waarde van de toekomstige op- 
brengst minus de variabele voorraadkosten. 
Naast de speculatieve voorraden kunnen de technische voorraden 
worden onderscheiden, die nodig zijn om onvoorzienbare variaties 
in productie en aflevering op te vangen. Het verkopen op een 
bepaalde deelmarkt impliceert vaak het aanhouden van een voor- 
raad. De afzet en de voorraadvorming moet men dan als een vaste 
combinatie opvatten. In het model kan hiermede op dezelfde wijze 


rekening worden gehouden als met de vaste combinatie van twee- 


afzetten: door introductie of uitbreiding van een vector v} en door 
introductie of uitbreiding van een vector u}. De technisch nood- 
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zakelijke voorraad heeft geen opbrengst; integendeel zal de ge- 
middelde bruto af fabriek opbrengst van u} moeten worden ver- 
minderd met de variabele opslagkosten van de technische voorraad. 
Samenvattend kan worden opgemerkt, dat het model niet wezen- 
lijk verandert door aankoop, voorraadintering en voorraadvorming. 
Er kan dus steeds worden uitgegaan van de in 4 weergegeven 
algemene formulering. 


7. Ter bepaling van het optimale verkoop-en produktieprogramma 
dienen het produktiemodel en het verkoopmodel met elkaar in 
verband te worden gebracht. Om na te gaan hoe dit kan geschie- 
den, wordt eerst het eenvoudige verkoopmodel van $ 2, 24 be- 
schouwd. Voorts wordt aangenomen, dat niet van gemeenschappe- 
lijke produktie (zie $ 1, 11) sprake is. 

In $ 1, 14 is opgemerkt, dat voor elke produktiewijze een produktie- 
activiteit a® wordt ingevoerd. Neemt men even aan, dat elk pro- 
dukt slechts op één manier kan worden gefabriceerd, dan zijn er 
evenveel produktieactiviteiten als produkten. 

Uit de produktieactiviteiten en de eenheidsvectoren ed uit het 
verkoopmodel vormt men nu nieuwe vectoren 


a® 
e(éî) 


Ee Ede r (aantal produkten) 
jm eee s (aantal deelmarkten). 


Het probleem van de vaststelling van het optimale verkoop- en 
produktieprogramma komt dan neer op de bepaling van de waar- 
den xj zodanig, dat 


EE fad) <[b 
td en v 
3E XE (0 — ki)Kij = max 
di 

ij > 0. 


Il 








Dit is een lineair programmeringsvraagstuk, dat met behulp van 
de simplex methode kan worden opgelost. Hiermede is dus een 
algemene oplossingsmethode verkregen voor planningvraagstukken 
van de hierboven omschreven soort. 
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Als per produkt verschillende basisactiviteiten bestaan doordat een- 
zelfde goed op verschillende manieren kan worden gemaakt, moet 
bovenstaand model worden uitgebreid. Dit komt erop neer dat 


a® 
het aantal vectoren en wordt vergroot, omdat aan eenzelfde 


els verschillende a® moeten worden toegevoegd. Het aantal 
variabelen xj neemt in dezelfde mate toe met zij, x;, enz. 


8. Is een verkoopmodel als dat in $ 2, 27 gegeven, dan dient eerst 


| 
te worden bepaald, op welke wijze(n) de hoeveelheden — v? in 
e | 
de ‘eenheidspakketten’ u} van het verkoopmodel kunnen worden 
gefabriceerd. Hiertoe bepaalt men alle lineaire combinaties van 


1 
de produktieactiviteiten a@® die —- wv? opleveren. Stel, bij elk 
er 


1 n 
eenheidspakket van het verkoopmodel behoort slechts één van 
deze combinaties. De respectievelijke combinaties worden voor- 
gesteld door 

l 
X0a) =— vt 

ï € 
Men vormt dan de vectoren 


5 67,a® 
kj ut 


De bepaling van het optimale verkoop- en produktieprogramma is 
dan wiskundig equivalent aan de oplossing van het volgende 


probleem: 
Eu f20)<fb 
hj ï 
ui v 





> 3, (% == 5 Olza) =—= max. 
î : i 


u =O. 
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Ook dit is weer een normaal lineair programmeringsvraagstuk, 
waarvan de oplossing in de voorgaande hoofdstukken is be- 
handeld. 


9. Bij gemeenschappelijke produktie wordt het model betrekkelijk 
ingewikkeld. In dat geval geeft een produktieactiviteit a® volgens 
$ 1, 11 de opeenvolgende bewerkingen weer, vereist om een bepaald 
‘eenheidspakket’ te fabriceren. Men zou dan alle combinaties van 


' 1 3 
de eenheidsvectoren el@ of de vectoren — vl uit het verkoop- 
€ 


] 
model moeten bepalen, die als uitkomst a® (of eventueel een 
lineaire combinatie van a® en een aantal andere produktieactivi- 


f at l 
teiten) opleveren. Elk van deze combinaties van eid of — v?, 
ë: 


zou men met a® (resp. met een lineaire combinatie van a@ he een 
aantal andere produktieactiviteiten) moeten samenvoegen tot 
vectoren van gelijke aard als die in 7 en 8. Met behulp hiervan zou 
het planningvraagstuk weer als een lineair programmeringsprobleem 
kunnen worden geformuleerd. 


10. Alvorens nog enkele opmerkingen te maken over de manier 
waarop het verkoop- en produktiemodel kan worden gebruikt na 
de bepaling van het optimale programma, wordt het voorgaande 
met een voorbeeld toegelicht. Daarbij wordt een kunstmestfabriek 
beschouwd welker produktieproces door het volgende schema kan 
worden voorgesteld (zie pag. 138). 

Daar door elke afdeling slechts één produkt wordt gefabriceerd 
kunnen de jaarcapaciteiten onmiddellijk worden uitgedrukt in een 
aantal eenheden produkt. Dit aantal wordt gevonden door het 
aantal effectieve draaiuren per jaar te vermenigvuldigen met de 
uurcapaciteit. Ook de variabele tijdkosten kunnen direct per een- 
heid produkt worden uitgedrukt door de uurkosten te delen door 
de uurcapaciteit. 

Uitgegaan wordt van de volgende (natuurlijk gefingeerde) gegevens. 
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Ammoniak Zwavelzuur 








Salpeterzuur 


Kalkammonsalpeter Kalksalpeter Zwavelzure ammoniak 
(KAS) (KS) (ZA) 





Opslagloodsen Opzakmachines 
Verladingen 


Jaarcapaciteit Variabele tijdskosten 


per ton 

ammoniakfabriek 75 000 t f 50 
zwavelzuurfabriek 80 000 „, ' - 30 
salpeterzuurfabriek 85 000 „, - 40 
KAS-fabriek 150 000 „ - 30 
KS-fabriek 60 000 „, - 20 
ZA-fabriek 120 000 „, - 30 
Loodsen, opzakmachines 

verladingen — - 20 


(per t, KAS, KS of ZA) 


Er wordt aangenomen, dat de capaciteiten van de opslagloodsen, 
opzakmachines en verladingen zo groot zijn, dat deze nimmer de 
produktie of de aflevering beperken. 

Uit het produktieschema blijkt, welke interne leveringen plaats 
vinden. In het nevenstaande zijn de verbruiken van eigen tussen- 
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produkten per geproduceerde eenheid weergegeven. Voorts zijn de 
overige grondstofkosten per geproduceerde eenheid vermeld. 


Verbruik van eigen Overige 
tussenprodukten grondstofkosten 
per ton per ton 
ammoniakfabriek — f 150 
zwavelzuurfabriek — - 60 
salpeterzuurfabriek 0,3 t ammoniak - 20 
’ 0,2 t ammoniak - 40 
knn es t salpeterzuur 
KS-fabriek 0,6 t salpeterzuur - 60 
ZA-tabriek 0,3 t ammoniak - 80 
0,8 t zwavelzuur 


Het afzetgebied wordt verdeeld in twee deelmarkten: ‘Bouwland’ 
en ‘Grasland’. Per deelmarkt is voor elk produkt één combinatie 
van afzet en totale bruto af fabriek opbrengst bepaald. Er wordt 
hier dus uitgegaan van de meest eenvoudige verkoopprognose 
(zie $ 2, 19 I). 

Uit de gegeven verkoopraming kunnen de volgende combinaties 
van afzetten en gemiddelde bruto af fabriek opbrengsten worden 
berekend. Het blijkt, dat KS niet voor de bemesting van grasland 
wordt verkocht. 







Grasland 


Gemiddelde 
Afzet | bruto af fabriek 
opbrengst per ton) 





Bouwland 


Gemiddelde 
Afzet bruto af fabriek 
opbrengst per ton 
























80 000 t 


60 000 t 







100 000 t 
50 000 t 
80 000 t 


f 340 
f 300 
f 370 
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Er wordt verondersteld, dat bovenstaande afzetten maxima aan- 
geven; bij de genoemde gemiddelde bruto af fabriek opbrengsten 
kan dus naar keuze ook een kleinere hoeveelheid worden verkocht. 
(Voor alle elementen van wv geldt het < teken). 


Uit de gegevens kunnen nu in de eerste plaats de productieactivi- 


teiten a® worden berekend. In dit geval zijn er drie: 


KAS KS ZA 
Ammoniak (0,38) (0,2+0,6x0,3) [0,18] (0,6x0,3) [0,3 
Zwavelzuur 0 0 0,8 
Salpeterzuur | 0,6 0,6 0 
KAS 1 0 
KS (0) | 6) 
ZA (0) °) 1 


De variabele kosten per basisactiviteit kunnen als volgt worden 


berekend : 
KAS KS ZA 

Variabele 0,38xf 50=f 19 0,18xf 50=f 9 O3Xxf SO=f 15 
tijdkosten 0,6X- 40=- 24 06X- 40=- 24 0,8Bx- 30=- 24 

1x- 30=- 30 Ix- 20=- 20 1X- 30=- 30 

1xX- 20=- 20 IxX- 20=- 20 1X- 20=- 20 
Variabele 0,38X-150=- 57 0,18X-150=- 27 0,3X-150—=- 45 
grondstof- 0,6X- 20=- 12 0,6X- 20=- 12 0,8x- 60=- 48 
kosten 1x- 40=- 40 Ix- 60=- 60 1X- W=- 80 
Totaal k; {202 {172 1262 


Er kunnen de volgende verschillen tussen de gemiddelde bruto af 
fabriek opbrengsten 04; en de variabele fabricagekosten k; worden 
bepaald. 
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KAS (1) KS (2) 
Bouwland (1) Grasland (2) Bouwland (1) 
Oi f 340 f 360 f 300 
ki - 202 - 202 - 172 
opk 1138 {158 1128 
ZA (3) 
Bouwland (1) Grasland (2) 
oi f 370 f 350 
ki - 262 - 262 
Oij — Ri 108 f 88 


Wordt de afzet van produkt 7 op deelmarkt j voorgesteld door #44, 
dan komt de bepaling van het optimale verkoop- en produktie- 
programma neer op de oplossing van het volgende lineaire pro- 


grammeringsvraagstuk. 


x1110,38| + x12(0,38| +4 #21[0,18| + #31/0,3| + #32)0,3| < | 75000 


0 0 0 0,8 
0,6 0,6 0,6 0 
Ì Î 0 0) 
0) 0 l 0 
0 0 0 1 
l 0 0) 0) 
0 1 0 0 
0 0 l 0 
0 0 0 | 
0 0 0) 0) 


138x11 + 158x12 + 128x21 + 108431 + 88432 —= max. 
X11, %12, %21, %31, %32 > 0. 


0,8 


0 


Ooo 


80 000 
85 000 
150 000 
60 000 
120 000 
100 000 
80 000 
50 000 
80 000 
60 000 


Na elf verschilvariabelen te hebben ingevoerd, kan men dit pro- 
bleem met behulp van de simplexmethode oplossen. Na een aantal 
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herhalingen verkrijgt men de volgende tabel. (De kolommen van 
de herleide vorm zijn evenals in het vorige hoofdstuk aangegeven 


door a al®) en b.) 


138 158 128 108 88 O O 0 0 0 


ab a) a) a al5) al) al) 


5 —1,9 —1,5 0 
0 0 0 0 
—-5 19 3,17 0 
0 O0 0 0 
0125 O0 0 
—-519 15 1 
5 —1,9 —3,170 
0—1,25 O0 O0 
—-519 15 0 
5 —1,9 —3,170 
01,25 O0 O0 


oooooooooo-=00O 
oooooooo=-O000 
ooooocoO=oooo 
oooo=-oooooo 


1 

0 
5 
k0) 
0 
0 
0) 
0 
0) 
0 
0) 


oooooooooooo-=oO 


138 158 128 108 88 50 90 200 O 
0 O0 O0 O0 0 —50—90—2000 O 


0 


ooo=ooooooo 


0 0 0 0 © 


a® al al1® all) al12) gl13) gl1) gl15) q(16) 


_—_ 


oo=oocoocoooo 
oo=ocooooo-l 

o=ooooocoooo 
=oooool=ocoo 
=ooocoocooooo 


Oo 20 0 
—20 O0 


0 20 
—20 0 





Uit deze tabel blijkt in de eerste plaats, dat het optimale verkoop- 


programma luidt: 


Bouwland Grasland 
KAS 18 000 t 80 000 t 
KS 44 000 t — 
ZA 80 000 t 20 000 t 


De verkooppolitiek moet dus m.a.w. zijn gericht op maximale afzet 

„van KAS in het grasland en ZA in het bouwland. Er moet verder 
naar worden gestreefd, de KS-afzet eveneens vrijwel tot het maxi- 
mum op te voeren. De afzet van KAS in het bouwland en ZA in 
het grasland dient minder te worden gestimuleerd; het is evenwel 
niet gewenst één van beide geheel te verwaarlozen. 
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Tenslotte kan uit de tabel het volgende produktieprogramma 
worden afgeleid: 


ammoniak 75 000 t 
zwavelzuur 80 000 t 
salpeterzuur 85 000 t 
KAS 98 000 t 
KS 44 000 t 
ZA 100 000 t 


De ammoniak-, de zwavelzuur- en de salpeterzuurfabriek zijn dus 
volbezet. Er bestaat een ondercapaciteit in de KAS-, de KS- en 
de ZA-fabriek. 

Het verschil tussen de totale bruto af fabriek opbrengst en de 
variabele fabricagekosten kan uit de maximeringsfunctie worden 
berekend: 


18 000 x f 138 = f 2 480 000 
80 000 Xx - 158 = - 12 640 000 
44 000 x - 128 = - 5 630 000 
80 000 Xx - 108 = - 8 640 000 
20000 x - 88 =- 1 760 000 


f 31 150 000 


Door van dit bedrag de vaste ‘out of pocket expenses’ (lonen, 
huur, enz.) en de afschrijvingen af te trekken vindt men de winst- 
begroting. 


11. Nadat het optimale verkoop- en produktieprogramma en de 
bijbehorende winst zijn bepaald, wordt vaak de vraag gesteld welke 
wijzigingen hierin zouden optreden als de uitgangspunten zouden 
worden veranderd. Deze vraag rijst in de eerste plaats bij de plan- 
ning op lange termijn, wanneer men de rentabiliteit van een ver- 
andering in de capaciteit der produktieafdelingen wil bepalen. Er 
wordt dan een verandering in een element 5; van b overwogen. 
In de tweede plaats kan zowel bij de planning op lange termijn 
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als bij die op korte termijn een verandering in de variabele kosten 
bj moeten worden beschouwd, bijv. wegens wijziging van de grond- 
stoffenprijzen. Voorts kan men in beide gevallen een gelijktijdige 
verandering van een gemiddelde bruto af fabriek opbrengst 044 
en een afzet vij van v beschouwen om te bepalen of stimulering 
van de verkoop in een bepaalde richting voordelen zou opleveren. 
Het bepalen van de invloed van een dergelijke gecombineerde ver- 
andering na berekening van het optimale programma is het alter- 
natief voor het invoeren van een nieuwe verkoopactiviteit in het 
model. Tenslotte is het denkbaar, dat veranderingen in de elemen- 
ten a van de produktieactiviteiten a@® optreden. Dergelijke ver- 
anderingen zouden kunnen worden veroorzaakt door verbetering 
van de efficiency der produktieafdelingen. 


12. De invloed van veranderingen in de gegevens waarop het 
planningmodel is gebaseerd, kan uiteraard altijd worden bepaald 
door een nieuw vergelijkingenstelsel op te stellen en dit opnieuw 
met behulp van de simplexmethode op te lossen. Soms is het 
echter ook mogelijk, het gewijzigde optimale verkoop- en pro- 
duktieprogramma en de daarbij behorende winst af te leiden uit 
de reeds bepaalde oplossing van het oorspronkelijke vraagstuk. 
Dit is o.m. het geval bij een verandering in de capaciteiten by, in 
de variabele kosten 4, in de afzetten vij en in gemiddelde bruto 
af fabriek opbrengsten 04. 
Een wijziging in een capaciteit b4 of in een af te zetten hoeveelheid 
vj van het planningmodel is een voorbeeld van een wijziging in 
een element van b in het algemene model 

Ax < of >b 

“=0 
€'$ = optimaal. 


Een wijziging in de variabele kosten k; of in een gemiddelde bruto 
af fabriek opbrengst 04j betekent daarentegen een verandering in 
een element van c’. 

Er wordt nu nagegaan, welke invloed beide veranderingen op de 
optimale oplossing hebben. 
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‚13. Stel, het oorspronkelijke vraagstuk is opgelost door toepassing 
van de simplexmethode. In dat geval is een optimale oplossing «* 
gevonden als laatste kolom van een tabel waarin alle 
(cj — zj) <0 (resp. > 0). 

Voor het gemak wordt aangenomen, dat de eerste mm elementen 
van x* ongelijk aan 0 zijn. De basisvectoren zijn dan ad) ...-. am, 
Een verandering van een element van b, bijv. een verandering 
van b1 tot (b1 + 461), heeft geen gevolgen voor de verschillen 
(c; — zj); de elementen van b komen immers niet in deze ver- 
schillen voor. Wel zal de optimale oplossing veranderen. Stelt 
men de matrix van de orde m Xx m, bestaande uit ad «--.- am), 
voor door A*, dan voldoet «* aan 


A*x* —= b 
«* — A*-1b 
Dus: 





(Volgens II $ 2, 9 is [Aj] de cofactor van het element a; van |A*].) 
Stelt men de oplossing van het nieuwe vraagstuk, behorende bij 
dezelfde vectoren ad -...- am) als #*, voor door «® dan moet 


gelden 
A*x0—=b + Ei 
0 


el 








Dus: 
TET | en O1 + Abi}. 
Wanneer nu x{ > 0, dan kan #® uit #* worden afgeleid, Immers: 
== + Di An 4b1 
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en de overige (m — 1) elementen van x@ zijn gelijk aan de over- 
eenkomstige elementen van x*. 

Het optimum van de maximerings-(resp. minimerings)functie van 
het tweede vraagstuk is gelijk aan 


hd 
10 Paak ee |A] 
c's0 = C'x* + Ci X 


Abi. 
ii |A*] 





Wanneer x1 < 0, kan de optimale oplossing van het nieuwe vraag- 
stuk niet zo gemakkelijk worden bepaald. In het algemeen zit er 
in dat geval niets anders op dan het nieuwe vraagstuk volgens de 
simplexmethode op te lossen. 


14. Bij overweging van investeringen is men primair geïnteres- 
seerd in veranderingen in de winst door verandering van b en pas 
op de tweede plaats in veranderingen in het optimale verkoop- en 
produktieprogramma. Ter bepaling van alleen de verandering in 
de optimale winst bestaat een eenvoudiger methode dan de in 13 
genoemde. Omdat de afleiding hiervan op een nog niet genoemde 
wiskundige eigenschap is gebaseerd, wordt deze methode in een 
afzonderlijke paragraaf 5 besproken. 


15. In tegenstelling tot een verandering in b heeft een verandering 
in c’ wel invloed op de verschillen (c; — 24). Stel bijv. dat het eerste 
element van c’ verandert van ci tot (c: + Aci). In dat geval ver- 
andert het eerste verschil cj — z1 tot 


(ci + Aci) — (5 mirct + ande) = 
el 


== (c1 — zi) + (1 — mii) Âc1. 


De overige verschillen (c; — zj), j= 1 +-"" n, worden: 


m 
y— (2 Hij: + mgdc.) = (& — 44) — maden. 
i=l 


Wanneer de nieuwe verschillen alle <0 (resp. > 0) zijn, is de 











Transportplanning 147 


oorspronkelijke oplossing x#* ook een optimale oplossing van het 
tweede vraagstuk. Alleen het optimum van de maximerings- 
(resp. minimerings)functie verandert en wel in 


Pak Ed 
c's* + Aex. 


Wanneer de nieuwe verschillen niet alle <0 (resp. >= 0) zijn, is 
«* wel een bruikbare maar niet een optimale oplossing van het 
nieuwe vraagstuk. Uitgaande van de gegeven tabel kunnen de 
berekeningen dan worden voortgezet totdat wel een optimale 
oplossing is gevonden. 


16. Hoe de gevolgen van gelijktijdige veranderingen in b en in c' 
kunnen worden bepaald, zal onmiddellijk duidelijk zijn: Hiertoe 
dient men de berekeningen volgens 13 en 15 achtereenvolgens uit 
te voeren. 


$ 4. Transportplanning 


1. In 1, 81,9 is een voorbeeld gegeven van een transportvraag- 
stuk. In het algemeen staat men bij de Eg voor het 
volgende probleem: 

Om m plaatsen zijn de hoeveelheden aj ---«- am van een zeker 
artikel beschikbaar. Deze hoeveelheden moeten worden vervoerd 
naar # andere plaatsen. Op de laatste plaatsen zijn resp. de hoeveel- 
heden bj ++-- bn nodig. Het totaal van de beschikbare hoeveel- 
heden is gelijk aan het totaal van de benodigde hoeveelheden. De 
transportkosten van plaats # naar plaats j bedragen 


kj (£ = 1 vee ee m;j=l esse n). 


Gevraagd wordt het transportprogramma zodanig te bepalen dat 
de totale transportkosten minimaal zijn. 
Wiskundig kan dit vraagstuk worden geformuleerd als de bepaling 
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van de te transporteren hoeveelheden 411 --"-: Xmn zodat: 
Kir deer Kin = 41 
Kar tret Kan == 42 
Xm1 + “+ Amn = Um 
X11 + X21 + eee eee eeeseesee + Xm1 == bi 
Xin + Xan + vee eee eee eee eee + Xmn —= On 
4D Xmn =O 
Ruz + *…*: + RmnAmn = min, 
waarbij 
a1 + eee + Am = b1 + ee + bm 


2. Aan de transportplanning zijn weinig economische problemen 
verbonden. Laat men de goederen door een derde transporteren, 
dan zijn de kosten kij eenvoudig gelijk aan het door de vervoerder 
berekende vrachttarief plus de kosten van verzekering en even- 
tueel andere bijkomstige kosten. Als men zelf transportmiddelen 
in eigendom heeft en geen andere aanwending hiervoor bestaat 
dan het vervoer van de beschouwde goederen, omvatten de kosten 
kij de variabele tijdkosten en grondstofkosten (olie, benzine, e.d). 
In het algemeen is het echter mogelijk, transportmiddelen te ver- 
huren. De kosten Z4j zullen dan moeten worden gelijkgesteld aan 
de alternatieve opbrengst bij verhuur. 


3. Zoals elk lineair programmeringsvraagstuk kan ook een trans- 
portprobleem door toepassing van de simplexmethode worden op- 
gelost. Wegens de specifieke vorm van de coëfficiëntenmatrix A 
bestaat er echter voor de oplossing van transportproblemen een 
veel eenvoudiger methode, welke hieronder in het kort wordt 
besproken. 


4. De beschikbare en benodigde hoeveelheden en de transport- 
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kosten uit het in 1 genoemde algemene transportprobleem, kunnen 
in een tabel als de onderstaande worden weergegeven. 





Voorbeeld: Stel, de volgende te transporteren hoeveelheden en 
transportkosten zijn gegeven (de hoofdletters stellen plaatsen voor). 


Transportschema X 1000 ton 





afteleveren | te ontvangen 


A 90 T 20 
B 50 18 40 
C 40 Vv 10 
D a0 W 70 
E 10 Xx 20 

Y 30 

Z 50 
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Deze gegevens zouden in de volgende tabel kunnen worden weer- 
gegeven. 





5. Stelling: Hoogstens m + n — 1 elementen van een optimale 
oplossing zijn ongelijk aan 0. 


Bewijs: Het in 1 weergegeven stelsel omvat m + „ vergelijkingen 
met mn onbekenden. De vergelijkingen zijn echter niet onafhanke- 
lijk omdat X4 a; = X4 bj. Telt men eerst de eerste m vergelijkingen 
en vervolgens de volgende n — 1 bij elkaar op en trekt men dan 
de tweede som van de eerste af, dan vindt men: 


Kin + Aan +11" + Amn = 
=a Hager Ham — bi — ba err — br 


Daar X4a; = Xyby is de term rechts van het gelijkteken gelijk 
aan bn. Op deze wijze kan dus de laatste vergelijking uit de eerste 
m + n — Ì worden bepaald. 

Van de coëfficiëntenmatrix A zijn m.a.w. slechts m + n — 1 rijen 
en kolommen orìafhankelijk. Omdat de elementen van een optimale 
oplossing volgens III, $ 2, 10 de gewichten moeten zijn van lineair 
_ onafhankelijke kolommen van A, kan ook een optimale oplossing 

slechts m + n — 1 elementen ongelijk aan O bevatten. 


6: Het eerste probleem is nu, in de hokjes van een tabel als in 3 
is weergegeven m + n — Ì waarden van de variabelen x4 in te 
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vullen, zodanig dat deze waarden de gewichten van lineair onaf- 
hankelijke kolommen van A zijn en dat aan de transporteisen 
is voldaan. Hierbij kan men als volgt tewerk gaan: 

In een willekeurig hokje, bijv. dat in de A-de rij en de s-de kolom, 
plaatst men een getal gelijk aan an als an < bs en gelijk aan bs 
als an > bs. (Het geval waarin an — bs wordt hierna genoemd.) 
Als an < bs en Xns — an, kan in alle andere hokjes van de h-de 
rij een O worden geplaatst. Als an > bs en xs — bs, kan in alle 
andere hokjes van de s-de kolom een nul worden geplaatst. Om 
een hierna te noemen reden vult men de nullen echter in feite niet 
in; de hokjes worden eenvoudig open gelaten. 

Door het eerste hokje in te vullen kan een rij of een kolom verder 
buiten beschouwing worden gelaten. Met de overgebleven tabel 
gaat men vervolgens op dezelfde wijze tewerk: Men vult weer een 
willekeurig hokje in met een waarde an als an < bs en met een 
waarde bs als a, > bs. Er valt dan opnieuw een rij, resp. een kolom 
weg. 

Dezelfde methode past men bij herhaling toe. De laatst overge- 
bleven ‘tabel’ bestaat uit één hokje. Door dit in te vullen vallen 
tegelijkertijd een ‘rij’ en een ‘kolom’ weg. Men heeft dan in totaal 
n + m — Ì variabelen ongelijk aan O bepaald. Hiermede a ook 
de daarbij behorende basisvectoren gevonden. 

Wanneer in enige etappe an — bs, stelt men uiteraard de betrokken 
variabele xps gelijk aan deze waarde. Men heeft dan in het alge- 
meen de keuze om een rij of een kolom te laten wegvallen. Wanneer 
evenwel slechts één rij en verschillende kolommen zijn overge- 
bleven, laat men een kolom wegvallen; wanneer slechts één kolom 
en verschillende rijen zijn overgebleven, laat men een rij weg- 
vallen. In elk geval houdt men na het laten vervallen van een rij 
of een kolom een tabel over, waarvan één kolom resp. één rij als 
totaaltelling O heeft. In een volgende etappe zal in deze kolom 
of rij een O moeten worden ingevuld. In tegenstelling tot de nul- 
waarden der variabelen die niet bij basisvectoren behoren, wordt 
deze O wel ingevuld. Op deze wijze geeft men aan, dat de bij het 
betrokken ‘hokje behorende kolom van A een basisvector is. 
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Voorbeeld: Hieronder is de in 4 weergegeven tabel op boven- 
genoemde wijze ingevuld. Zoals is opgemerkt is de keuze van 
het in elke etappe in te vullen hokje in principe willekeurig. Men 
beperkt de latere berekeningen evenwel zo veel mogelijk door 
steeds het hokje of één van de hokjes met de laagste kosten- 
component het eerst in te vullen. Past men deze methode toe, 
dan kunnen achtereenvolgens in de hokjes ij de volgende waarden 
worden ingevuld: 


. 


== NRN 
AUDNONU RR ANNN == 
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7. Door volgens 6 te werk te gaan vult men de hokjes van de tabel 
op een speciale manier in: Alle rijen en kolommen van de tabel 
bevatten minstens één ingevuld hokje en er is minstens één 
rij of kolom die niet meer dan één ingevuld hokje bevat. Laat 
men een rij of een kolom met slechts één ingevuld hokje weg, 
dan blijft een tabel over waarvoor hetzelfde geldt. Laat men ook 
hiervan een rij of een kolom met niet meer dan één ingevuld 
hokje weg, dan houdt men weer een tabel over waarvoor 
hetzelfde geldt, enz. Uiteindelijk houdt men één ingevuld hokje 
over. 


8. Stelling: Een stel van (m +» — 1) kolommen van het in 
1 weergegeven vergelijkingenstelsel vormen een basis, wanneer 
en alleen wanneer deze corresponderen met hokjes, die op de 
in 7 beschreven wijze zijn verdeeld over de in 3 weergegeven 
tabel. 


Bewijs: 

1. Omdat in elke rij en elke kolom minstens één hokje is ingevuld, 
kan aan elke verandering in a of b worden voldaan door aan de 
in deze hokjes geplaatste variabelen andere waarden toe te kennen. 


Dit betekent, dat elke kolom in de met de gekozen hokjes 


a 
b 
corresponderende vectoren kan worden uitgedrukt. Deze m J- n — 1 
vectoren vormen dus een basis. 


IL. Uit I volgt, dat alleen dan een verzameling basisvectoren kan 
worden verkregen als in elke rij en elke kolom minstens één hokje 
is ingevuld. Voorts kan onmogelijk een stel basisvectoren worden 
verkregen als in alle rijen en kolommen twee of meer hokjes zouden 
worden ingevuld. In dat geval zou het aantal ingevulde hokjes en 
dus het aantal vectoren >> 2m en >= 2n dus >= m + n zijn. Een 
basis bestaat echter slechts uit m +» — 1 vectoren. 

Ook de rijen en kolommen van de tabel, die uit de oorspronkelijke 
ontstaat door daarvan een rij of een kolom met slechts één ingevuld 
hokje weg te laten, kunnen niet alle twee of meer ingevulde hokjes 
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bevatten. Stel bijv, dat een rij met één ingevuld hokje uit de 
oorspronkelijke tabel wordt weggelaten. Zouden dan alle rijen of 
kolommen van de gereduceerde tabel twee of meer ingevulde hokjes 
bevatten, dan zou het aantal vectoren hierin > 2(m — 1) en > 2n, 
dus >= m — 1 + # zijn. Telt men hierbij de vector uit de wegge- 
laten rij van de oorspronkelijke tabel, dan vindt men dat het 
totaal aantal vectoren > m + #1 is. 

Een zelfde redenering kan worden gevolgd voor de tabellen die 
men verkrijgt door verschillende rijen en kolommen van de oor- 
spronkelijke tabel weg te laten. 


9. Uit 8 volgt, dat een volgens 6 bepaalde oplossing een bruikbare 
basisoplossing is. 


10. Stel de kostencoëfficiënten in de ingevulde hokjes voor door 
kj}. Volgens 7 komt in elke rij en in elke kolom van de tabel minstens 


één ingevuld hokje voor. Elke waarde 1 -:--: m wordt dus door 
minstens één van de indices £ van de Ke aangenomen; evenzo wordt 
elke waarde 1 ----- n door minstens één van de indices j van de 
kj aangenomen. 
Definieer nu de grootheden Aj ----- Am, H1****" Hm als de op- 
lossing van het vergelijkingenstelsel 

kj + A + uy =0 izle: m 

jl ” 


Er zijn m + n — 1 kostencoëfficiënten ki en dus ook m + — Ì 
dergelijke vergelijkingen. Er zijn daarentegen m + n variabelen 
Às, Hij. 

_ Volgens II, $ 4, 11 is dit homogene stelsel dus oplosbaar als men 
hetzij aan één A; hetzij aan één mj een bepaalde waarde toekent. 
Zoals in 12 blijkt, moet in dit geval een willekeurige A; of uj gelijk 
aan O worden gesteld. Voorlopig wordt dit zonder bewijs aange- 
nomen. 
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Voorbeeld: Bij de in 6 weergegeven tabel behoort het volgende 
vergelijkingenstelsel: 


44 A1 + ua=0 
S + Â1 + us =0 
9 + A1 + ue =O 
1 + A2 + ui =0 
6 + 2 + us =O 
24 A2 +ur=0 
2 4 Â3 + ur =0 
2 + a + ua =0 
6 + Ja + ug =O 
8 + A4 + u6 =O 
3 + A5 + us =O 


De oplossing van dit stelsel (bij een willekeurige A; of u; gelijk aan 0) 
kan met behulp van de tabel worden bepaald. Daarom pleegt men 
dit stelsel bij feitelijke berekeningen niet expliciet op te schrijven. 
Hieronder is de tabel uit 6 nog eens weergegeven. Men stelt nu een 
willekeurige A; of A; gelijk aan O. Om de berekeningen zo beknopt 
mogelijk te houden kiest men hiervoor die 4 of uy welke bij de 
rij of kolom met de meeste ingevulde hokjes behoort. In dit geval 
is dat As. De waarde Â2 — O schrijft men achter de tweede rij. 
Onmiddellijk blijkt nu uit de in bovenstaande vegelijken weer- 
gegeven definities, dat: 


mi —l 
ue= —6 
ua= —2. 


De waarden schrijft men resp. onder de kolommen 1, 6 en 7. 


Uit ue = —6 leidt men af dat Aj = —3 
A= —2 
en uit uz = —2 dat Âg3 — 0. 


Deze waarden schrijft men resp. achter de eerste, de vierde en de 
derde rij. 








156 De toepassing in het bedrijf 


Uit A1 = —3 volgt ua=—l 
ps —?2 

Uit A4 = —2 volgt ua = 
ua= —4 

Tenslotte volgt uit u5= —2 As = —l. 





11. Beschouw nu weer het vergelijkingenstelsel in de in 1 weer- 
gegeven vorm. Schrijf de a; en bj voor de gelijktekens; elk van deze 
constanten krijgt dan een min-teken. Na het gelijkteken staat in 
elke vergelijking een nul. 


Vermenigvuldig de eerste m vergelijkingen resp. met A1 *"" Am 
en de laatste # vergelijkingen resp. met gi: tn. Men ver- 
krijgt dan 


n 
nl Ems — as) =0 iz laren pn 
jl 


j= 


* . 
(Eu — 6) =0 j=l herdr zi n. 
iet 


Telt men deze produkten, die alle gelijk aan O zijn, op bij de mini- 
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meringsfunctie, dan blijft deze onveranderd. Dus: 


mm Nn m 


YE Ekym=D Lkyjxg 


i=l j=1 i=l j=l 


+ Ae ĳ— a) +3 ml > zij — o) 
i=1 = j= 


m n me n 
=DE D (kij + Ae + 9) — Zhai — Di sbs. 
=1 j=1 i=1 j=l 
Omdat de factoren A; en uj zijn bepaald uit de vergelijkingen 
hi + Me + us =0 


volgens 10, zijn de termen van 


r: (lij + À + Hy) wij 


„Ms 


1 


voor de basisvariabelen gelijk aan O. In deze som zijn dus m.a.w. 
de termen, behorende bij de basisvectoren, weggevallen. De ter- 
men Xi; Mas en Y-_, MJD; vormen tezamen een constante. Stelt 
men deze voor door C, dan kan men schrijven: 


EE hijxij = E Rp) (ig + Ja + 149) vig + C. 

iel j=l 
Hiermede is de nn geschreven uitsluitend als 
functie van de niet-basisvariabelen. 


12. Zoals in 5 is opgemerkt, zijn slechts m + n — 1 vergelijkingen 
van het in 1 weergegeven stelsel lineair onafhankelijk. Bij de in 11 
uitgevoerde optelling kan derhalve één van de produkten 
MEL. Xij — 44) of één van de produkten u;(X7*, %ij — bj) worden 
weggelaten. Dit betekent, dat een willekeurige A; of uy gelijk aan OQ 
kan worden gesteld. Hiervan is in 10 reeds uitgegaan. 


13. Uit de in 11 bepaalde vorm van de minimeringsfunctie volgt: 
— Wanneer een term (kij + As + uy) < O is, betekent een vergro- 
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ting van de bijbehorende variabele xj, een nintanie van de 
funktiewaarde. 

— Wanneer alle termen (Rij + As + wy) =O zijn kan de functie- 
waarde niet verder worden verlaagd door verandering van de 
variabelen. De gevonden basisoplossing is dus m.a.w. een opti- 
male oplossing. 

Hiermede is dus een methode gevonden om te bepalen of het voor- 

deel heeft een nieuwe vector in de basis te brengen. Tevens is het 

criterium gevonden waaraan een optimale oplossing moet voldoen. 

De eerstvolgende stap in de calculaties na de berekening van de 

grootheden À; en us volgens 10 zal dus zijn het bepalen van de som- 

men (kij + As + uy). Constateert men dat alle (kij + À; + uy) > O 

dan kan het vraagstuk als opgelost worden beschouwd. Is hieraan 

niet voldaan, dan moeten verdere stappen worden gedaan. 

Voorbeeld: In de in 10 weergegeven tabel zijn hieronder de waar- 

den Zy + A; + uy ingevuld linksonder in de verschillende hokjes. 

Dit behoeft alleen te worden gedaan voor de hokjes waarin geen 

basisvariabele voorkomt. Per definitie is immers iedere 


hij + MH u =0. 








4 9 

70 10 —_3 

OERERENE 
50 ò 

3 

SEEKER 
eee 
EEE R 





ENE EEEN 
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14. Stel, dat niet alle (kij + A: + uy) =O zijn. In dat geval be- 
paalt men min(ky + Â: + uy). Wanneer verschillende van deze 
sommen minimaal en gelijk zijn, kiest men een willekeurige. 

Stel, min (kij + As + uy) — (hens + Àn + Hs). Als men dan aan xps 
een positieve waarde toekent, zal de functiewaarde dalen. Er 
wordt daarom besloten a) in de basis te brengen. 

In het hokje As plaatst men nu een letter 0 als symbool voor de 
te bepalen waarde van xps. Zou men aan deze Ò zonder meer een 
positieve waarde toekennen, dan zou de tabel niet meer ‘in even- 
wicht’ zijn, dus de som van de getallen in de h-de rij zou niet 
meer gelijk zijn aan an en de som van de getallen in de s-de kolom 
zou niet meer gelijk zijn aan bs. Om de tabel in evenwicht te bren- 
gen volgt men de volgende procedure: 

Na een 9 in het bokje As te hebben ingevuld, trekt men van de 
waarde in een reeds ingevuld hokje van bijv. de h-de rij af. (Men 
kan ook met de s-de kolom beginnen). De h-de rij is dan in evenwicht. 
De kolom waarin het hokje met — 0 staat, is evenwel uit evenwicht 
gebracht. Om deze onevenwichtigheid op te heffen telt men bij een 
ander, reeds ingevuld hokje van deze kolom 6 op. Is de betrokken 
kolom toevallig de s-de, dan is hiermede de gehele tabel in even- 
wicht. Is deze kolom een andere dan de s-de, dan moet opnieuw 
in een reeds ingevuld hokje van een nieuwe rij +0 worden ge- 
plaatst, waarna in een volgende kolom — 9 moet worden ingevuld. 
Als de laatstbedoelde kolom de s-de is, is het evenwicht na deze 
stap bereikt. Anders gaat men door. Telkens worden in elke rij of 
kolom die niet in evenwicht is, één en niet meer dan één +6 en 
één en niet meer dan één — 0 geplaatst in reeds ingevulde hokjes. 
Hiermede houdt men op zodra de tabel als geheel in evenwicht is. 
Nadat de tabel in evenwicht is gebracht, kan de waarde van 6 
worden bepaald. Hiervoor gelden de volgende richtlijnen: 

a) De daling van de functiewaarde is volgens 11 evenredig aan de 
grootte van xps. Derhalve moet 0 zo groot mogelijk worden gekozen. 
b) De variabelen moeten niet-negatief zijn. Derhalve kan 6 niet groter 
zijn dan de kleinste waarde in de hokjes waarin — 6 is geplaatst. 
c) Bij introductie van a» in de basis moet een andere kolom van 
het in 1 weergegeven stelsel eruit worden verdreven. De variabele 
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in één van de ingevulde hokjes moet dus gelijk aan O worden. 
De conclusie hieruit is, dat 9 moet worden gelijkgesteld aan de 
kleinste waarde in de hokjes waarin — 0 is geplaatst. Substitueert 
men deze waarde voor 0, dan wordt een nieuwe tabel gevonden 
met hoogstens m + n — 1 variabelen ongelijk aan O. De nieuwe 
basisvariabelen leveren een lagere functiewaarde op dan de eerste. 
Indien twee of meer waarden in de hokjes met — 6 minimaal en 
gelijk zijn, worden twee of meer basisvariabelen gelijk aan 0. 
Niettemin laat men in dat geval slechts één variabele en bijbeho- 
rende vector vervallen; in de andere hokjes plaatst men in de 
nieuwe tabel nullen. 


Voorbeeld: Als eerste is hieronder de in 13 weergegeven tabel nog 
eens overgenomen. Nu zijn evenwel de grootheden +60 en — 6 
hierin geplaatst. 

Min (kij + As + uy) is gelijk aan —2 in het hokje h‚s = 1,7. In 
dit hokje plaatst men dus de eerste +0. De posities van de overige 
+ 0's en — 0’s kunnen verder door trial and error worden bepaald. 
Dit is in dit geval zeer eenvoudig. 

De maximumwaarde van 6 blijkt gelijk aan 10 te zijn. (De waarden 
van de variabelen in de hokjes 7,7 = 1, 6 en 4,7 —= 2,7). De waar- 
den der nieuwe basisvariabelen, die men verkrijgt door 10 voor 
te substitueren, zijn in de tweede tabel weergegeven. 





=jejejrle eel 


—4 —l —2 —6 —2 
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15. Stelling: Het is altijd mogelijk de tabel in evenwicht te bren- 
gen door hierin 6’s te plaatsen op de in 14 beschreven manier. 
Het evenwicht is uniek. 


Bewijs: Stel de oorspronkelijke basisvariabelen in de eerste tabe 
voor door «* en de oorspronkelijke basisvectoren door A*. Er 
geldt dan 


A*s* = [a 
b 


Volgens 14 voegt men nu een nieuwe kolom a» van A en een 
waarde Ò in. Bij het in evenwicht brengen van de tabel bepaalt 
men tevens nieuwe waarden van de oorspronkelijke basisvaria- 
belen. Stelt men de laatste waarden voor door «® dan geldt 


A*x0 + 0alhs) — H 
b 


x= A*-1f fa) — baths 
|b En 


= g* — OA*-L, ah 


Daar de oorspronkelijke vectoren lineair onafhankelijk zijn is 











162 De toepassing in-het bedrijf 


|A*| #0. Derhalve levert de vermenigvuldiging A*-1,a(h® een 
unieke uitkomst op. 


16. Stelling: De volgens 14 bepaalde m + n — 1 variabelen en 
vectoren zijn weer basisvariabelen en basisvectoren. 


Bewijs: Volgens 14 wordt in elke rij en kolom hetzij niets hetzij 
eenmaal +60 en eenmaal — 6 geplaatst. Als in elke rij en kolom 
van de oorspronkelijke tabel een basisvector voorkomt, zal der- 
halve ook nadat aan 0 een bepaalde waarde is toegekend in elke 
rij en kolom van de nieuwe tabel een vector voorkomen. Dit bete- 


kent, dat elke 5 in de nieuwe vectoren kan worden uitgedrukt 


(zie 8). Deze vectoren vormen dus m.a.w. een basis. 


17. Indien alle (kij + A4 + uy) van de nieuwe, volgens 14 bepaalde 
tabel groter dan of gelijk aan O zijn, is de nieuwe oplossing volgens 
13 een optimale. Is hieraan niet voldaan, dan moet uit de gevonden 
tabel weer een nieuwe volgens 14 worden bepaald. Men gaat door 
met de in 14 beschreven methode totdat een optimale oplossing 
is gevonden. 


Voorbeeld: De waarden van de variabelen in onderstaande tabel 
zijn gelijk aan die in de laatste tabel van 14. Bij deze tabel zijn de 
grootheden As, uy en de sommen (Rij + Às + uy) bepaald, uitgaande 
van Â2—=0. Het blijkt dan, dat min(Zij + A + uy) = —l in 
hokje h‚s = 3,4. In dit hokje plaatst men 9 en men brengt de 
tabel verder in evenwicht. 

Het blijkt dan, dat aan 0 de waarde 10 moet worden toegekend. 
Substitueert men deze waarde in de eerste tabel dan verkrijgt 
men de tweede. 

Ook van de tweede tabel zijn de grootheden A4, uy en de sommen 
(kig + Àt + uy) berekend, uitgaande van Â4 = 0. Het blijkt, dat 
alle (kij + As + uy) 20. De gevonden oplossing is dus een op- 
timale. 
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Het optimale transportschema kan als volgt worden weergegeven : 





A B C D E 


T __U VW _X Y Z 
20.000t. 40.000t. 410.000t. 70.000t. 10.000t. 30.000 t. 20.000 t. 
10.000 t. 30.000 t. 


De minimale transportkosten bedragen: 


20000 Xx f 1 =f 20 000 
40 000 X -2==- 80 000 
10 000 x -6 = - 60 000 
70 000 x - 4 — - 280 000 
10 000 X -3 =—- 30000 
10 000 X -3 == - 30 000 
30 000 x -6 = - 180 000 
20 000 x -3 =- 60 000 
30 000 x - 2 —= - 60 000 


f 800 000 


18. Wanneer geen degeneratie optreedt leidt de hierboven be- 
schreven oplossingsmethode in een eindig aantal stappen tot een 
oplossing. Bij elke stap verkrijgt men dan immers een basisoplos- 
sing met een hogere functiewaarde dan de voorgaande, terwijl het 
aantal basisoplossingen volgens III, $ 2, 6 eindig is. 

Degeneratie zou kunnen optreden als min(kij + A + uy) in ver- 
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schillende hokjes tegelijk voorkomt. In dat geval verkrijgt men een 
nieuwe basisoplossing waarvan een' variabele gelijk aan 0 is (of zelfs 
verschillende variabelen gelijk aan O zijn). Dit heeft geen gevolgen 
voor de verdere berekeningen als geen of geen negatieve 0 moet 
worden geplaatst in de hokjes met een O. Als een — 6 in een der- 
gelijk hokje moet worden geplaatst, is min (kij + A4 + uy) = O en 
wordt dus de O in het hokje van de oorspronkelijke tabel vervangen 
door een O in een ander hokje. De functiewaarde verandert uiter- 
aard niet door deze substitutie. In het algemeen heeft ook dit geen 
gevolgen voor de verdere berekeningen. Doorgaans doet zich in een 
volgende fase van de berekeningen een situatie voor, waarbij geen 
of geen negatieve 0 in een hokje met een O moet worden geplaatst. 
Men verkrijgt dan weer een oplossing met een kleinere functie- 
waarde dan de voorgaande. 

Theoretisch bestaat evenwel de mogelijkheid, dat de laatstgenoemde 
situatie niet optreedt en men een vroegere oplossing terugkrijgt. 
De berekeningen zijn dan gedegenereerd. 

Evenals dit in III, $5 voor de simplexmethode is aangetoond, 
kan men voor de in deze paragraaf besproken methode bewijzen, 
dat altijd minstens bij één van de uit de basis te verdrijven vec- 
toren geen degeneratie ontstaat. Er bestaat ook weer een methode 
om te bepalen, welke keuze hiertoe in elke etappe moet worden 
gedaan. Deze methode kan men afleiden door naast het gegeven 
probleem een tweede te beschouwen, waarin één van de gegeven 


Ì 
ais is Seen door Ke + Â en alle bj’ s door bj Epe — Â. De in de 


oplossing van het tweede vraagstuk Nd. À laat men 
dan tot O naderen. Hierop wordt echter niet verder ingegaan. 


s 5. Dualiteit 


1. Bij verkoop- en produktieplanning moet na de bepaling van het 
optimale‘ programma vaak worden nagegaan, welke veranderingen 
in'het verschil tussên de bruto af fabriek opbrengst en de variabele 
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fabricagekosten zouden optreden als de capaciteit der produktie- 
afdelingen of de afzetmogelijkheden zouden veranderen. De be- 
slissing tot investering is doorgaans gebaseerd op een afweging 
van het te investeren bedrag tegen de uit de capaciteitsvergroting 
voortvloeiende toeneming van het verschil tussen bruto af fabriek 
opbrengst en variabele kosten minus de toeneming van de vaste 
jaarlijkse ‘out of pocket expenses’. Als de op een bepaalde deel- 
markt af te zetten hoeveelheid bij constante gemiddelde bruto 
af fabriek opbrengst toeneemt, zal men willen nagaan welk effect 
dit op de winst heeft. 

In $3, 13 is een methode genoemd om zowel de veranderingen in 
het optimale verkoop- en produktieprogramma als die in het daar- 
bij behorende verschil tussen bruto af fabriek opbrengst en variabele 
kosten te bepalen. In het volgende wordt — zoals reeds in $ 3, 14 is 
aangekondigd — een eenvoudiger methode besproken om alleen het 
verschil tussen de bruto af fabriek opbrengst en de variabele kosten 
te berekenen. Teneinde deze methode te kunnen aangeven, moet 
in het kort worden ingegaan op het verschijnsel van dualiteit dat 
bij lineaire programmeringsvraagstukken optreedt. 


2. Stel, gegeven is het vraagstuk 
Ax <b 
« =>0 
C's = max. 


Hieruit leidt men een tweede vraagstuk af, luidende 
y A >= c’ 
y z=0 
y'b —= min. 


Het tweede vraagstuk noemt men het duale van het gegeven 
probleem. 

Op deze wijze kan uit elk maximum vraagstuk een minimum- 
vraagstuk worden afgeleid. Het omgekeerde is eveneens waar: Uit 
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elk gegeven minimum-vraagstuk kan een maximum-vraagstuk 
worden afgeleid. 

Teneinde niet in herhalingen te moeten vervallen wordt in het vol- 
gende verondersteld, dat het maximumvraagstuk is gegeven en het 
minimum-vraagstuk als duale hieruit is afgeleid. Voor deze situatie 
staat men bij de verkoop- en produktieplanning volgens $3. In 
vele andere gevallen is echter het minimumprobleem gegeven en 
het maximumprobleem hieruit afgeleid. Mutatis Mutandis geldt 
het onderstaande eveneens voor de laatste situatie. 


3. In het volgende wordt ervan uitgegaan, dat in alle voorwaarden 
Ax <b het < geldt en dus geen enkele een gelijkheid is. Dit 
betekent geen verlies van algemeenheid, aangezien een gelijkheid 
altijd door twee gelijkheids/ongelijkheidsrelaties kan worden ver- 
vangen. In plaats van 


a=b 
kan men immers schrijven 
a< b of a> b 
—as<-—b —a > —b. 


Als nu in een gegeven stelsel Ax < b een gelijkheid voorkomt, 
vervangt men door twee betrekkingen van bovenstaande vorm. 
Het oorspronkelijke en het nieuwe stelsel zijn identiek. Alles wat 
voor het tweede wordt bewezen geldt derhalve ook voor het eerste. 
Het uiterlijke verschil tussen het oorspronkelijke stelsel en het 
nieuwe is, dat het eerste een vergelijking minder heeft dan het 
tweede. Derhalve is ook het aantal variabelen in het duale van 
het oorspronkelijke vraagstuk 1 kleiner dan het aantal variabelen 
in het duale van het nieuwe vraagstuk. 


4. Stelling: De functiewaarden c'x van alle bruikbare oplossin- 
gen « van het maximumvraagstuk zijn kleiner dan of gelijk aan 
de functiewaarden y'b van alle bruikbare oplossingen y van het 
minimumvraagstuk. 
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Bewijs: Als # een bruikbare oplossing van het maximum-vraag- 
stuk is en y van het minimum-vraagstuk, dan is 
Ax < b 
yvAr=c'. 
Vermenigvuldiging van de eerste ongelijkheden met y’ en van de 
tweede met « levert op: 
“_y'Ar <y'b 
v'Ar >= C's. 
Dus: 
c's < y'b 


5. Er kan o.a. op grond van stelling 4 een tweetal stellingen tegelijk 
worden bewezen. 


Stelling: Als het oorspronkelijke vraagstuk een eindige oplossing 
heeft, heeft het duale ook een eindige oplossing en 


min y'b — max C'x. 


Stelling: De optimale oplossing y’ van het duale vraagstuk is 
gelijk aan de waarden van de verschilvariabelen in het oorspron- 
kelijke vraagstuk. 


Bewijs: De herleide vorm van het maximum-vraagstuk volgens 
III, 8 1 kan worden voorgesteld door 


BE 


zj 
OT 


Stel, dat een optimale oplossing ne | rj >= 0 van dit herleide vraag- 


stuk is gevonden door toepassing van de simplexmethode. De ‘als. 











| 
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laatste gevonden tabel van coëfficiënten wij wordt voorgesteld 
door de matrix van de orde m Xx (n + m) waarin 


Ui = (uu: Hin U2= [ina 07007 Hiin+m) 


Hm1**"*" Hmn Hm(n+1) leken Hn(n+m) 


De kolommen van Ui stellen dus voor: de gewichten in de lineaire 
* 
combinaties der bij pl behorende basisvectoren, welke als 


resultaat de kolommen van A opleveren. De kolommen van U2 
stellen voor: de gewichten in de lineaire combinaties dezer basis- 
vectoren, welke als resultaat de eenheidsvectoren in I opleveren. 
Stelt men de mm basisvectoren voor door A*, dan is dus 


[AI] = A*{U, Us} 
zodat 
A= A*U, 
en I= A*U, dus Uz = A*1, 
De verschillen (c‚ — 24) behorende bij de optimale oplossing zijn 
alle niet-positief. Dus: 


(or == ij eenen ‚(en — Zn), (O — Zn) 117" ‚(0 — Zintm)) SO. 
Stel 
v, = [en — zi) «ere (en — zn)l SO 
vz = —[(O — Br) 279 (O — Zenm))] 
= [4Znai) "rr FZintmjl = 0. 


Stel voorts de m coëfficiënten van c,‚ die bij de basisvariabelen 
behoren, voor door c*. 
Daar wv,=c —c*U; en dus c*'U; =c —v 
en daar vs —= C*Uz = c*'A*-1 
is vaA —= c*'A*-1, A*U, 
== c*’U: 


== —Uv EE. 
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Het blijkt dus, dat và een bruikbare oplossing van het duale mini- 
mum-vraagstuk is. 
Voorts is 
vsb —= Cc*'A*L, A*y* 
= e*'x*, 
Volgens 4 zijn de functiewaarden y'b van het minimum-vraagstuk 
voor alle bruikbare oplossingen y groter dan of gelijk aan de 


functiewaarden c'x van het maximum-vraagstuk voor alle bruik- 
bare oplossingen «. Dit geldt dus ook voor #*; zodat: 


yb <ct'x* voor alle 7’. 
Dus: 
vb <vab voor alle y’. 


Hiermede is gezegd, dat y’ een optimale oplossing van het mini- 
mum-vraagstuk is, terwijl min y'b —= vab — c*'x* — max C'%. 


Voorbeelden: Uit het in III, $ 4, 3 gegeven vraagstuk kan het vol- 
gende duale probleem worden afgeleid: 


yi + 2ya + 33 >= | 
Wi Hat Y3 2 | 
yi + Ya + 293 ZI 
Yu Ya, Va =O 
6y1 + 6y2 + 6y3 = min. 
In III, $ 4, 20 is het oorspronkelijke vraagstuk opgelost. Uit de 
laatste tabel blijkt, dat de oplossing van het duale vraagstuk luidt: 


Vi= ta 5 


y= == tE 
Vz = 6 = 5 


In III, $ 4, 20 werd als maximum van de functiewaarde 3 gevonden. 
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Het minimum van het duale vraagstuk blijkt eveneens gelijk aan 3 
te zijn: 
68 + GB +64 =3 


Als tweede voorbeeld kan het in $3, 10 genoemde planningmodel 
worden beschouwd. Hieruit kan het volgende duale vraagstuk 
worden afgeleid: 


0,38y1 + 0,6y3 + ya + ya > 138 
0,381 + 0,6y3 + ya + y8 > 158 
0,181 + 0,8y3 +5 + y9 = 128 
0,31 + 0,8y2 +6 +10 > 108 
0,3y1 + 0,8y2 +Y6 +yu 2 88 


75 000y; + 80 000y2 + 85 000yz +- 150 OOOya + 60 00Oys + 
+ 120 000yg + 100 000yz + 80 O0Oyg + 50 000yg + 
+ 80 00Oy10 + 60 00Oy11 = min 

Yu Ya, Y3, Va, Y5, Ve, V7, Ye, Vo, Y10, Vui =O. 


Uit de oplossing van het maximum-vraagstuk kan de volgende 
oplossing van het minimum-vraagstuk worden bepaald: 


Y1 == = 50 
ya =247 = 90 
y3 == 28 == 200 
Ya = 29 
Y5 =Z0o= OO 
Ye =Zu= O 
yr == O 
Ys —Z1i38 = 20 
Yo =zia= O 
Y1o = 215 — 20 


Il 
o 


Vui = 21e = O 
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Het minimum van het duale vraagstuk is weer gelijk aan het 
maximum van het oorspronkelijke: 


75 000 x 50 == 3750 000 
80 000 Xx 9 = 7 200 000 
85 000 x 200 —= 17 000 OOO 
80 000 Xx 20 == 1 600 000 
80 000 Xx 20 == 1 600 000 


31 150 000 


6. De in 1 aangekondigde methode ter bepaling van een winst- 
stijging door verandering van een produktiecapaciteit kan nu ge- 
makkelijk worden aangegeven : 

Beschouw het in $3 besproken planningmodel en het daarbij be- 
horende duale minimum-vraagstuk. Voor het gemak worden beide 
nog eens in hun algemene vorm naast elkaar weergegeven. 


Ax <b VAC 
«=0 y =0 
c'x = max y'b = min 


Als «* een optimaal verkoop-en produktieprogramma is en y*’ een 
optimale oplossing van het minimum-vraagstuk, dan is volgens 5 


c's* = y*'D. 


Zoals in $ 3, 13 is opgemerkt, heeft een verandering van b in het 
maximum-vraagstuk een verandering van de optimale oplossing tot 
gevolg. Blijkens $ 3, 15 leidt een verandering van b in het mini- 
mum-vraagstuk binnen zekere grenzen evenwel niet tot een ver- 
andering van y*. 

Stel nu, dat een element bj van b toeneemt met Ab, bijv. door in- 
vestering in een produktieafdeling. In dat geval zal het minimum 
van het duale vraagstuk volgens $3,15 toenemen (mits alle 
(cy — 24) =O blijven) met 


Abiyt. 
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Aangezien echter max c'x —= min y’b, zal ook het maximum van 
het oorspronkelijke planningmodel toenemen met 4b;y;. 

Dus: De verandering in het verschil tussen de bruto af fabriek op- 
brengst en de variabele fabricagekosten door toeneming van een 
element van b in het planningmodel is gelijk aan het produkt van 
deze verandering en de waarde van de overeenkomstige variabele 
in het duale probleem. De waarden van de variabelen in het duale 
vraagstuk geven m.a.w. aan, welke verandering optreedt in het ver- 
schil tussen de bruto af fabriek opbrengst en de variabele kosten. 
als één van de elementen van b met een eenheid verandert. Om deze 
reden kunnen de elementen van y*’ interne rekenprijzen worden 
genoemd. 


Voorbeelden: Als de capaciteit van de ammoniakfabriek in het in 
$3, 10 genoemde voorbeeld met 1000 t/jaar zou kunnen worden 
vergroot, zou het verschil tussen de bruto af fabriek opbrengst 
volgens de in 5 van deze paragraaf berekende yj toenemen met 
1000 x f 50 = f 50 000/jaar. Een vergroting van de zwavelzuur- 
fabriek of de salpeterzuurfabriek met 1000 ton zou dit verschil 
toenemen met 1000 x f 90 —= f 90 OOO/jaar, resp. 1000 x f 200 = 
== f 200 000O/jaar. Een vergroting van de capaciteit van de KAS- 
fabriek, de KS-fabriek of de ZA-fabriek zou geen enkele verande- 
ring in het verschil tussen bruto af fabriek opbrengst en variabele 
kosten teweeg brengen. Dit is ook te verwachten, daar elk van deze 
fabrieken onderbezet is. 

Een vergroting van de afzetmogelijkheid van KAS en KS in het 
bouwland en van ZA in het grasland heeft geen gevolgen voor de 
winst. Ook dit is begrijpelijk, daar in deze markten nog niet het 
maximum wordt verkocht. Zou daarentegen 1000 t KAS meer in 
het grasland of 1000 t ZA meer in het bouwland kunnen worden 
afgezet, dan zou de winst met 1000 x f 20 = f 20 000/jaar toe- 
nemen. 

De berekening van bovengenoemde vergrotingen van de ‘cash flow’ 
ter grootte van f 50 000, f 90 000 en f 200 000 kan het uitgangs- 
punt vormen voor een investeringsanalyse. Door deze bedragen te 
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plaatsen tegenover de investering, vereist om de betrokken capaci- 
teitsvergroting te verwezenlijken, kan een rentabiliteitsberekening 
worden gemaakt. Op grond van een onderzoek naar de financierings- 
mogelijkheden kan voorts worden beoordeeld of in feite tot inves- 
tering zal worden overgegaan. 

Het laatstgenoemde bedrag van f 20 000 zou de aanleiding voor 
het uitvoeren van een marktstudie kunnen zijn. Er zou kunnen 
worden onderzocht op welke wijze de KAS-afzet in het grasland 
of de ZA-afzet in het bouwland met 1000 t zou kunnen worden 
vergroot. Men zou kunnen overwegen de verkooppolitiek op een 
dergelijke afzetvergroting te richten indien de kosten hiervan lager 
dan f 20 000 zouden zijn. 
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